Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,283,383 G→A D89N (GAC→AAC)  agaD → PTS system glucosamine enzyme IID component 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,283,3830GA96.2% 531.3 / 3.1 157D89N (GAC→AAC) agaDPTS system glucosamine enzyme IID component 1
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base A (79/72);  total (82/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

Only 100 of 163 total aligned reads displayed.

GATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTAT  >  NC_000913/3283293‑3283471
                                                                                          |                                                                                        
gATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGac                                                                                        >  2:1256796/1‑93 (MQ=255)
gATAACCTCGAATTTATTAATACCCACCCGAATCTGGGCGGATTCCTGCTGGGGTTATTAATTTCGATGGAAGAAAAAGGAGCAAACAGCAac                                                                                        >  2:1068699/1‑93 (MQ=255)
 aTAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAaca                                                                                       <  1:3182528/93‑1 (MQ=255)
  tAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAacac                                                                                      <  1:4412160/93‑1 (MQ=255)
   aaCCTCGAATTTATTAATCCCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACAcc                                                                                     <  2:4180802/93‑1 (MQ=255)
   aaCCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACAcc                                                                                     >  2:2190676/1‑93 (MQ=255)
      cTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGTTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCAtt                                                                                  <  2:1977749/93‑1 (MQ=255)
      cTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCAtt                                                                                  >  1:1765507/1‑93 (MQ=255)
       tCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGGTGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTa                                                                                 >  2:2585120/1‑93 (MQ=255)
       tCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTa                                                                                 >  1:396977/1‑93 (MQ=255)
        cGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTaa                                                                                >  2:83616/1‑93 (MQ=255)
         gAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTaaa                                                                               <  1:3369700/93‑1 (MQ=255)
         gAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTaaa                                                                               >  2:2871530/1‑93 (MQ=255)
          aaTTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAg                                                                              <  1:603073/93‑1 (MQ=255)
            tttATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGc                                                                            >  1:1909078/1‑93 (MQ=255)
             ttattaATACCCCCCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGcc                                                                           <  2:72628/93‑1 (MQ=255)
             ttattaATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGcc                                                                           <  1:809339/93‑1 (MQ=255)
                ttaATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCa                                                                        >  1:1814942/1‑93 (MQ=255)
                 taATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCaa                                                                       >  1:1299600/1‑93 (MQ=255)
                   aTACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAg                                                                     >  1:3621088/1‑93 (MQ=255)
                    tACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGt                                                                    <  2:775463/93‑1 (MQ=255)
                    tACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGt                                                                    >  1:3413600/1‑93 (MQ=255)
                    tACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGt                                                                    <  1:484833/93‑1 (MQ=255)
                    tACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGt                                                                    >  2:3623533/1‑93 (MQ=255)
                     acccacccGAATCTGGTCGGATTCCTGATGGGTTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTg                                                                   <  1:4444187/93‑1 (MQ=255)
                      cccacccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTgg                                                                  <  1:4161088/93‑1 (MQ=255)
                        cacccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCa                                                                <  2:3309444/93‑1 (MQ=255)
                         ccccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCAc                                                               <  1:599510/92‑1 (MQ=255)
                         ccccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCAc                                                               <  2:3954817/92‑1 (MQ=255)
                            cGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGAGGCACTGt                                                            >  2:3535724/1‑93 (MQ=255)
                               aTCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTg                                                         >  1:1939561/1‑93 (MQ=255)
                               aTCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTg                                                         <  1:3755075/93‑1 (MQ=255)
                                tCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTgg                                                        >  1:4257206/1‑93 (MQ=255)
                                  tGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGcc                                                      >  1:3287823/1‑93 (MQ=255)
                                    gTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCa                                                    >  1:2826806/1‑93 (MQ=255)
                                     tCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCaa                                                   <  2:2797575/93‑1 (MQ=255)
                                      cGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAAt                                                  <  2:2185508/93‑1 (MQ=255)
                                       ggATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATc                                                 <  2:4248806/93‑1 (MQ=255)
                                        gATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCg                                                >  2:3459305/1‑93 (MQ=255)
                                         aTTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGc                                               >  2:3564920/1‑93 (MQ=255)
                                          ttCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGcc                                              <  1:100581/93‑1 (MQ=255)
                                           tCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCg                                             <  1:3309252/93‑1 (MQ=255)
                                            ccTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCgg                                            <  1:429362/93‑1 (MQ=255)
                                                aTGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGAtt                                        >  1:4124531/1‑93 (MQ=255)
                                                aTGGGGCTATTAATTTCGCTGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGAtt                                        >  2:4499876/1‑93 (MQ=255)
                                                 tGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTg                                       >  2:466690/1‑93 (MQ=255)
                                                 tGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTg                                       >  1:2448878/1‑93 (MQ=255)
                                                 tGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTg                                       <  2:195784/93‑1 (MQ=255)
                                                 tGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTg                                       >  1:1697598/1‑93 (MQ=255)
                                                  ggggTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTgg                                      >  2:132971/1‑93 (MQ=255)
                                                  ggggTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTgg                                      <  2:2370408/93‑1 (MQ=255)
                                                   gggTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGgc                                     <  2:1560160/93‑1 (MQ=255)
                                                    ggTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTG‑TTGGCCCAATCGCCGGGATTGgcga                                   <  2:148176/93‑1 (MQ=255)
                                                      ttattaATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGgcgat                                  <  1:30051/93‑1 (MQ=255)
                                                      ttattaATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGgcgat                                  >  2:2855198/1‑93 (MQ=255)
                                                       tattaATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGgcgatg                                 >  2:2826274/1‑93 (MQ=255)
                                                       tattaATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGgcgatg                                 >  2:681525/1‑93 (MQ=255)
                                                        attaATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGgcgatgc                                >  1:1543819/1‑93 (MQ=255)
                                                         ttaATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGgcgatgcg                               >  1:910319/1‑93 (MQ=255)
                                                          taATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGgcgatgcga                              >  1:1388096/1‑93 (MQ=255)
                                                          taATTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGgcgatgcga                              >  1:2117554/1‑93 (MQ=255)
                                                           aaTTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAAGCCTCAAAGTGGCACTGTTTGGCCCAATCGCTGAGAGTGgcgctgtcat                             >  2:2751118/1‑89 (MQ=255)
                                                            aTTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGAtt                            <  2:579661/93‑1 (MQ=255)
                                                            aTTTCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGAtt                            >  1:558595/1‑93 (MQ=255)
                                                             tttCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGAttt                           <  1:2187319/93‑1 (MQ=255)
                                                             tttCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGAttt                           >  2:3597943/1‑93 (MQ=255)
                                                              ttCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGAtttt                          >  1:1153547/1‑93 (MQ=255)
                                                               tCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTc                         <  2:4280677/93‑1 (MQ=255)
                                                               tCGATGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTc                         >  1:3471435/1‑93 (MQ=255)
                                                                  aTGGAAGAAAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTgg                      >  2:566385/1‑93 (MQ=255)
                                                                     gaagaaAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGttt                   >  1:2613973/1‑93 (MQ=255)
                                                                       agaaAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTAc                 <  1:4169106/93‑1 (MQ=255)
                                                                       agaaAAAGGAGAAAACCGCAACACCATTAAAGGCATCAAAGTGGCACTGTTTGGCCCAATAGCCGGGATTGGCGATGCGATTTTCTGGTTTAc                 >  1:1878879/1‑93 (MQ=255)
                                                                        gaaAAAGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACt                <  2:2938560/93‑1 (MQ=255)
                                                                         aaaaaGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACtt               <  1:2166564/93‑1 (MQ=255)
                                                                          aaaaGGAGAAAACCGCAACACCATTAAAGGCCTCAGAGTGGCGCTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACttt              >  2:2947369/1‑93 (MQ=255)
                                                                          aaaaGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACttt              <  2:1749278/93‑1 (MQ=255)
                                                                          aaaaGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACttt              <  1:2185599/93‑1 (MQ=255)
                                                                          aaaaGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACttt              >  1:320168/1‑93 (MQ=255)
                                                                           aaaGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTttg             >  1:922558/1‑93 (MQ=255)
                                                                           aaaGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTttg             >  1:4580793/1‑93 (MQ=255)
                                                                            aaGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTttgt            >  1:2654900/1‑93 (MQ=255)
                                                                            aaGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTttgt            <  1:3341442/93‑1 (MQ=255)
                                                                            aaGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTttgt            >  1:3905641/1‑93 (MQ=255)
                                                                            aaGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTttgt            >  1:1361557/1‑93 (MQ=255)
                                                                             aGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTttgtt           >  2:348254/1‑93 (MQ=255)
                                                                             aGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTttgtt           >  2:3941726/1‑93 (MQ=255)
                                                                             aGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTttgtt           <  2:3951205/93‑1 (MQ=255)
                                                                             aGGAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTttgtt           <  1:3905462/93‑1 (MQ=255)
                                                                              ggAGAAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTttgttg          <  1:4676203/93‑1 (MQ=255)
                                                                               gagaAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGc         <  2:2413926/93‑1 (MQ=255)
                                                                                agaAAACCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGcc        >  1:2166132/1‑93 (MQ=255)
                                                                                    aaCCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGAtt    >  2:3695499/1‑93 (MQ=255)
                                                                                     aCCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTa   <  1:2857015/93‑1 (MQ=255)
                                                                                     aCCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTa   <  1:1017809/93‑1 (MQ=255)
                                                                                     aCCGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTa   <  1:2316266/93‑1 (MQ=255)
                                                                                      ccGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTAt  <  1:2035410/93‑1 (MQ=255)
                                                                                      ccGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTAt  <  1:1464558/93‑1 (MQ=255)
                                                                                      ccGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTAt  >  2:4455539/1‑93 (MQ=255)
                                                                                      ccGCAACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTAt  <  2:2613642/93‑1 (MQ=255)
                                                                                          |                                                                                        
GATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTAT  >  NC_000913/3283293‑3283471

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑