Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 560,123 A→G T402A (ACG→GCG)  sfmD → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913560,1230AG95.9% 644.8 / 9.9 196T402A (ACG→GCG) sfmDputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (5/3);  new base G (100/88);  total (105/91)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 213 total aligned reads displayed.

GTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCATG  >  NC_000913/560034‑560213
                                                                                         |                                                                                          
gTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGgcgc                                                                                         >  2:2372084/1‑93 (MQ=255)
gTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGgcgc                                                                                         <  1:5146052/93‑1 (MQ=255)
 tAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCa                                                                                        >  2:5391096/1‑93 (MQ=255)
   gCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAAt                                                                                      <  1:1678360/93‑1 (MQ=255)
   gCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAAt                                                                                      >  2:4075211/1‑93 (MQ=255)
   gCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAAt                                                                                      <  1:3233449/93‑1 (MQ=255)
    cGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAAtt                                                                                     <  2:2056664/93‑1 (MQ=255)
     ggTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTa                                                                                    <  1:1470296/93‑1 (MQ=255)
      gTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTAt                                                                                   >  1:569996/1‑93 (MQ=255)
       tAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATc                                                                                  >  1:4377456/1‑93 (MQ=255)
       tAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATc                                                                                  <  2:4760590/93‑1 (MQ=255)
         aTAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTg                                                                                <  1:2969050/93‑1 (MQ=255)
          tAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGGAATTATCttc                                                                               >  1:230090/1‑91 (MQ=255)
          tAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGc                                                                               >  2:1133860/1‑93 (MQ=255)
          tAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGc                                                                               >  1:231807/1‑93 (MQ=255)
           aGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGcc                                                                              <  2:4590763/93‑1 (MQ=255)
            gTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCa                                                                             >  2:3850958/1‑93 (MQ=255)
               aGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAAtt                                                                          <  1:4060811/93‑1 (MQ=255)
                gCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTa                                                                         >  1:4258596/1‑93 (MQ=255)
                gCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTa                                                                         >  2:4294675/1‑93 (MQ=255)
                  aaTCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCGAATTacc                                                                       >  2:5679261/1‑92 (MQ=255)
                  aaTCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTaca                                                                       <  1:198546/93‑1 (MQ=255)
                   atcatcGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCCATTacac                                                                      >  1:2504429/1‑93 (MQ=255)
                   atcatcGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTacac                                                                      <  2:1537107/93‑1 (MQ=255)
                    tcatcGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACAcc                                                                     >  2:703574/1‑93 (MQ=255)
                      atcGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGc                                                                   >  2:1951696/1‑93 (MQ=255)
                      atcGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGc                                                                   >  2:2664967/1‑93 (MQ=255)
                       tcGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTTACAATTACACCGcc                                                                  >  2:1920515/1‑93 (MQ=255)
                       tcGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGcc                                                                  >  2:2980512/1‑93 (MQ=255)
                            tttctttTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCttttt                                                             >  1:4633136/1‑93 (MQ=255)
                            tttctttTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCttttt                                                             <  1:2821195/93‑1 (MQ=255)
                            tttctttTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCttttt                                                             >  1:5443592/1‑93 (MQ=255)
                              tctttTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTttat                                                           >  2:4205431/1‑93 (MQ=255)
                                ttttttCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTttatta                                                         <  1:374961/93‑1 (MQ=255)
                                ttttttCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTttatta                                                         >  2:1792485/1‑93 (MQ=255)
                                ttttttCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTttatta                                                         <  2:232926/93‑1 (MQ=255)
                                ttttttCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTttatta                                                         <  1:6005871/93‑1 (MQ=255)
                                ttttttCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTttatta                                                         <  2:4548605/93‑1 (MQ=255)
                                 tttttCAGGGTACGGCACTCGGCTGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAg                                                        <  1:3500814/93‑1 (MQ=255)
                                 tttttCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAg                                                        <  1:2013856/93‑1 (MQ=255)
                                 tttttCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCCAATTACACCGCCTTTTTATTAg                                                        <  1:2340160/93‑1 (MQ=255)
                                  ttttCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAgg                                                       >  1:3402943/1‑93 (MQ=255)
                                  ttttCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAgg                                                       <  1:50086/93‑1 (MQ=255)
                                    ttCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGc                                                     <  2:4140176/93‑1 (MQ=255)
                                     tCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCt                                                    <  2:3826491/93‑1 (MQ=255)
                                     tCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCt                                                    <  1:4573887/93‑1 (MQ=255)
                                       aGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTgg                                                  >  2:3789796/1‑93 (MQ=255)
                                        gggTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCTCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTggg                                                 >  1:5581506/1‑93 (MQ=255)
                                        gggTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTggg                                                 >  2:5860580/1‑93 (MQ=255)
                                           tACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGgcg                                              >  2:5577001/1‑93 (MQ=255)
                                              ggCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCaa                                           <  2:3335869/93‑1 (MQ=255)
                                                 aCTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAAtct                                        >  2:1477964/1‑93 (MQ=255)
                                                 aCTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAAtct                                        >  1:5131281/1‑93 (MQ=255)
                                                  cTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAAtctc                                       <  2:1031492/93‑1 (MQ=255)
                                                  cTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAAtctc                                       <  2:2714275/93‑1 (MQ=255)
                                                  cTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAAtctc                                       <  2:1244182/93‑1 (MQ=255)
                                                  cTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAAtctc                                       >  1:2731643/1‑93 (MQ=255)
                                                  cTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAAtctc                                       >  1:199859/1‑93 (MQ=255)
                                                    cggcggTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCgg                                     <  2:6007455/93‑1 (MQ=255)
                                                    cggcggTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCgg                                     <  1:1340211/93‑1 (MQ=255)
                                                    cggcggTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCgg                                     <  1:4801691/93‑1 (MQ=255)
                                                     ggcggTTTACCACAGGAATTTACTGCCTACGGCGGGGGGCAATTATCTGCCAATTACACCGCCTTGTTATTAGTGCTGGGGCGCAATCTCgtg                                    >  2:3144246/1‑91 (MQ=255)
                                                     ggcggTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCggg                                    >  1:4119317/1‑93 (MQ=255)
                                                      gcggTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGa                                   >  2:5312286/1‑93 (MQ=255)
                                                        ggTTTACCACAGGAATTTACTGCCTACGGCGGGGTGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAAc                                 >  1:776513/1‑93 (MQ=255)
                                                         gTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACt                                >  1:95981/1‑93 (MQ=255)
                                                          tttACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTg                               <  1:3559469/93‑1 (MQ=255)
                                                           ttACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTgg                              <  1:4526976/93‑1 (MQ=255)
                                                            tACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTATGGCTGGGGCGCAATCTCGGGCACTggg                             >  1:3920971/1‑93 (MQ=255)
                                                            tACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTggg                             <  2:501467/93‑1 (MQ=255)
                                                             aCCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTgggg                            >  1:728895/1‑93 (MQ=255)
                                                              ccACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGgc                           <  2:2875177/93‑1 (MQ=255)
                                                               cacaGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGgcg                          >  1:4873699/1‑93 (MQ=255)
                                                                acaGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGgcgc                         >  2:5590250/1‑93 (MQ=255)
                                                                 caGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGGGCa                        >  2:1121852/1‑93 (MQ=255)
                                                                  aGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGGGCAg                       >  1:1788811/1‑93 (MQ=255)
                                                                   ggAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGGGCAgg                      >  1:3477544/1‑92 (MQ=255)
                                                                   ggAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAgt                      >  2:5935490/1‑93 (MQ=255)
                                                                   ggAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAgt                      >  2:1647373/1‑93 (MQ=255)
                                                                    gAATTTACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAgtg                     <  2:4402252/93‑1 (MQ=255)
                                                                     aaTGTACTGCCTACGGCTGGGCGCAAGTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAgtgt                    <  1:3810528/93‑1 (MQ=255)
                                                                       tttACTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCg                  >  1:5422630/1‑93 (MQ=255)
                                                                           cTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTgg              <  1:2882024/93‑1 (MQ=255)
                                                                           cTGCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTgg              <  2:3528290/93‑1 (MQ=255)
                                                                             gCCTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGAt            >  1:5881395/1‑93 (MQ=255)
                                                                              ccTACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATg           <  2:5204896/93‑1 (MQ=255)
                                                                                tACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCCATCTCGGGGACTGGGGCGCAGTGTCGCTGGATgga         >  2:2821880/1‑91 (MQ=255)
                                                                                tACGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTa         >  2:5239494/1‑93 (MQ=255)
                                                                                tACGGCGGGACGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTa         >  1:1150136/1‑93 (MQ=255)
                                                                                tACGGCGGGACGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTa         >  1:4228355/1‑93 (MQ=255)
                                                                                 aCGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTaa        <  2:4870314/93‑1 (MQ=255)
                                                                                 aCGGCGGGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTaa        <  2:2931173/93‑1 (MQ=255)
                                                                                  cggcggGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAAc       <  1:6139987/93‑1 (MQ=255)
                                                                                  cggcggGGCGC‑CTTATCTGGCAATTACACCGCATTTTTTTTAGGGCTGGGGCGCACTCTCGGGGACTGGGGCGGCGTGGCGGTGGATGTATCg      >  2:4414180/1‑93 (MQ=255)
                                                                                  cggcggGACGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAAc       <  2:3124004/93‑1 (MQ=255)
                                                                                   ggcggGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGGGCAGTGTCGCTGGATGTAACg      >  2:3966746/1‑93 (MQ=255)
                                                                                   ggcggGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACg      >  2:5544534/1‑93 (MQ=255)
                                                                                     cggGGCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCa    <  1:392950/93‑1 (MQ=255)
                                                                                      ggggCGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCAt   <  1:2610306/93‑1 (MQ=255)
                                                                                      ggggCGCAATTATCAGCCAATTACACCGCCTTTCTATTAGGGCTGGGGCGCAAACTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCAt   <  2:4365473/93‑1 (MQ=255)
                                                                                         |                                                                                          
GTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCCAATTACACCGCCTTTTTATTAGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCATG  >  NC_000913/560034‑560213

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