Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,289,094 T→A 5.4% I444N (ATT→AAT)  nrfA → cytochrome c552 nitrite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,289,0940TA 5.4% 1043.7 / 28.9 313I444N (ATT→AAT) nrfAcytochrome c552 nitrite reductase
Reads supporting (aligned to +/- strand):  ref base T (126/170);  new base A (9/8);  total (135/178)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 316 total aligned reads displayed.

CCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGTG  >  NC_000913/4288946‑4289241
                                                                                                                                                    |                                                                                                                                                   
ccGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCAtt                                                                                                                                                    >  2:1556660/1‑150 (MQ=255)
ccGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCAtt                                                                                                                                                    >  1:990759/1‑150 (MQ=255)
  ggAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTgg                                                                                                                                                  >  2:462072/1‑150 (MQ=255)
    aagaagGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGcc                                                                                                                                                >  1:704300/1‑150 (MQ=255)
    aagaagGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGcc                                                                                                                                                <  1:1601723/150‑1 (MQ=255)
       aagGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGa                                                                                                                                             <  2:844506/150‑1 (MQ=255)
       aagGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGa                                                                                                                                             >  1:198820/1‑150 (MQ=255)
       aagGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGa                                                                                                                                             <  1:1394506/150‑1 (MQ=255)
         ggTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAAc                                                                                                                                           >  2:236021/1‑150 (MQ=255)
         ggTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAAc                                                                                                                                           <  1:502156/150‑1 (MQ=255)
           tttACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCAATGGCCTgctggg                                                                                                                                         <  2:1339988/150‑6 (MQ=255)
            ttACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATg                                                                                                                                        <  1:1326472/150‑1 (MQ=255)
                ggATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAAc                                                                                                                                    >  2:362284/1‑150 (MQ=255)
                 gctGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACa                                                                                                                                   <  2:1044848/148‑1 (MQ=255)
                    gCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAAt                                                                                                                                >  2:1849166/1‑150 (MQ=255)
                     cTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCt                                                                                                                                               <  2:303333/134‑1 (MQ=255)
                          tACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGc                                                                                                                          <  2:587896/150‑1 (MQ=255)
                              gCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCgag                                                                                                                      <  1:1811303/150‑1 (MQ=255)
                               cGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCgaga                                                                                                                     <  2:1303642/150‑1 (MQ=255)
                                   ggATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCa                                                                                                                 >  1:938275/1‑150 (MQ=255)
                                   ggATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCa                                                                                                                 >  1:266615/1‑150 (MQ=255)
                                   ggATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCa                                                                                                                 >  1:835884/1‑150 (MQ=255)
                                     aTAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAgg                                                                                                               <  2:1882465/150‑1 (MQ=255)
                                     aTAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAgg                                                                                                               <  2:1125094/150‑1 (MQ=255)
                                       aaaGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGAc                                                                                                             >  2:1597406/1‑150 (MQ=255)
                                             gcggATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTtcatc                                                                                                       <  1:1592312/150‑1 (MQ=255)
                                                    cacgcacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACgg                                                                                                >  1:666229/1‑150 (MQ=255)
                                                       gcacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGa                                                                                             >  2:257077/1‑150 (MQ=255)
                                                        cacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCa                                                                                                                              >  1:1147590/1‑116 (MQ=255)
                                                           cAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATccc                                                                                         <  2:1165512/150‑1 (MQ=255)
                                                             aaCTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGc                                                                                       <  2:198820/150‑1 (MQ=255)
                                                             aaCTGCCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGc                                                                                       <  1:566512/150‑1 (MQ=255)
                                                               cTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATcc                                                                                          >  1:1755790/1‑145 (MQ=255)
                                                               cTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAg                                                                                     <  2:1090734/150‑1 (MQ=255)
                                                                 ggCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTg                                                                                   <  1:682568/150‑1 (MQ=255)
                                                                 ggCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTg                                                                                   <  1:890674/150‑1 (MQ=255)
                                                                  gcgcgcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTgg                                                                                  <  1:462072/150‑1 (MQ=255)
                                                                   cgcgcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGaa                                                                                                                    >  1:1291201/1‑115 (MQ=255)
                                                                    gcgcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGga                                                                                <  1:1849166/150‑1 (MQ=255)
                                                                      gcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGgaag                                                                              <  2:1406033/150‑1 (MQ=255)
                                                                       ccTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGgaaga                                                                             >  2:378129/1‑150 (MQ=255)
                                                                        ctgctgGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATcc                                                                                          >  1:1386521/1‑136 (MQ=255)
                                                                        ctgctgGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATcc                                                                                          <  2:1386521/136‑1 (MQ=255)
                                                                          gctgGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACa                                                                          <  1:236021/150‑1 (MQ=255)
                                                                           ctgGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAg                                                                         <  1:771514/150‑1 (MQ=255)
                                                                                gACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACg                                                                    >  2:756656/1‑150 (MQ=255)
                                                                                gACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACg                                                                    <  2:1802888/150‑1 (MQ=255)
                                                                                 aCCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGGCTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGt                                                                   >  1:118155/1‑150 (MQ=255)
                                                                                  ccAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTa                                                                  >  1:1444020/1‑150 (MQ=255)
                                                                                    aaaGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTaaa                                                                >  1:1420519/1‑150 (MQ=255)
                                                                                      aGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTaaaaa                                                              >  2:540628/1‑150 (MQ=255)
                                                                                       ggCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAAc                                                             <  1:1910393/150‑1 (MQ=255)
                                                                                         cATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACgg                                                           <  1:1038745/150‑1 (MQ=255)
                                                                                           tCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTc                                                         <  1:756656/150‑1 (MQ=255)
                                                                                           tCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTc                                                         <  2:1613147/150‑1 (MQ=255)
                                                                                            cACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTcg                                                        >  1:1644937/1‑149 (MQ=255)
                                                                                              cccATGAAATCCAGATCGCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGATCATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGGCTGt                                                      >  2:1784114/1‑150 (MQ=255)
                                                                                               ccATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGTTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGtt                                                     <  2:1285212/150‑1 (MQ=255)
                                                                                               ccATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGtt                                                     <  1:796545/150‑1 (MQ=255)
                                                                                               ccATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGtt                                                     <  2:1261074/150‑1 (MQ=255)
                                                                                               ccATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGtt                                                     <  1:1827625/150‑1 (MQ=255)
                                                                                                cATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGTGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTa                                                    <  1:1471347/150‑1 (MQ=255)
                                                                                                cATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTa                                                    >  1:1674514/1‑150 (MQ=255)
                                                                                                 aTGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTaa                                                   >  1:383622/1‑150 (MQ=255)
                                                                                                    aaaTCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGcc                                                <  1:1556660/150‑1 (MQ=255)
                                                                                                        ccAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGGAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATa                                            >  2:1483578/1‑150 (MQ=255)
                                                                                                        ccAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATa                                            <  2:775508/150‑1 (MQ=255)
                                                                                                          aGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAAc                                          >  2:810625/1‑150 (MQ=255)
                                                                                                           caTCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAAcc                                         <  1:1260996/149‑1 (MQ=255)
                                                                                                              cccGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTaaa                                                                >  1:1530690/1‑124 (MQ=255)
                                                                                                              cccGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTaaa                                                                <  2:1530690/124‑1 (MQ=255)
                                                                                                              cccGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCg                                      >  2:1100728/1‑150 (MQ=255)
                                                                                                                cGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTT                                    <  2:1425614‑M1/150‑1 (MQ=255)
                                                                                                                     atCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCC                               >  2:1135046‑M1/1‑150 (MQ=255)
                                                                                                                      tctcAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCt                              >  1:1743228/1‑150 (MQ=255)
                                                                                                                       ctcAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGcc                                                >  2:710967/1‑131 (MQ=255)
                                                                                                                       ctcAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTc                             <  2:1008341/150‑1 (MQ=255)
                                                                                                                        tcAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCcct                            <  2:651847‑M1/150‑4 (MQ=255)
                                                                                                                        tcAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCg                            >  1:362774/1‑150 (MQ=255)
                                                                                                                            ccAAAGAGAGAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAg                        <  2:1741529/150‑1 (MQ=255)
                                                                                                                            ccAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAg                        >  1:1103166/1‑150 (MQ=255)
                                                                                                                              aaaGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGg                      <  1:332172‑M2/150‑2 (MQ=255)
                                                                                                                                   gaAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCccttgcgaggcgga                 >  2:892035‑M1/1‑136 (MQ=255)
                                                                                                                                     aaGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCccttgcgaggcggaac               <  1:1069751‑M1/150‑17 (MQ=255)
                                                                                                                                     aaGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGaaa               >  2:1195157/1‑150 (MQ=255)
                                                                                                                                       gCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAac             <  1:1002609/150‑1 (MQ=255)
                                                                                                                                        cccAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGggcgg                  <  1:1356969‑M2/144‑6 (MQ=255)
                                                                                                                                        cccAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAaca            <  2:923389/150‑1 (MQ=255)
                                                                                                                                        atgttcagGCCAATGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAaca            >  2:1661544/6‑150 (MQ=255)
                                                                                                                                         tgttcagGCCAATGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAacac           >  2:96236/5‑150 (MQ=255)
                                                                                                                                         tgttcagGCCAATGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAacac           >  1:626071/5‑150 (MQ=255)
                                                                                                                                          cagcagGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCccttgcgaggcggaacggggt          <  1:917782‑M1/150‑22 (MQ=255)
                                                                                                                                          cagcagGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAAcaca          >  1:85213/1‑150 (MQ=255)
                                                                                                                                          cagcagGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAAcaca          >  1:113026/1‑150 (MQ=255)
                                                                                                                                           agcagGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACaa         <  1:1485967/150‑1 (MQ=255)
                                                                                                                                             cagGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATg       >  1:1868930/1‑150 (MQ=255)
                                                                                                                                             cagGCCAATGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTg                                                       >  2:442377/1‑102 (MQ=255)
                                                                                                                                             cagGCCAATGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTg                                                       <  1:442377/102‑1 (MQ=255)
                                                                                                                                                 ccATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGt   <  2:1674514/150‑1 (MQ=255)
                                                                                                                                                  cATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGTg  <  2:743416/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGTG  >  NC_000913/4288946‑4289241

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: ‑