Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,543,231 A→C 7.2% M277R (ATG→AGG) ‡ narU ← nitrate/nitrite transporter NarU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,543,2310AC 7.2% 287.8 / 13.5 98M277R (ATG→AGG) ‡narUnitrate/nitrite transporter NarU
Reads supporting (aligned to +/- strand):  ref base A (49/42);  new base C (4/3);  total (53/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

CATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCG  >  NC_000913/1543083‑1543380
                                                                                                                                                    |                                                                                                                                                      
cATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCCCCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCAt                                                                                                                                                       >  1:789860/1‑150 (MQ=255)
 aTAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATg                                                                                                                                                      <  2:415556/150‑1 (MQ=255)
  tAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATgg                                                                                                                                                     >  1:620937/1‑150 (MQ=255)
  tAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAAAATTCACATCCGGGAACTGGGTTTTTGCCAGCATgg                                                                                                                                                     >  1:656914/1‑150 (MQ=255)
   taaaaTAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGc                                                                                                                                                    <  2:974905/149‑1 (MQ=255)
          aaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGGTTTTGCCATGCTGGCAAAAacc                                                                                                                                            <  1:1157467/150‑1 (MQ=255)
          aaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTGCCATGCTGGCAAAAaccc                                                                                                                                           <  2:323532/150‑2 (MQ=255)
          aaGTTGATCAACGTCACCCGCACGCCGCCGAACCTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTGCCATGCTGGCAAAAaccc                                                                                                                                           <  1:1031318/150‑2 (MQ=255)
           aGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAAcccc                                                                                                                                           >  2:812599/1‑149 (MQ=255)
            gTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGc                                                                                                                                          <  2:419475/150‑1 (MQ=255)
             ttGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCa                                                                                                                                         >  1:1014813/1‑150 (MQ=255)
             ttGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGCGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCa                                                                                                                                         <  1:243204/150‑1 (MQ=255)
              tGATCAACGTCAACCGCACGCCGCCGAACTTATCGGAATTAGCACCACCAACCGACCGCGCGATGGCCCCGATACATGGGCCAAAGAACGCCAGGCGCTGAATTATCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAACCGCGCAg                                                                                                                                        >  2:1240926/1‑150 (MQ=255)
               gATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCCCCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGa                                                                                                                                       >  1:139440/1‑150 (MQ=255)
               gATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGa                                                                                                                                       >  1:1324693/1‑150 (MQ=255)
               gATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGa                                                                                                                                       >  2:756957/1‑150 (MQ=255)
                aTCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGaa                                                                                                                                      >  1:289153/1‑150 (MQ=255)
                   aaCGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGaaaaa                                                                                                                                   >  2:1007598/1‑150 (MQ=255)
                    aCGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAAc                                                                                                                                  >  1:1267474/1‑150 (MQ=255)
                     cGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAAcc                                                                                                                                 <  1:933614/150‑1 (MQ=255)
                     cGTCACCAGCACGCCGCCGAACTTATCGGAAATACCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAAcc                                                                                                                                 <  2:404569/150‑1 (MQ=255)
                        cACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGAt                                                                                                                              >  2:740731/1‑150 (MQ=255)
                          cccGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATaa                                                                                                                            <  1:54133/150‑1 (MQ=255)
                            cGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAAc                                                                                                                          >  2:952043/1‑150 (MQ=255)
                            cGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAAc                                                                                                                          <  1:416438/150‑1 (MQ=255)
                            cGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAAc                                                                                                                          <  1:417828/150‑1 (MQ=255)
                                        aaCTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCAGCGATAAATGGGCCTAAGAAGGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAATCCGATAAACGAACCGAAGGTg                                                                                                              <  2:552750/150‑1 (MQ=255)
                                         aCTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTgg                                                                                                             <  1:740731/150‑1 (MQ=255)
                                         aCTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTgg                                                                                                             <  2:647117/150‑1 (MQ=255)
                                           ttATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCa                                                                                                           >  1:496662/1‑150 (MQ=255)
                                            tATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCaa                                                                                                          >  2:136407/1‑150 (MQ=255)
                                            tATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGCCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCaa                                                                                                          >  2:1155698/1‑150 (MQ=255)
                                             aTCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAg                                                                                                         <  1:412393/150‑1 (MQ=255)
                                             aTCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAg                                                                                                         <  2:579898/150‑1 (MQ=255)
                                               cGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGt                                                                                                       <  2:44025/150‑1 (MQ=255)
                                                ggAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCATAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTa                                                                                                      <  1:1313208/150‑1 (MQ=255)
                                                ggAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTa                                                                                                      <  2:1248019/150‑1 (MQ=255)
                                                ggAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTa                                                                                                      <  1:1007598/150‑1 (MQ=255)
                                                ggAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTa                                                                                                      >  1:511184/1‑150 (MQ=255)
                                                       gCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGACCGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAACCAgg                                                                                               >  1:638799/1‑150 (MQ=255)
                                                          ccaccaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTc                                                                                            <  1:140219/150‑1 (MQ=255)
                                                          ccaccaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTc                                                                                            >  1:1393089/1‑150 (MQ=255)
                                                           caccaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTca                                                                                           >  2:253993/1‑150 (MQ=255)
                                                               aaCCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTcagcag                                                                                       >  2:1240642/1‑150 (MQ=255)
                                                               aaCCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGAAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTcagcag                                                                                       >  2:1473078/1‑150 (MQ=255)
                                                                 ccgaccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGcc                                                                                     <  2:458393/150‑1 (MQ=255)
                                                                   gatcgcgcgATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCag                                                                                   <  1:660314/147‑1 (MQ=255)
                                                                   gaccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCag                                                                                   >  1:638731/1‑150 (MQ=255)
                                                                    accgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCaga                                                                                  <  1:756957/150‑1 (MQ=255)
                                                                     ccgCGCGATGGCACCGATAAATGGGCCAAAGAACGACAGGCGCAGAATATTGACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCagag                                                                                 >  2:329404/1‑150 (MQ=255)
                                                                       gcgcgATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGAt                                                                               >  1:1404002/1‑150 (MQ=255)
                                                                         gcgATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGc                                                                             <  2:927343/150‑1 (MQ=255)
                                                                         gcgATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGc                                                                             <  1:1216241/150‑1 (MQ=255)
                                                                             tGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGc                                                                         >  1:391666/1‑150 (MQ=255)
                                                                             tGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGc                                                                         >  1:1518406/1‑150 (MQ=255)
                                                                              ggCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCg                                                                        >  1:1222609/1‑150 (MQ=255)
                                                                              ggCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCg                                                                        <  2:511184/150‑1 (MQ=255)
                                                                                cACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGtt                                                                      >  1:13981/1‑150 (MQ=255)
                                                                                  ccGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGt                                                                    <  1:52028/150‑1 (MQ=255)
                                                                                   cGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTa                                                                   <  1:812599/150‑1 (MQ=255)
                                                                                    gATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                                  >  1:615411/1‑150 (MQ=255)
                                                                                    gATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGCTGGCAAGGTAAAGCAGGCTCAGCAGCCAGCGATGCAGGCGTTGTaa                                                                  >  2:52702/1‑150 (MQ=255)
                                                                                      tAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGTAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGa                                                                <  1:647274/150‑1 (MQ=255)
                                                                                        aaTGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACa                                                              <  2:1267474/150‑1 (MQ=255)
                                                                                           gggCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAggg                                                           <  2:1122837/150‑1 (MQ=255)
                                                                                             gCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAg                                                         >  1:623885/1‑150 (MQ=255)
                                                                                              ccAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGc                                                        <  2:638731/150‑1 (MQ=255)
                                                                                              ccAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGc                                                        <  1:927931/150‑1 (MQ=255)
                                                                                              ccAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGc                                                        >  2:937760/1‑150 (MQ=255)
                                                                                               cAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCt                                                       >  2:1325547/1‑150 (MQ=255)
                                                                                               cAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCt                                                       >  2:911610/1‑150 (MQ=255)
                                                                                                           ggCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTg                                           <  2:151372/150‑1 (MQ=255)
                                                                                                               cAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGgc                                       >  2:1280964/1‑150 (MQ=255)
                                                                                                                aGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGgcg                                      >  2:740494/1‑150 (MQ=255)
                                                                                                                aGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGgcg                                      <  2:1014813/150‑1 (MQ=255)
                                                                                                                    tatTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTg                                  <  2:1404002/150‑1 (MQ=255)
                                                                                                                        cacaTCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAAcg                              >  2:1282961/1‑149 (MQ=255)
                                                                                                                        cacaTCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACt                              <  2:1324693/150‑1 (MQ=255)
                                                                                                                        cacaTCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACt                              >  2:743063/1‑150 (MQ=255)
                                                                                                                          caTCCGGGAACTGGGTTTTTGCCAGAATCGCAAAACCCGCAGAAAAACCGATAAACGAAACGAACCTGGCAAGCTAAAGAAGGCTCAGCAGCCAGAGATGCAGGCGTAGTAAGACACGGAGCTGGTCGGCAATTGAGGCGCGTGAACTgg                            >  2:116606/1‑150 (MQ=255)
                                                                                                                           aTCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGc                           <  1:952043/150‑1 (MQ=255)
                                                                                                                            tCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCTGCAATTGAGGCGCGTGAACTGGCg                          <  2:13981/150‑1 (MQ=255)
                                                                                                                              cGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGat                        >  1:127685/1‑150 (MQ=255)
                                                                                                                                ggAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGatat                      <  1:566424/150‑1 (MQ=255)
                                                                                                                                 taaCTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtc                     <  2:642624/149‑1 (MQ=255)
                                                                                                                                     tGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattc                 >  2:132172/1‑150 (MQ=255)
                                                                                                                                      gggTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattca                >  2:460056/1‑150 (MQ=255)
                                                                                                                                      cgggttttGCCATGCTGGCAAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattc                 >  2:220985/5‑150 (MQ=255)
                                                                                                                                       gggttttGCCATGCTGGCAAAACCCGCAGAAAAACCGATAAACGCACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattcat               >  1:1409367/4‑150 (MQ=255)
                                                                                                                                       gggttttGCCATGCTGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattcat               >  1:1353403/4‑150 (MQ=255)
                                                                                                                                        ggttttGCCATGCTGGCAAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattcat               >  2:461163/3‑150 (MQ=255)
                                                                                                                                          ttttGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATccc            <  1:1155698/150‑1 (MQ=255)
                                                                                                                                             tGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGa         >  2:454116/1‑150 (MQ=255)
                                                                                                                                              gCCAGCATGGCAAAACCCGCAGAAAACCCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCNANAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGNACTGGCGATATCATTCATNCCTGAc        >  2:1212877/1‑150 (MQ=255)
                                                                                                                                              gCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAc        >  1:1423673/1‑150 (MQ=255)
                                                                                                                                              gCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAc        >  1:1217510/1‑150 (MQ=255)
                                                                                                                                               ccAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAcc       >  1:710342/1‑150 (MQ=255)
                                                                                                                                                  gCATGGCAAAACCCGCAGAAAACCCGATCAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTCAGACAGTGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCTTTCATCCCTGACCATg    <  2:656914/150‑1 (MQ=255)
                                                                                                                                                    aTGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATgcg  <  2:615411/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCG  >  NC_000913/1543083‑1543380

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