Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 378,400 G→A 5.9% intergenic (‑32/+62) frmB ← / ← frmA S‑formylglutathione hydrolase/alcohol dehydrogenase class III; glutathione‑dependent formaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913378,4000GA 5.9% 606.1 / 15.6 185intergenic (‑32/+62)frmB/frmAS‑formylglutathione hydrolase/alcohol dehydrogenase class III; glutathione‑dependent formaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (78/96);  new base A (6/5);  total (84/101)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

CGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGT  >  NC_000913/378252‑378548
                                                                                                                                                    |                                                                                                                                                    
cgcgGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGTCATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCACCGGACg                                                                                                                                                      >  3:1485409/1‑149 (MQ=255)
           gtggGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGGCAGCCGTCCggtg                                                                                                                                           <  2:1280893/149‑3 (MQ=255)
            tggGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGga                                                                                                                                          <  1:853041/149‑1 (MQ=255)
               gagaTAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGagag                                                                                                                                       >  2:642265/1‑149 (MQ=255)
               gagaTAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGagag                                                                                                                                       <  1:416118/149‑1 (MQ=255)
                 gaTAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTa                                                                                                                                     >  3:1112937/1‑149 (MQ=255)
                   tAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTAtt                                                                                                                                   <  3:79757/149‑1 (MQ=255)
                      aCGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAg                                                                                                                                <  3:630059/149‑1 (MQ=255)
                           gACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAAccc                                                                                                                           >  1:1018676/1‑149 (MQ=255)
                              attcattTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCg                                                                                                                        >  3:466108/1‑149 (MQ=255)
                              attcattTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCg                                                                                                                        >  3:1297626/1‑149 (MQ=255)
                               ttcattTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCggg                                                                                                                                            <  3:940291/128‑1 (MQ=255)
                               ttcattTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCggg                                                                                                                                            >  4:940291/1‑128 (MQ=255)
                                tcattTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGcc                                                                                                                      <  3:1064523/149‑1 (MQ=255)
                                  atttcacatttcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGGCAGCCGTCCgggg                                                                                                                                           >  1:29587/1‑126 (MQ=255)
                                  atttcacatttcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGGCAGCCGTCCgggg                                                                                                                                           <  2:29587/126‑1 (MQ=255)
                                   tttcacatttcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGGCAGCCGTCCggggga                                                                                                                                         <  1:282675/127‑3 (MQ=255)
                                   tttcacatttcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGa                                                                                                                   >  3:549726/1‑149 (MQ=255)
                                     tcacatttcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGGCAGCCGTCCgggg                                                                                                                                           <  3:1450864/123‑1 (MQ=255)
                                             caGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACggg                                                                                                         <  3:1438527/149‑1 (MQ=255)
                                             caGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACggg                                                                                                         <  4:454274/149‑1 (MQ=255)
                                              ggTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACgggg                                                                                                        <  2:599284/148‑1 (MQ=255)
                                                  aTTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTaaa                                                                                                    <  4:1526204/149‑1 (MQ=255)
                                                   ttGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAAc                                                                                                   <  2:69558/149‑1 (MQ=255)
                                                     gggAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACgggg                                                                                                        <  1:1622599/142‑1 (MQ=255)
                                                     gggAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACgggg                                                                                                        >  2:1622599/1‑142 (MQ=255)
                                                       gAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGc                                                                                               <  2:1398461/149‑1 (MQ=255)
                                                          taatGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTa                                                                                                                                  <  3:960193/111‑1 (MQ=255)
                                                                cGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCa                                                                                      <  4:1613148/149‑1 (MQ=255)
                                                                  aTACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGt                                                                                    >  2:476964/1‑149 (MQ=255)
                                                                  aTACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGt                                                                                    >  2:1557471/1‑149 (MQ=255)
                                                                    aCATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTaa                                                                                  <  2:579200/149‑1 (MQ=255)
                                                                      atatTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACg                                                                                >  2:1445097/1‑149 (MQ=255)
                                                                       tatTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGc                                                                                               <  4:502344/133‑1 (MQ=255)
                                                                        atTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTTAACCTGCGGGAAATCAGTATCGaa                                                                              <  1:1491306/149‑1 (MQ=255)
                                                                          tttGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAAtt                                                                            <  3:369991/149‑1 (MQ=255)
                                                                           ttGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTa                                                                           <  3:1284783/149‑1 (MQ=255)
                                                                               cAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGt                                                                       >  2:857849/1‑149 (MQ=255)
                                                                                aGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTAt                                                                                                                                    <  2:1008167/87‑1 (MQ=255)
                                                                                aGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTAt                                                                                                                                    >  1:1008167/1‑87 (MQ=255)
                                                                                aGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCg                                                                                              >  1:1467820/1‑125 (MQ=255)
                                                                                aGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCg                                                                                              <  2:1467820/125‑1 (MQ=255)
                                                                                 gccgccAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTc                                                                     >  3:776299/1‑149 (MQ=255)
                                                                                    gccAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGaa                                                                  <  2:1235589/149‑1 (MQ=255)
                                                                                     ccAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCATTAACGAATTACGGTTCGAAt                                                                 <  2:114015/149‑1 (MQ=255)
                                                                                     ccAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCCCAGAGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAAt                                                                 <  4:1439008/149‑1 (MQ=255)
                                                                                      cAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATg                                                                >  3:791155/1‑149 (MQ=255)
                                                                                         aGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGAt                                                             <  2:1310053/149‑1 (MQ=255)
                                                                                            tGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCgg                                                                                             >  2:655899/1‑114 (MQ=255)
                                                                                            tGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCgg                                                                                             <  1:655899/114‑1 (MQ=255)
                                                                                             gACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGc                                                         <  4:1236952/149‑1 (MQ=255)
                                                                                              aCATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGcc                                                        >  4:1260899/1‑149 (MQ=255)
                                                                                              aCATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGcc                                                        >  2:1052297/1‑149 (MQ=255)
                                                                                                aTGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCtt                                                      >  4:449102/1‑149 (MQ=255)
                                                                                                 tGTTTTACAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTc                                                     >  4:1449993/1‑149 (MQ=255)
                                                                                                  gTTTTTCAATGAGATCCATCAGTCATTAGCTCAGGGTTAAGCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCttta                                                    <  1:1527025/149‑3 (MQ=255)
                                                                                                   tttttCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGTAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCAt                                                   <  2:541069/149‑1 (MQ=255)
                                                                                                   tttttCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCAt                                                   >  1:541069/1‑149 (MQ=255)
                                                                                                        cAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATc                                              <  4:783541/149‑1 (MQ=255)
                                                                                                        cAATGAGTTCCATCAGTCATTATCTCAGCGTTAGTCACCCGGACGGCTGCCGGGGGAGAGTATTAAGAACCCCCGCCTGAAGAGGACGGGGGAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTTCCTTCATGCATc                                              >  4:1317733/1‑149 (MQ=255)
                                                                                                          aTGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACgggg                                                                                                        >  1:1121985/1‑89 (MQ=255)
                                                                                                          aTGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAg                                            >  3:1467581/1‑149 (MQ=255)
                                                                                                          aTGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAg                                            >  4:258278/1‑149 (MQ=255)
                                                                                                           tGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAgg                                           <  1:1718231/149‑1 (MQ=255)
                                                                                                                tCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGaa                                      >  2:111056/1‑149 (MQ=255)
                                                                                                                tCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGaa                                      <  4:1398636/149‑1 (MQ=255)
                                                                                                                tCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGaa                                      <  3:690697/149‑1 (MQ=255)
                                                                                                                   aTCAGTCATTATCTCAGGGTTAGTCCCCCGGACGGCTGCGCGGGGAGAGTATTAAGAACCCGCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGc                                   <  4:974889/149‑1 (MQ=255)
                                                                                                                       gTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCa                               <  2:994677/149‑1 (MQ=255)
                                                                                                                       gTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCa                               >  3:170255/1‑149 (MQ=255)
                                                                                                                       gTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGATAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCa                               >  4:1128909/1‑149 (MQ=255)
                                                                                                                        tCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCAt                              >  2:1546203/1‑149 (MQ=255)
                                                                                                                        tCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCAt                              >  2:171201/1‑149 (MQ=255)
                                                                                                                         cATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCAtt                             >  1:1526984/1‑149 (MQ=255)
                                                                                                                             aTCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAAtt                         <  4:1041213/149‑1 (MQ=255)
                                                                                                                             aTCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAAtt                         <  4:211333/149‑1 (MQ=255)
                                                                                                                                tcAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCa                                                                                      <  1:608156/85‑1 (MQ=255)
                                                                                                                                tcAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGt                                                                       >  1:1077489/1‑100 (MQ=255)
                                                                                                                                tcAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGc                                   <  2:193115/136‑1 (MQ=255)
                                                                                                                                tcAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGc                                   >  1:193115/1‑136 (MQ=255)
                                                                                                                                tcAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCa                      >  2:1654405/1‑149 (MQ=255)
                                                                                                                                 cAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCAt                     <  4:171872/149‑1 (MQ=255)
                                                                                                                                   gggTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGAttt                                                           >  4:663579/1‑109 (MQ=255)
                                                                                                                                     gTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAAt                                                                                        <  1:343394/78‑1 (MQ=255)
                                                                                                                                     gTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAg                 <  1:67075/149‑1 (MQ=255)
                                                                                                                                          tCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCa            >  3:1574776/1‑149 (MQ=255)
                                                                                                                                           tccccGGGCAGCCGTCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTACTTTCATCCAGGCTCAt           >  3:66076/2‑149 (MQ=255)
                                                                                                                                           cccccGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCAt           <  1:79944/149‑1 (MQ=255)
                                                                                                                                           cccccGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCAt           >  1:1308908/1‑149 (MQ=255)
                                                                                                                                            ccccGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATg          <  4:549726/149‑1 (MQ=255)
                                                                                                                                            ccccGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATg          <  1:1654405/149‑1 (MQ=255)
                                                                                                                                            ccccGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATg          >  3:15260/1‑149 (MQ=255)
                                                                                                                                            ccccGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATg          <  3:258278/149‑1 (MQ=255)
                                                                                                                                             cccGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGtt                                                                      >  2:487674/1‑88 (MQ=255)
                                                                                                                                             cccGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGtt                                                                      <  1:487674/88‑1 (MQ=255)
                                                                                                                                             cccGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGAGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATgg         >  3:33189/1‑149 (MQ=255)
                                                                                                                                              ccGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGt        >  2:216352/1‑149 (MQ=255)
                                                                                                                                               cGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTa       <  3:181125/149‑1 (MQ=255)
                                                                                                                                                 gACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATg     >  1:906137/1‑149 (MQ=255)
                                                                                                                                                    ggCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGt  <  4:1456434/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGT  >  NC_000913/378252‑378548

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: