Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,003,383 (G)7→8 coding (184/681 nt) ycbF → putative chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,003,3761.G98.5% 456.0 / ‑0.5 134coding (177/681 nt)ycbFputative chaperone
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (114/18);  total (115/19)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 167 total aligned reads displayed.

GCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCT  >  NC_000913/1003284‑1003464
                                                                                              |                                                                                         
gCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGcc                                                                                            <  2:3897956/93‑1 (MQ=255)
 cagacagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑g                                                                                         <  2:2771793/93‑1 (MQ=255)
 cagacagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑g                                                                                         <  2:3981754/93‑1 (MQ=255)
  agacagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑gg                                                                                        >  1:3877228/1‑93 (MQ=255)
   gacagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑ggg                                                                                       >  1:1728262/1‑93 (MQ=255)
   gacagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑ggg                                                                                       <  2:3614968/93‑1 (MQ=255)
   gacagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑ggg                                                                                       <  2:4324921/93‑1 (MQ=255)
   gacagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑ggg                                                                                       >  1:4673553/1‑93 (MQ=255)
     cagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑ggggg                                                                                     >  1:2597950/1‑93 (MQ=255)
       gagaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑ggggggg                                                                                   >  1:5701294/1‑93 (MQ=255)
       gagaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑ggggggg                                                                                   >  1:1775661/1‑93 (MQ=255)
       gagaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑ggggggg                                                                                   >  2:5944524/1‑93 (MQ=255)
         gaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑gggggggga                                                                                 >  1:4640290/1‑91 (MQ=255)
          aaTCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑ggggggggaa                                                                                >  2:1056019/1‑90 (MQ=255)
           atctatatctatTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑ggggggggaag                                                                               <  1:2192721/93‑5 (MQ=255)
             cGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCg                                                                              >  1:4094688/1‑93 (MQ=255)
              gATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGt                                                                             >  2:4852740/1‑93 (MQ=255)
                tatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCg                                                                           >  2:5098489/1‑93 (MQ=255)
                 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGa                                                                          >  2:6188995/1‑93 (MQ=255)
                   cTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACa                                                                        >  1:6252647/1‑93 (MQ=255)
                   cTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACa                                                                        >  2:1485894/1‑93 (MQ=255)
                    tATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAg                                                                       >  2:4441099/1‑93 (MQ=255)
                    tATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAg                                                                       >  2:4876170/1‑93 (MQ=255)
                    tATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAg                                                                       >  1:5456596/1‑93 (MQ=255)
                     aTTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGt                                                                      >  2:2386551/1‑93 (MQ=255)
                          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAg                                                                 >  1:1153253/1‑93 (MQ=255)
                          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAg                                                                 >  1:4431273/1‑93 (MQ=255)
                            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGCGCCAgga                                                               >  2:5809408/1‑91 (MQ=255)
                             gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCAc                                                              >  2:1129744/1‑93 (MQ=255)
                             gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCAc                                                              >  2:5981754/1‑93 (MQ=255)
                               cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCa                                                            >  1:2819622/1‑93 (MQ=255)
                                 aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                                                          >  2:2203558/1‑93 (MQ=255)
                                   tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTAc                                                        >  1:3680817/1‑93 (MQ=255)
                                   tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTAc                                                        >  1:98596/1‑93 (MQ=255)
                                     ccAGTCGGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGATGCGTCGACAGTGCCAGCACCATTACTg                                                      >  1:5568901/1‑93 (MQ=255)
                                     ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                                                      >  1:4591870/1‑93 (MQ=255)
                                     ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                                                      >  1:2296026/1‑93 (MQ=255)
                                     ccAGTCAGGAACCCTGGCTTATAAACAGTAAAACCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                                                      <  2:1666831/93‑1 (MQ=255)
                                     caaGTCAGGAATCCTGTCTTATAAACAGAAAAATCACCCCCCAAACGCGTGGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                                                      <  1:6068141/91‑1 (MQ=255)
                                      cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg                                                     >  1:45169/1‑93 (MQ=255)
                                       aGTCAGAATCCCTGTTTTATAAAAAGAAAAACCACCCGCCCACCGGGTGGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGc                                                    <  1:5918657/93‑1 (MQ=255)
                                       aGTAAGGAATCCTGGCTTATAAACAGAAAAACCAACCGCCCAAGGGGTGGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGc                                                    <  2:736441/93‑1 (MQ=255)
                                        gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTggac                                                   >  1:4477237/1‑91 (MQ=255)
                                         tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCg                                                  >  1:4658101/1‑93 (MQ=255)
                                           aGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCt                                                <  2:5250578/93‑1 (MQ=255)
                                           aGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCt                                                >  1:4972915/1‑93 (MQ=255)
                                             gAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                              >  2:202936/1‑93 (MQ=255)
                                             gAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                              >  1:2289376/1‑93 (MQ=255)
                                             gAATCCTGGCTTATAAACAGTAAAATCACCCGCCCAGCGTTTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                              <  2:2939651/93‑1 (MQ=255)
                                                tCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGTAGCGTCGACAGTGCCAGCACCATTACTGGCGGCTCCGcc                                           >  2:6087931/1‑93 (MQ=255)
                                                tCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc                                           >  2:6024196/1‑93 (MQ=255)
                                                tCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc                                           >  2:5681178/1‑93 (MQ=255)
                                                 ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCAc                                         >  2:2063701/1‑93 (MQ=255)
                                                  cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCAc                                         >  1:5330622/1‑93 (MQ=255)
                                                  cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCAc                                         >  1:3575356/1‑93 (MQ=255)
                                                  cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCAc                                         >  2:5200298/1‑93 (MQ=255)
                                                     gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTa                                      >  1:1155176/1‑93 (MQ=255)
                                                       ttattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAtt                                    >  2:5390815/1‑93 (MQ=255)
                                                       ttattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAtt                                    >  2:5955941/1‑93 (MQ=255)
                                                        tattaACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtc                                   <  1:4009837/93‑1 (MQ=255)
                                                        tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAtt                                    >  1:4146015/1‑93 (MQ=255)
                                                        tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtc                                   >  1:1959933/1‑93 (MQ=255)
                                                         attaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtct                                  >  2:3358158/1‑93 (MQ=255)
                                                         attaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGATCCCCCACTTATtct                                  >  1:3764741/1‑93 (MQ=255)
                                                          ttaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtctc                                 >  2:1748008/1‑93 (MQ=255)
                                                          ttaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtctc                                 >  2:1600460/1‑93 (MQ=255)
                                                           taACAGTAAAACAAACCCCCCAAGGCGTTGGGGCG‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTcc                                <  2:4111904/93‑1 (MQ=255)
                                                             aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                              >  2:2323873/1‑93 (MQ=255)
                                                             aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                              >  1:515579/1‑93 (MQ=255)
                                                             aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                              >  1:3252652/1‑93 (MQ=255)
                                                               aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGaa                            >  2:6447724/1‑93 (MQ=255)
                                                               aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGaa                            >  1:2998050/1‑93 (MQ=255)
                                                                agaaaaTGCACCCCCCAACGGATTGGGGCAGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGaa                            <  1:499978/91‑1 (MQ=255)
                                                                  aaacTCACCCGCCCAACGCGTTGGCCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                         <  1:863168/93‑1 (MQ=255)
                                                                  aaaaTCAACCGCCCAACGCGTTGGGCCGGGGGGGGCAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                         >  1:2064231/1‑93 (MQ=255)
                                                                  aaaaTCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                         >  2:5543430/1‑93 (MQ=255)
                                                                  aaaaTCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                         >  1:922067/1‑93 (MQ=255)
                                                                   aaaTCAACCGCCCACCGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGccc                        <  1:2014685/93‑1 (MQ=255)
                                                                    aaTCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCg                       >  1:2844723/1‑93 (MQ=255)
                                                                      tCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGt                     >  1:2314956/1‑93 (MQ=255)
                                                                        aaCCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTag                   >  2:6064521/1‑92 (MQ=255)
                                                                        aaCCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTAc                   >  2:3708643/1‑93 (MQ=255)
                                                                         aCCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACg                  >  2:2495351/1‑93 (MQ=255)
                                                                          ccGCCCAACGGGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGa                 <  1:6448040/93‑1 (MQ=255)
                                                                          ccGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGa                 >  1:5685666/1‑93 (MQ=255)
                                                                          ccGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGa                 >  2:2100366/1‑93 (MQ=255)
                                                                             cccAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAg              >  1:6053249/1‑93 (MQ=255)
                                                                             cccAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAg              >  1:1629431/1‑93 (MQ=255)
                                                                              ccAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg             >  2:3075501/1‑93 (MQ=255)
                                                                              ccAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg             >  2:2608838/1‑93 (MQ=255)
                                                                                aaCGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa           >  2:3765279/1‑93 (MQ=255)
                                                                                aaCGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa           >  2:2034351/1‑93 (MQ=255)
                                                                                aaCGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa           >  2:3037493/1‑93 (MQ=255)
                                                                                aaCGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa           >  1:6448948/1‑93 (MQ=255)
                                                                                 aCGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa           <  1:1804054/92‑1 (MQ=255)
                                                                                  cgcgTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACg         >  2:2833180/1‑93 (MQ=255)
                                                                                     gTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTg      >  1:662974/1‑93 (MQ=255)
                                                                                      ttGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgc     >  2:621954/1‑93 (MQ=255)
                                                                                       tGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcg    >  1:22303/1‑93 (MQ=255)
                                                                                         ggccc‑‑GGGGGGGAAGCGTCGACAGTCCCAGCACCATTACTGGCCGCTCCGACACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCt  <  2:4025797/90‑1 (MQ=255)
                                                                                              |                                                                                         
GCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCT  >  NC_000913/1003284‑1003464

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