Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,436,956 2 bp→TT coding (202‑203/741 nt) cysQ → 3'(2'),5'‑bisphosphate nucleotidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,436,9560CT93.1% 195.4 / 4.6 73P68F (CCC→TTC) cysQ3'(2'),5'‑bisphosphate nucleotidase
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base T (25/42);  total (29/44)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-02
*NC_0009134,436,9570CT94.3% 175.5 / 4.2 73P68F (CCC→TTC) cysQ3'(2'),5'‑bisphosphate nucleotidase
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base T (25/41);  total (29/44)
Fisher's exact test for biased strand distribution p-value = 2.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-02

AAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTCCGCGATG  >  NC_000913/4436721‑4437193
                                                                                                                                                                                                                                           ||                                                                                                                                                                                                                                            
aaatgaaccccGCTCACAGAGACGAAGTGGAGAAATGTTAGATCCAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGGTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTc                                                                                                                                                                                                                                 <  2:100046/241‑1 (MQ=255)
 aGTTAACACCGCTCACCGAGACGCGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGCATGCAGGCGATGCCATTATGCAGGTCTCCGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGCCAAATCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCGTGGGCGGTTTAATGACGCTGAAAACGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAGGTcc                                                                                                                                                                                                                                >  2:262835/1‑250 (MQ=255)
   ttAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCCCACCGTTATCATGGACGGTTTACGTCCGCTGAAACCGGATGTTCCTGTCCTTTCTGAATGAGATCCTTTCGGTTGGGAAGTCCGTc                                                                                                                                                                                                                             >  2:143973/1‑251 (MQ=255)
   ttAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTc                                                                                                                                                                                                                             >  1:119673/1‑251 (MQ=255)
              cacaGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGc                                                                                                                                                                                                                    >  2:45199/1‑249 (MQ=255)
                 agagACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGt                                                                                                                                                                                                               >  2:172979/1‑251 (MQ=255)
                  gagaCGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCATGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTGGGGAAGTCCGTCAGCACTGGCAGCGtt                                                                                                                                                                                                              <  2:317057/251‑1 (MQ=255)
                           tGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAAATCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTActggctgg                                                                                                                                                                                                     >  1:254943/1‑251 (MQ=255)
                            ggAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCGCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTActggctgg                                                                                                                                                                                                     <  2:313575/250‑1 (MQ=255)
                             gagaAATGTTAGATCAAGTATTCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTACGTCAGCACTGGCAGCGTTACTGGCTGGTa                                                                                                                                                                                                   <  1:514772/251‑1 (MQ=255)
                               gaAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGa                                                                                                                                                                                                 <  1:293274/251‑1 (MQ=255)
                                aaaTGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGAc                                                                                                                                                                                                <  2:308861/251‑1 (MQ=255)
                                 aaggTTAGATCAAGTATGGCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGAGCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGAcc                                                                                                                                                                                               <  1:12054/248‑1 (MQ=255)
                                 aaTGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCGGCACTGGCAGCGGTACTGGCTGGGAGAcc                                                                                                                                                                                               >  2:166800/1‑251 (MQ=255)
                                  aTGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGAAGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGTGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGAccc                                                                                                                                                                                              <  2:517660/251‑1 (MQ=255)
                                     ttAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCt                                                                                                                                                                                           >  1:237242/1‑251 (MQ=255)
                                     ttAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACATACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGTCCCGCt                                                                                                                                                                                           >  1:161758/1‑251 (MQ=255)
                                      tAGATCAAGTTAGCCAGCTTGCACTGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGGTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCt                                                                                                                                                                                           <  2:159304/250‑1 (MQ=255)
                                      tAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTACTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCt                                                                                                                                                                                           <  1:131656/250‑1 (MQ=255)
                                         aacAAGGATGCCAGCTTGCACGGAAGGCCGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGGCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGAt                                                                                                                                                                                       <  2:336612/249‑1 (MQ=255)
                                                aTGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACt                                                                                                                                                                                 <  1:506200/250‑1 (MQ=255)
                                                  gCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGGAGACCCGCTGggtggtcctta                                                                                                                                                                               >  2:104880/1‑245 (MQ=255)
                                                   ccAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTaa                                                                                                                                                                               >  1:457499/1‑249 (MQ=255)
                                                              ggaagtccggCGATGCCCGTATGCAGGTCTACGTCCGGACGAAACCGATGGAGGTCGGCAGCAAAGCGGCAAATTCTCCGGTAACGGCATCGGATTATGCCGCTCACACCGTTACAATGGTCGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAACTCCGTCAGCACAGGCAGCGTTACTGGCTGGTAGTCCCGCTGTATGGTACTAAAGAGATTATTaaa                                                                                                                                                                  <  2:535523/243‑1 (MQ=255)
                                                              ggAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTaaa                                                                                                                                                                  <  1:45199/251‑1 (MQ=255)
                                                              ggAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTaa                                                                                                                                                                   >  1:126155/1‑250 (MQ=255)
                                                                       gCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAACGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGc                                                                                                                                                         <  2:120455/251‑1 (MQ=255)
                                                                        cGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGGCCGTCAGCACTGGCCGCGTTACTGGCTGGGTGACCCGCTGGGTGGTACTAAAGAGTTTATTAAACGTAATGGCg                                                                                                                                                        >  2:209319/1‑251 (MQ=255)
                                                                           tGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGCCGTCGTCAGCTAAGCGGACACTTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATCGACGGTTTACGTCCGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAATATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTTGGCCCACTGGTTGGACCTTAAGGGTTTTTTATACGTAGTGGCGAAt                                                                                                                                                     >  2:86050/1‑251 (MQ=255)
                                                                             ccATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGGTGGTACTAAAGAGTTTATTAACCGACATGGCGAATTc                                                                                                                                                   >  1:286622/1‑251 (MQ=255)
                                                                              cATTATGCTGGACTACGACGGGACGAAACCGATTGACGTCGTCAGCAAAGCGTACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGGTTACTTACGCTGACACCGGATGTTCCTGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAATACCGTCAGCACTGGCAGCGTTACTTGCTGGTAGACCCGCTGGATGGTTCTAAAGAGTTTATTAAACGTAATGGCGAATTCa                                                                                                                                                  <  1:348377/251‑1 (MQ=255)
                                                                              cATTATGCAGGTCTACGTCGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACCCCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCa                                                                                                                                                  <  2:184155/251‑1 (MQ=255)
                                                                                    gCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACTCCGGATGTTCCGGTCATTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGtt                                                                                                                                             <  1:370370/250‑1 (MQ=255)
                                                                                      gggTCTACGACGGGACGAAACCGATGGACGTTGTCAGCAAAGCGGACAATTCTCCGGTAACTGCAGCGGATATTGCCGCTCACACCGTTCTCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTGTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGGAATGGCGAATTCACCGTTAAc                                                                                                                                          <  2:43596/250‑1 (MQ=255)
                                                                                        gTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACAt                                                                                                                                        <  1:404616/251‑1 (MQ=255)
                                                                                         gcTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCGCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCGTTCGGTTGGGAAGTCCTTCAGCACTTGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTGATTAAACGTAATGGCGAATTCACCGTTAACAtt                                                                                                                                       <  2:464974/250‑1 (MQ=255)
                                                                                              gacgGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGCCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGGTGGGACTAAAGAGTTTTTTAAACGTAAAGGCGGATTTACCGTTTACAtttatgt                                                                                                                                  >  2:8737/1‑246 (MQ=255)
                                                                                                                gACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAAcc                                                                                                                >  1:439319/1‑251 (MQ=255)
                                                                                                                        cAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCAGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTa                                                                                                         <  1:417207/250‑1 (MQ=255)
                                                                                                                            aaaGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGCGTTTATTAAACGTAATGGCGAATTCACCGGTAACATTGCGCTCATTGAACATGGCaaa                                                                                                                  >  2:146821/1‑237 (MQ=255)
                                                                                                                            aaaGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCaaa                                                                                                                  <  1:146821/237‑1 (MQ=255)
                                                                                                                            aaaGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGATTATTAAAAGTCATGGCGAATTACCCGTTAACATTGAGCTCATTGACCATGGCAAAACGGTTTTAGGCgt                                                                                                    >  1:185396/1‑251 (MQ=255)
                                                                                                                             aaGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCgtg                                                                                                   <  1:172979/251‑1 (MQ=255)
                                                                                                                              agaggACAATGCGCCGGGAACGGCAGCGGATATTGCCGCTCACACCGGTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGAGCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATGGACCATGGCAAACCGAATTTAGGCgtgg                                                                                                  <  2:196800/248‑1 (MQ=255)
                                                                                                                               gCGGACAATTCTCCGGGAACGGCAGCGGATATTGCCGCTCACACCGTTATCAGTGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCgtggt                                                                                                 <  2:220026/251‑1 (MQ=255)
                                                                                                                                cGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGGTGGTACTAAAGAGTTTATTAAACGTAAGGGCGACTTCACCGTTAACATTGCGCTCATTGACCATGTCAAACCGATTTTAGGCgtgggg                                                                                                >  1:151442/1‑249 (MQ=255)
                                                                                                                                  gACAATTCTCCTGTAACGGCAGCGGATATTGCCGCTCACTCCCTAATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCGTTCTGAAGAAGAGCCTTTCGGTGGGGAAGTCAGTCAGCACTGGCAGCGTTACTGGATGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTa                                                                                              <  2:205411/251‑1 (MQ=255)
                                                                                                                                  gACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTa                                                                                              <  1:4962/251‑1 (MQ=255)
                                                                                                                                  gACAATTCTCCGGTAACGGCAGCGGAAATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTACGCCCGGTGGCTGGTACTAAAGAATTTATTAAACGTAATGGCGACATCACCGTTAACATTTCCATCTATGACCATGGCAAACAAATTTTaagagtggggt                                                                                               >  2:340007/1‑240 (MQ=255)
                                                                                                                                           ccGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCg                                                                                       <  1:145414/249‑1 (MQ=255)
                                                                                                                                           ccGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCg                                                                                       <  2:506966/249‑1 (MQ=255)
                                                                                                                                            cGGTAACGGCTGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGAAACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATTCTCCGGTa                                                                                    <  2:440262/251‑1 (MQ=255)
                                                                                                                                                 aCGGCAGCGGATATTGCCGCTCACTCCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTGCGGTTTGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGaa                                                                               <  1:542965/251‑1 (MQ=255)
                                                                                                                                                  cGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAAc                                                                              <  1:351178/251‑1 (MQ=255)
                                                                                                                                                   ggCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACg                                                                             <  1:220181/251‑1 (MQ=255)
                                                                                                                                                      aGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTa                                                                           >  1:355112/1‑250 (MQ=255)
                                                                                                                                                          gATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTGCGGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATGGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTa                                                                      <  2:477125/251‑1 (MQ=255)
                                                                                                                                                                      cacaCCTTTCTCATGGACGGTTTACGTACGCTGACACCGGAAGTTCCGGTCCTTTCTGAAGAAGATCCTTGCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGGCCCGCTGGATGGTACTAAATAGTTTATTAAACGTAATGGCGAAGTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTGTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCa                                                            <  1:240635/249‑1 (MQ=255)
                                                                                                                                                                            gTTATCATGGACGGTTTACGTACGCTGCCACCGGATGTTCCGGTCCTTTATGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGAAACGCTGGATGGTACTAAAGAGTTTATTAAACGGAATGGCGAATTCCCCGTTAACATTTGGCTCATTTGCCATAGCAAACCGATTTTAGGCGTGGGGTATGCCGCGGTAATGAACGGAATGTACacccggggcgaaggca                                                     >  2:233719/1‑235 (MQ=255)
                                                                                                                                                                                 cATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGATCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGt                                               <  1:211632/251‑1 (MQ=255)
                                                                                                                                                                                  ctGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGTGCTCATTGACCATGGCAAAACGATTTTAGGCGTGGTGGATGCGCCGGGAATTAACGTAATGAACAGCGCGGCAGACGGCAAAGCggg                                              >  2:113256/2‑249 (MQ=255)
                                                                                                                                                                                          tttACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGACGGCAAAGCGTGGCaaga                                        >  1:397633/1‑249 (MQ=255)
                                                                                                                                                                                                  cgCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGgt                              <  1:104880/251‑1 (MQ=255)
                                                                                                                                                                                                          ccGGATGTTCCGGGCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCaa                      <  2:126155/251‑1 (MQ=255)
                                                                                                                                                                                                            gtctGTTCCGGTCCTTTCTGAATAAGATCCTCTCGTTTGGGCAGGCCGTCAGCAATGTCAGCGTTAATGGTTGGTAGACCCGCTTGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTATCATTGCACACATGGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGTCAGAAGGCAAAGCGTGGAACGAAGAGTGCGGTGTgcgccagg                    <  2:397866/248‑5 (MQ=255)
                                                                                                                                                                                                                gTTCCGGTCCTTTCTTAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGc                    >  1:6939/1‑247 (MQ=255)
                                                                                                                                                                                                                  tcaggTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGGTTTTAGGCGTGGTGTATGCGCCGGTAATGCACGTAATGTACAGCGCGGCGGAAGGCCAAGCGGGGAAAGAAGAGTGCGGTGTGCGCAAACAGATTc              >  2:237284/4‑251 (MQ=255)
                                                                                                                                                                                                                   ccGGTCCTTTCTGAAGAAGATCCTTGCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTc              <  1:295958/250‑1 (MQ=255)
                                                                                                                                                                                                                    cGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCCGAAGGCAAAGCGGGGGAAGAAGAGTGCGGGGTGCGCAAGCAGGTTc              >  2:541041/1‑249 (MQ=255)
                                                                                                                                                                                                                     ggTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAgg           <  1:504752/251‑1 (MQ=255)
                                                                                                                                                                                                                           ttcggaataagaTCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTTGTAGACCCGCTGGATGGTACTAAAGAGTTGAGTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCTCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAATCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTCcgc      <  2:119673/242‑1 (MQ=255)
                                                                                                                                                                                                                              tGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTCCGCGATg  >  1:90280/1‑251 (MQ=255)
                                                                                                                                                                                                                              ggaagaaGATCCTTTCGGGTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCATAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTCCGCGATg  <  1:38406/250‑1 (MQ=255)
                                                                                                                                                                                                                                           ||                                                                                                                                                                                                                                            
AAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTCCGCGATG  >  NC_000913/4436721‑4437193

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑