Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,627,704 A→C 5.6% E63A (GAG→GCG)  ydfG → 3‑hydroxy acid dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,627,7040AC 5.6% 226.6 / 5.0 72E63A (GAG→GCG) ydfG3‑hydroxy acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (28/40);  new base C (2/2);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

TGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGTATATA  >  NC_000913/1627555‑1627853
                                                                                                                                                     |                                                                                                                                                     
tGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAaga                                                                                                                                                       <  2:1332180/150‑1 (MQ=255)
 gAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAagag                                                                                                                                                      <  2:943208/150‑1 (MQ=255)
  aaTGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATGATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAgaga                                                                                                                                                     >  1:1003083/1‑150 (MQ=255)
  aaTGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAgaga                                                                                                                                                     >  1:282598/1‑150 (MQ=255)
  aaTGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAgaga                                                                                                                                                     >  1:147536/1‑150 (MQ=255)
      cATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGCGATGcc                                                                                                                                                 >  1:1198299/1‑149 (MQ=255)
      cATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCt                                                                                                                                                 >  1:1012239/1‑150 (MQ=255)
         tACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGc                                                                                                                                              >  2:1058459/1‑150 (MQ=255)
            tcgtcgTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATc                                                                                                                                           <  1:1386419/150‑1 (MQ=255)
             cgtcgtTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCg                                                                                                                                          >  2:7332/1‑150 (MQ=255)
              gtcgtTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGc                                                                                                                                         <  1:106274/150‑1 (MQ=255)
                 gtTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGCGATGCCAGCATCTCTTc                                                                                                                                      <  2:1198299/150‑1 (MQ=255)
                 gtTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTc                                                                                                                                      >  1:1439026/1‑150 (MQ=255)
                      aTTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGcc                                                                                                                                 >  2:1086574/1‑150 (MQ=255)
                       ttCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGCGATGCCAGCATCTCTTcaatagc                                                                                                                                <  1:1523608/150‑7 (MQ=255)
                       ttCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCg                                                                                                                                <  2:1079954/150‑1 (MQ=255)
                          aacaaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgt                                                                                                                             >  1:646757/1‑150 (MQ=255)
                           acaaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAAGTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtg                                                                                                                            <  2:1039126/150‑1 (MQ=255)
                           acaaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtg                                                                                                                            <  1:470066/150‑1 (MQ=255)
                            caaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtgg                                                                                                                           <  2:370375/150‑1 (MQ=255)
                            caaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtgg                                                                                                                           <  2:24167/150‑1 (MQ=255)
                                  cATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAat                                                                                                                     >  2:1515199/1‑150 (MQ=255)
                                     aaaGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACTAACTGGGAGATAATCTGTATATCGCCCAACTTGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATAtt                                                                                                                  <  2:336178/150‑1 (MQ=255)
                                      aaGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTg                                                                                                                 <  1:866392/150‑1 (MQ=255)
                                      aaGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTg                                                                                                                 <  2:1012239/150‑1 (MQ=255)
                                       aGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGa                                                                                                                <  1:280879/150‑1 (MQ=255)
                                       aGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGa                                                                                                                <  2:147536/150‑1 (MQ=255)
                                             cGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCt                                                                                                          <  2:829864/150‑1 (MQ=255)
                                                 aCTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTa                                                                                                      <  1:367123/150‑1 (MQ=255)
                                                           gCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATGATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCg                                                                                            <  2:1003083/150‑1 (MQ=255)
                                                           gCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCg                                                                                            >  1:1073968/1‑150 (MQ=255)
                                                               ggAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCt                                                                                        <  1:983441/150‑1 (MQ=255)
                                                                 aaCGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTgg                                                                                      >  1:1161260/1‑150 (MQ=255)
                                                                  aCGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGc                                                                                     <  2:59802/150‑1 (MQ=255)
                                                                   cGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCg                                                                                    <  1:23119/150‑1 (MQ=255)
                                                                    gggtGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAAGCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGt                                                                                   <  2:1360038/147‑1 (MQ=255)
                                                                       tGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTgg                                                                                <  2:282598/150‑1 (MQ=255)
                                                                         cAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGc                                                                              >  2:218955/1‑150 (MQ=255)
                                                                             tgTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATgg                                                                          <  2:397712/149‑1 (MQ=255)
                                                                             aGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATgg                                                                          <  1:719599/150‑1 (MQ=255)
                                                                                aaaaaGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCCCTATTGAAGAGATCCAGGCATCCCTTCCTGCAGAGTGGTCCAATAATGATATCCTGGTAAATAATGACGGCCTGGCGTTGGGCATGGAGc                                                                       >  1:1138279/2‑150 (MQ=255)
                                                                                       cGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCa                                                                <  2:547538/150‑1 (MQ=255)
                                                                                           cTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATaaa                                                            >  2:987515/1‑150 (MQ=255)
                                                                                            tGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAg                                                           >  2:879870/1‑150 (MQ=255)
                                                                                            tGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAg                                                           <  2:448718/150‑1 (MQ=255)
                                                                                            tGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAg                                                           <  2:334973/150‑1 (MQ=255)
                                                                                               gagaTAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCa                                                        >  2:1148355/1‑150 (MQ=255)
                                                                                                agaTAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATCCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGCCCTGGCGTTGGGCATGGACCCAGCGCATAAAGCCAg                                                       >  2:962271/1‑150 (MQ=255)
                                                                                                    atcttttataTCGCCCAACTGGACGTTCGCAACCGCGCCTCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTTGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGtt                                                   <  2:776259/144‑1 (MQ=255)
                                                                                                     aTCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTg                                                  <  2:1439026/150‑1 (MQ=255)
                                                                                                        tGTATATCGCCCAACTTGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCATCGTTGAAg                                               <  1:1032571/150‑1 (MQ=255)
                                                                                                        tGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAg                                               >  1:1470403/1‑150 (MQ=255)
                                                                                                        tGTATATCGCCCAACTGGACGTTCGCAACAGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGttagat                                               >  1:776259/1‑146 (MQ=255)
                                                                                                                gCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGaa                                       <  1:1004262/150‑1 (MQ=255)
                                                                                                                gCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGaa                                       >  1:301905/1‑150 (MQ=255)
                                                                                                                gCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGATGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGaa                                       >  1:1412007/1‑150 (MQ=255)
                                                                                                                   cAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACg                                    >  2:464675/1‑150 (MQ=255)
                                                                                                                     aCTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACgat                                  <  1:621696/150‑1 (MQ=255)
                                                                                                                          aCGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGAtgattg                             <  1:1299782/150‑1 (MQ=255)
                                                                                                                           cGTTCGCAACCGCGCCCCTATTGAAGAGAAGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGAtgattga                            >  2:48168/1‑150 (MQ=255)
                                                                                                                            gTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGAtgattgat                           >  2:79089/1‑150 (MQ=255)
                                                                                                                             ttCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGCGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATa                          <  1:1348776/150‑1 (MQ=255)
                                                                                                                             ttCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCGGCGTTGAAGACTGGGAAACGATGATTGATa                          >  2:273289/1‑150 (MQ=255)
                                                                                                                             ttCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATa                          >  2:739671/1‑150 (MQ=255)
                                                                                                                                gCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTgcg                                                                  >  1:680480/1‑107 (MQ=255)
                                                                                                                                gCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTgcg                                                                  <  2:680480/107‑1 (MQ=255)
                                                                                                                                  aaCCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACcaac                     <  2:701961/150‑1 (MQ=255)
                                                                                                                                   aCCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACcaaca                    <  1:1507483/150‑1 (MQ=255)
                                                                                                                                    ccGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACcaacaa                   <  1:1079714/150‑1 (MQ=255)
                                                                                                                                     cgcgCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACcaacaac                  <  2:113648/150‑1 (MQ=255)
                                                                                                                                         ccgcTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAg              <  1:1058459/150‑1 (MQ=255)
                                                                                                                                           gcTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGc            <  2:1470403/150‑1 (MQ=255)
                                                                                                                                           cggcaggAAGCGATGCCAGCATCTCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGc            >  1:431595/7‑150 (MQ=255)
                                                                                                                                              aTTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTg         >  1:988116/1‑150 (MQ=255)
                                                                                                                                                 gAAGAGATGCTGTCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTTTGCATGGAGAATGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCATGGta      >  2:349033/1‑150 (MQ=255)
                                                                                                                                                 gAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGta      <  1:218955/150‑1 (MQ=255)
                                                                                                                                                 gAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGta      >  2:1131119/1‑150 (MQ=255)
                                                                                                                                                   agagATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGATGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGtata    <  2:1412007/150‑1 (MQ=255)
                                                                                                                                                     agaTGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGtatata  >  2:517474/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGTATATA  >  NC_000913/1627555‑1627853

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