Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,134,669 T→A 6.6% intergenic (+112/‑60) flgF → / → flgG flagellar component of cell‑proximal portion of basal‑body rod/flagellar component of cell‑distal portion of basal‑body rod

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,134,6690TA 6.6% 312.0 / 7.2 104intergenic (+112/‑60)flgF/flgGflagellar component of cell‑proximal portion of basal‑body rod/flagellar component of cell‑distal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base T (50/47);  new base A (4/3);  total (54/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 104 total aligned reads displayed.

AGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCAAACGTC  >  NC_000913/1134524‑1134815
                                                                                                                                                  |                                                                                                                                                  
aGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                                <  1:1106093/150‑1 (MQ=255)
aGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGAGGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGAATTTGCGAGCACGTGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                                <  1:275991/150‑1 (MQ=255)
  gCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAGGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                              >  2:591108/1‑150 (MQ=255)
  gCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                              >  2:1674143/1‑150 (MQ=255)
  gCCGTGCAAACCAACTGCTGTCGAAGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATTCGCAAGTTTATCGGGTCTAT‑GGTTGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                              <  1:944044/150‑1 (MQ=255)
      tGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATGCGGCGTAAAcgccttatc                                                                                                                                          <  2:1608900/150‑6 (MQ=255)
       gCCAACCAACTGCAGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                         <  2:275429/150‑1 (MQ=255)
        ccaaccaaCTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                        >  2:1725852/1‑150 (MQ=255)
         caaccaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATAGGGGGCAATCGCAATTTATCGATTTTGCGAGCACATGTAGGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                        <  2:1746725/150‑1 (MQ=255)
         caaccaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                       >  2:511775/1‑150 (MQ=255)
         caaccaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                       <  1:1400415/150‑1 (MQ=255)
          aaccaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                      >  2:758954/1‑150 (MQ=255)
          aaccaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                      <  2:810551/150‑1 (MQ=255)
           accaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                     >  2:282016/1‑150 (MQ=255)
             caacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaag                                                                                                                                   >  2:218561/1‑150 (MQ=255)
             caacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaag                                                                                                                                   <  1:846776/150‑1 (MQ=255)
              aacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGATCACTTGTAGGCCGGATAAGGAGTTTACGCCGCATCCGGCaaga                                                                                                                                  >  1:496271/1‑150 (MQ=255)
                ctgctgTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaagaag                                                                                                                                >  2:12343/1‑150 (MQ=255)
                ctgctgTCGATGAGTTAATTGAAAGGATACATGACAAGTAAAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaagaag                                                                                                                                <  1:39730/150‑1 (MQ=255)
                  gctgTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGAc                                                                                                                              <  1:68070/150‑1 (MQ=255)
                   ctgTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACa                                                                                                                             >  2:1213551/1‑150 (MQ=255)
                   ctgTCGATCACTTAATTGAAAGGATACATGACAAGAATAAGTTGCCAGATCCGCAAGTTTATCGGGTCTAT‑GGGTGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACa                                                                                                                             <  2:868232/150‑1 (MQ=38)
                    tttCGATGAGTTAATTGAAAGTATACATGACAAGGATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACat                                                                                                                            <  2:1708247/148‑1 (MQ=255)
                     gTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGCCata                                                                                                                           >  2:1583956/1‑150 (MQ=255)
                       cGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATAGGGTCTAG‑GGGGGCAATCTAAATTTATCGAGATTGAAAGCACTTTTAGCCCGGATAAGACTTTTCCCACGCATAAGAAAAAAAGACattgg                                                                                                                         >  1:702633/1‑147 (MQ=255)
                            aGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACtt                                                                                                                    >  2:928367/1‑150 (MQ=255)
                             gTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACttt                                                                                                                   >  1:1088229/1‑150 (MQ=255)
                             gTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACttt                                                                                                                   >  1:1806294/1‑150 (MQ=255)
                               tAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGt                                                                                                                 >  2:506882/1‑150 (MQ=255)
                                aaTTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGta                                                                                                                >  1:967496/1‑149 (MQ=255)
                                 aTTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCa                                                                                                               >  1:153563/1‑150 (MQ=255)
                                 aTTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCa                                                                                                               >  2:538590/1‑150 (MQ=255)
                                  ttGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCAc                                                                                                              <  1:418364/150‑1 (MQ=255)
                                  ttGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCAc                                                                                                              >  2:1442019/1‑150 (MQ=255)
                                  ttGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCAc                                                                                                              >  2:1221768/1‑150 (MQ=255)
                                  ttGAAAGGATACATGACAAGTATAAGTTGCCCGATGAGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCAc                                                                                                              >  2:1538056/1‑150 (MQ=255)
                                   tGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACt                                                                                                             <  1:378940/150‑1 (MQ=255)
                                    gAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTa                                                                                                            >  1:979471/1‑150 (MQ=255)
                                          aTACATGACAAGTATAAGTTTCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCAc                                                                                                      <  1:295025/150‑1 (MQ=255)
                                          aTACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCAc                                                                                                      <  2:1710397/150‑1 (MQ=255)
                                           tACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACt                                                                                                     >  2:858257/1‑150 (MQ=255)
                                              aTGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACa                                                                                                  >  2:487879/1‑150 (MQ=255)
                                               tGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGCGGCTAT‑GGGGGCAATGGCAATTTATCGATTTTGCGAGCACTTGTAGCCCGGATAAGGCGTTTACTCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAg                                                                                                 <  1:347014/150‑1 (MQ=255)
                                                gACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGTCAATCGCAATTTATCGATTTTGCGAGGACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTATCACTAATCCACTACAgg                                                                                                <  2:1533147/150‑1 (MQ=255)
                                                    aGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACAt                                                                                            >  1:851527/1‑150 (MQ=255)
                                                     gTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACAtt                                                                                           >  2:480782/1‑150 (MQ=255)
                                                        taAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTa                                                                                        <  2:1164783/150‑1 (MQ=255)
                                                         aaGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTAt                                                                                       <  1:1725852/150‑1 (MQ=255)
                                                          aGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATg                                                                                      <  2:718802/150‑1 (MQ=255)
                                                          aGTTGACCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATg                                                                                      <  1:994506/150‑1 (MQ=255)
                                                           gTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCAGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGa                                                                                     >  1:1291844/1‑150 (MQ=255)
                                                             tGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATc                                                                                   <  2:1052508/150‑1 (MQ=255)
                                                             tGCCCGAGGCGCAAGTTTTTCGGGTCTAT‑GGGGGTAATCGCAATTTATCGATTTTGCGAGCACTTGTTGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTACTCCACTACAGGACATTTTATGATc                                                                                   <  2:1663109/150‑1 (MQ=39)
                                                                ccGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGt                                                                                >  1:127968/1‑150 (MQ=255)
                                                                  gATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTc                                                                              <  1:1271394/150‑1 (MQ=255)
                                                                  gATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTc                                                                              <  1:1636991/150‑1 (MQ=255)
                                                                      cgcAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTa                                                                          >  2:554635/1‑150 (MQ=255)
                                                                         aaGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATgg                                                                       <  2:1368075/150‑1 (MQ=255)
                                                                          aGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGa                                                                      <  1:416331/150‑1 (MQ=255)
                                                                           gTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGAt                                                                     <  2:1359218/150‑1 (MQ=255)
                                                                            tttATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATc                                                                    >  2:1120685/1‑150 (MQ=255)
                                                                             ttATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCg                                                                   >  2:510387/1‑150 (MQ=255)
                                                                              tATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGc                                                                  <  1:1559356/150‑1 (MQ=255)
                                                                              tATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGc                                                                  <  1:1558205/150‑1 (MQ=255)
                                                                                  gggTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCaaa                                                              >  1:1499962/1‑150 (MQ=255)
                                                                                    gTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAAc                                                            >  1:476156/1‑150 (MQ=255)
                                                                                    gTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAAc                                                            <  1:1641558/150‑1 (MQ=255)
                                                                                      cTAT‑GGTGGCAATCGCAAGTTATCGATTTTGCGATCACTTGTAGGCCGGATAAGGCTTTTACGCCGCATCCGGAAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACgg                                                          <  1:294191/150‑1 (MQ=255)
                                                                                       tAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACggg                                                         >  1:1516776/1‑150 (MQ=255)
                                                                                       tAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACggg                                                         <  1:1047729/150‑1 (MQ=255)
                                                                                       tAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACggg                                                         <  2:851527/150‑1 (MQ=255)
                                                                                       tAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACggg                                                         <  2:800231/150‑1 (MQ=255)
                                                                                       tAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACggg                                                         >  1:556222/1‑150 (MQ=255)
                                                                                            ggggCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCtt                                                     >  2:695476/1‑150 (MQ=255)
                                                                                             gggCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTg                                                    <  2:140959/150‑1 (MQ=255)
                                                                                                 aaTCGCAATTTATCGATTTTGCGAGCAGTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGc                                                <  2:675227/150‑1 (MQ=255)
                                                                                                   tCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGccc                                              <  1:1200374/150‑1 (MQ=255)
                                                                                                   tCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGccc                                              >  2:1086778/1‑150 (MQ=255)
                                                                                                    cGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCa                                             >  2:1396900/1‑150 (MQ=255)
                                                                                                      cAATTTATCGATTTTGCTAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGGCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGc                                           <  2:116211/150‑1 (MQ=255)
                                                                                                       aaTTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCa                                          <  1:1642165/150‑1 (MQ=255)
                                                                                                           tATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAAcc                                      <  2:1428045/150‑1 (MQ=255)
                                                                                                               gATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAata                                  >  1:371907/1‑150 (MQ=255)
                                                                                                               gATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAata                                  <  1:1538056/150‑1 (MQ=255)
                                                                                                                    tGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGAc                             <  1:1213551/150‑1 (MQ=255)
                                                                                                                      cGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCTAGCAAACCAATATGGACGt                           >  2:231564/1‑150 (MQ=255)
                                                                                                                          cACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCAtt                       >  2:158744/1‑150 (MQ=255)
                                                                                                                            cTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGCCGCCCAGCAAACCAATATGGACGTCATTGc                     >  1:178188/1‑150 (MQ=255)
                                                                                                                            cTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGc                     >  1:1594591/1‑150 (MQ=255)
                                                                                                                            cTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGc                     >  2:263317/1‑150 (MQ=255)
                                                                                                                            cTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGc                     >  2:1100961/1‑150 (MQ=255)
                                                                                                                               gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCcaa                  <  2:737782/150‑1 (MQ=255)
                                                                                                                                cttgCCGGATGCGGCGTAAACGCCTTATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCcaac                 >  1:1749413/4‑150 (MQ=255)
                                                                                                                                 ttgCCGGATGCGGCGTAAACGACTTATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCAATACAGGACATTTTATGATCAGTTCAGTATGGATCGCCAAAACGGGCCTTGGCGCCCAGCAAACCAATATGGACGTCATTGCcatca                >  1:1382412/3‑150 (MQ=255)
                                                                                                                                 aGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCcaaca                <  1:506882/150‑1 (MQ=255)
                                                                                                                                  tgCCGGATGCGGCGTAAACGCCTTATCCGGCAAGAAGCCAGATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTGCATTATGGATCGCCACAACGGGCCTTGCCGCCCAGCAAACCAATATGGACGTCATTGCcaacaa               >  1:99635/2‑150 (MQ=255)
                                                                                                                                      ggATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTg           <  1:471723/150‑1 (MQ=255)
                                                                                                                                         tAAGGCGTTTACGCCGCATCCGGCAAGAAGCCATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAAAATGGACGTCATTGCCAACAACCTGGCa        >  2:335171/1‑150 (MQ=255)
                                                                                                                                            ggCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCAAAc     >  2:1333137/1‑150 (MQ=255)
                                                                                                                                               gTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCAAACGTc  <  1:511775/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAT‑GGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCAAACGTC  >  NC_000913/1134524‑1134815

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑