Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,751,176 G→A L227F (CTT→TTT)  b1673 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,751,1760GA91.5% 686.5 / 30.5 223L227F (CTT→TTT) b1673hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (10/9);  new base A (112/92);  total (122/101)
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.29e-01

Only 100 of 230 total aligned reads displayed.

TACGTGCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAGTTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCTTTC  >  NC_000913/1751086‑1751264
                                                                                          |                                                                                        
tACGTGCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAAtt                                                                                        >  1:1987094/1‑93 (MQ=255)
 aCGTGCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGCCATGGTTGTTATTCGCACCAATAAATTc                                                                                       >  2:7576115/1‑93 (MQ=255)
  cGTGCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCa                                                                                      <  2:2656826/93‑1 (MQ=255)
  cGTGCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCa                                                                                      <  2:8934423/93‑1 (MQ=255)
  cGTGCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCa                                                                                      <  2:7428778/93‑1 (MQ=255)
      cTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTCCTTGGCACGACATGGTTGTTATTCGCACAAATAAATTCAGtca                                                                                  >  2:2318500/1‑93 (MQ=255)
      cTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGtca                                                                                  <  1:6993185/93‑1 (MQ=255)
       tGTCCAGCGTGACTTGGGATCTGCATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGtcat                                                                                 <  2:7201679/93‑1 (MQ=255)
       tGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGtcat                                                                                 >  2:4921356/1‑93 (MQ=255)
       tGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGtcat                                                                                 >  2:1118808/1‑93 (MQ=255)
        gTCCAGCGTGACTTGGGATCTGAGTACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGtcatc                                                                                >  2:7049392/1‑93 (MQ=255)
        gTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGtcatc                                                                                >  1:7214142/1‑93 (MQ=255)
        gTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGtcatc                                                                                <  1:9093074/93‑1 (MQ=255)
        gTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGtcatc                                                                                <  2:4991173/93‑1 (MQ=255)
         tCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGtcatca                                                                               <  2:2062192/93‑1 (MQ=255)
         tCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGtcatca                                                                               <  1:2208529/93‑1 (MQ=255)
          ccAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTcatcat                                                                              >  2:2224428/1‑93 (MQ=255)
          ccAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTcatcat                                                                              <  1:631841/93‑1 (MQ=255)
             gCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAGTTTAGTCATCATGTa                                                                           >  1:7967342/1‑93 (MQ=255)
              cGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTaa                                                                          <  1:1057169/93‑1 (MQ=255)
              cGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTaa                                                                          <  1:1702813/93‑1 (MQ=255)
                 gACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCa                                                                       <  2:8256090/93‑1 (MQ=255)
                   cTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCAtt                                                                     <  1:1329827/93‑1 (MQ=255)
                    ttGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTg                                                                    >  1:2716908/1‑93 (MQ=255)
                     tGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGATATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTga                                                                   <  1:6786931/93‑1 (MQ=255)
                      gggATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAGTTCAGTCATCATGTAATCATTgag                                                                  >  2:8355860/1‑93 (MQ=255)
                       ggATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTgaga                                                                 <  1:5131438/93‑1 (MQ=255)
                       ggATCTGAATACTCTGCCCACGACTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTgaga                                                                 <  2:5504107/93‑1 (MQ=255)
                       ggATCTGAATACTATGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGAACCAATAAGTTCAGTCATCATGTAATCATTgaga                                                                 <  1:4231420/93‑1 (MQ=255)
                        gATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAGTTCAGTCATCATGTAATCATTGAGAc                                                                <  1:503277/93‑1 (MQ=255)
                          tCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGc                                                              >  2:4332939/1‑93 (MQ=255)
                          tCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGc                                                              <  2:7260359/93‑1 (MQ=255)
                            tGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCtt                                                            >  2:2979026/1‑93 (MQ=255)
                            tGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCtt                                                            >  1:7676645/1‑93 (MQ=255)
                              aaTACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCa                                                          <  2:6100179/93‑1 (MQ=255)
                              aaTACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCa                                                          <  2:9031643/93‑1 (MQ=255)
                               aTACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCAt                                                         >  2:5578106/1‑93 (MQ=255)
                                  ctctGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTTAGTCATCATGTAATCATTGAGACGCTTCAactc                                                      <  1:2913073/93‑1 (MQ=255)
                                  ctctGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCAtctc                                                      <  2:1584174/93‑1 (MQ=255)
                                  ctctGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCAtctc                                                      >  1:2563405/1‑93 (MQ=255)
                                  ctctGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGGTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCAtctc                                                      >  2:6207959/1‑93 (MQ=255)
                                   tctGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAGTTCAGTCATCATGTAATCATTGAGACGCTTCATCTcc                                                     <  1:7651190/93‑1 (MQ=255)
                                   tctGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTcc                                                     <  2:7157523/93‑1 (MQ=255)
                                   tctGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTcc                                                     <  2:2027197/93‑1 (MQ=255)
                                    ctGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCt                                                    <  2:9630040/93‑1 (MQ=255)
                                       cccACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGAc                                                 >  1:4783156/1‑93 (MQ=255)
                                         cACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGa                                               <  1:4626016/93‑1 (MQ=255)
                                          aCGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGAt                                              >  1:9873291/1‑93 (MQ=255)
                                            gATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCg                                            <  1:2161599/93‑1 (MQ=255)
                                            gATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCg                                            >  1:9916675/1‑93 (MQ=255)
                                            gATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCg                                            >  1:7255535/1‑93 (MQ=255)
                                             aTTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCgg                                           >  1:2848699/1‑93 (MQ=255)
                                             aTTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCgg                                           >  2:4149628/1‑93 (MQ=255)
                                               tgtgGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCa                                         <  2:1087764/93‑1 (MQ=255)
                                                 tgGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAGTTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAAt                                       <  2:9516353/93‑1 (MQ=255)
                                                  ggCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATg                                      >  2:6760288/1‑93 (MQ=255)
                                                   gCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGt                                     <  2:3491199/93‑1 (MQ=255)
                                                    catACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGtt                                    <  1:535030/91‑1 (MQ=255)
                                                       aCTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGAc                                 <  2:1736240/93‑1 (MQ=255)
                                                        cTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGAcc                                <  2:6638615/93‑1 (MQ=255)
                                                         ttGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGAccc                               >  1:3395668/1‑93 (MQ=255)
                                                          tGGCACGACATGGTTGTTATTCGCACCAATAAGTTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGAcccc                              <  1:8851277/93‑1 (MQ=255)
                                                          tGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGAcccc                              <  2:8226820/93‑1 (MQ=255)
                                                           tgCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCt                             <  1:3546683/92‑1 (MQ=255)
                                                              aCGACATGGTTGTTATTCGCACCAATAAGTTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTc                          >  1:3843186/1‑93 (MQ=255)
                                                              aCGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTc                          >  2:3612800/1‑93 (MQ=255)
                                                              aCGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTc                          >  1:811048/1‑93 (MQ=255)
                                                              aCGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTc                          <  1:5901659/93‑1 (MQ=255)
                                                              aCGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTc                          >  1:6446685/1‑93 (MQ=255)
                                                              aCGACATGGATGATATTCGCACAAATAAATTCAGTCATCATGTAATCATAGAAAAACTTCAACTACTGACGAACGGCAATGTTGACCCCTTTc                          <  2:9244883/93‑1 (MQ=255)
                                                                gACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGt                        >  2:6741926/1‑93 (MQ=255)
                                                                  cATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTcc                      >  1:5063733/1‑93 (MQ=255)
                                                                   aTGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCTCTTTCGTccc                     >  2:615573/1‑93 (MQ=255)
                                                                    tGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCt                    >  1:3246746/1‑93 (MQ=255)
                                                                      gttgttATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTc                  <  1:9739623/93‑1 (MQ=255)
                                                                       ttgttATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCa                 >  1:458/1‑93 (MQ=255)
                                                                        tgttATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCaa                >  1:956027/1‑93 (MQ=255)
                                                                        tgttATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCaa                >  2:3152540/1‑93 (MQ=255)
                                                                         gttATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAAc               >  2:945118/1‑93 (MQ=255)
                                                                         gttATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAAc               <  2:6538382/93‑1 (MQ=255)
                                                                         gttATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAAc               <  1:5263135/93‑1 (MQ=255)
                                                                         gttATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAAc               >  1:1559507/1‑93 (MQ=255)
                                                                         gttATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAAc               >  1:5690040/1‑93 (MQ=255)
                                                                         gttATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAAc               >  2:2819322/1‑93 (MQ=255)
                                                                          ttATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAAcc              >  1:8774544/1‑93 (MQ=255)
                                                                             ttCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCa           >  2:9091742/1‑93 (MQ=255)
                                                                              tCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCaa          >  2:9291336/1‑93 (MQ=255)
                                                                                 cACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCg       >  1:7420748/1‑93 (MQ=255)
                                                                                  aCCAATAAGTTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGc      >  2:3637198/1‑93 (MQ=255)
                                                                                  aCCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCCATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGc      >  1:6385303/1‑93 (MQ=255)
                                                                                  aCCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGc      >  2:7624886/1‑93 (MQ=255)
                                                                                   ccAATAAGTTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCt     >  2:3605934/1‑93 (MQ=255)
                                                                                   ccAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCt     >  2:95459/1‑93 (MQ=255)
                                                                                    cAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCtt    >  2:8629982/1‑93 (MQ=255)
                                                                                    cAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCtt    >  2:4461412/1‑93 (MQ=255)
                                                                                    cAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCtt    >  2:4692275/1‑93 (MQ=255)
                                                                                    cAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCtt    >  2:6512434/1‑93 (MQ=255)
                                                                                    cAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCtt    >  1:8196921/1‑93 (MQ=255)
                                                                                     aaTAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCttt   >  1:8680514/1‑93 (MQ=255)
                                                                                      aTAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCTTTc  >  1:3751481/1‑93 (MQ=255)
                                                                                          |                                                                                        
TACGTGCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAGTTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCTTTC  >  NC_000913/1751086‑1751264

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