Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,438,883 C→T 8.1% R312C (CGT→TGT)  flk → putative flagella assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,438,8830CT 8.1% 218.6 / 7.0 74R312C (CGT→TGT) flkputative flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base C (31/36);  new base T (4/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAACCA  >  NC_000913/2438737‑2439026
                                                                                                                                                  |                                                                                                                                               
cgcAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGTGTc                                                                                                                                              <  1:1007316/150‑1 (MQ=255)
          gCTTTCTTCTGCCCAGGCCGAGTATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCtgtg                                                                                                                                    <  1:260169/150‑1 (MQ=255)
          gCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATGGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTTCAAAATATGGCAACGCGTCCCGCGCtgtg                                                                                                                                    <  2:1006725/150‑1 (MQ=255)
             ttcttcTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGAt                                                                                                                                 <  1:1392497/150‑1 (MQ=255)
              tcttctGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATa                                                                                                                                >  1:1262291/1‑150 (MQ=255)
                 tctGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTg                                                                                                                             <  2:1262291/150‑1 (MQ=255)
                  ctGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGt                                                                                                                            <  1:841136/150‑1 (MQ=255)
                   tGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGtt                                                                                                                           <  1:717510/150‑1 (MQ=255)
                   tGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGtt                                                                                                                           <  2:585883/150‑1 (MQ=255)
                    gCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTa                                                                                                                          >  1:709817/1‑150 (MQ=255)
                     cccAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAg                                                                                                                         <  1:1326743/150‑1 (MQ=255)
                      ccAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGt                                                                                                                        >  2:720906/1‑150 (MQ=255)
                                      cTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATg                                                                                                        <  1:1127603/150‑1 (MQ=255)
                                        tAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGAt                                                                                                                  <  2:954103/138‑1 (MQ=255)
                                        tAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGAt                                                                                                                  >  1:954103/1‑138 (MQ=255)
                                        tAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCt                                                                                                      <  2:842543/150‑1 (MQ=255)
                                         aaTCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTg                                                                                                     >  2:534218/1‑150 (MQ=255)
                                          aTCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTgg                                                                                                    >  1:462835/1‑150 (MQ=255)
                                          aTCAGCTTTATCAACGCCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATACTGATGCTgg                                                                                                    >  2:613571/1‑150 (MQ=255)
                                           tCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGt                                                                                                   >  1:598391/1‑150 (MQ=255)
                                              gCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTctg                                                                                                >  1:1112299/1‑150 (MQ=255)
                                               cTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTCGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTctgg                                                                                               <  2:530616/150‑1 (MQ=255)
                                                 ttATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTctggct                                                                                             <  1:1035849/150‑1 (MQ=255)
                                                 ttATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGGTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTctggct                                                                                             <  2:1155146/150‑1 (MQ=255)
                                                  tATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTctggctg                                                                                            >  2:310778/1‑150 (MQ=255)
                                                  tATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTctgg                                                                                               >  2:229743/1‑147 (MQ=255)
                                                  tATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTctgg                                                                                               <  1:229743/147‑1 (MQ=255)
                                                   aTCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTctggctgg                                                                                           >  2:590238/1‑150 (MQ=255)
                                                      aaCGGCGGGTTAAAGGGATTGATCCGCGTGTTATGAAACCGCTGCTTAATACTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGTTCTGGCTGGTTc                                                                                        >  2:341728/1‑150 (MQ=255)
                                                       aCGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCg                                                                                       <  2:801920/150‑1 (MQ=255)
                                                         ggcggGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGtt                                                                                     <  2:392243/150‑1 (MQ=255)
                                                             ggTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAAcc                                                                                 >  2:375353/1‑150 (MQ=255)
                                                                 aaaGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCcga                                                                             >  2:769100/1‑150 (MQ=255)
                                                                 aaaGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCcga                                                                             >  2:1223747/1‑150 (MQ=255)
                                                                  aaGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCcgac                                                                            <  1:1013371/150‑1 (MQ=255)
                                                                      gATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGaaa                                                                        >  2:594977/1‑150 (MQ=255)
                                                                       aTTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAAt                                                                       >  1:1020087/1‑150 (MQ=255)
                                                                       aTTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGGGGATACTGTAGTCGCGATTATCCTGAGGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATg                                                                      <  2:774159/150‑1 (MQ=255)
                                                                           aTCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGGGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACa                                                                   <  1:679969/150‑1 (MQ=255)
                                                                             ccGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGt                                                                 >  1:1372783/1‑150 (MQ=255)
                                                                                   gtTATGCAACCGCTGCTTAATCCTTTTCCACCGATGCTGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACACCGACGAAATGACAGTATCGTg                                                           >  2:76672/1‑150 (MQ=255)
                                                                                   gtTATGCAACCGCTGCATAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTg                                                           <  2:684085/150‑1 (MQ=255)
                                                                                    ttATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGa                                                          <  1:1120386/150‑1 (MQ=255)
                                                                                     tATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATTTGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTTTCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGAc                                                         >  2:1299727/1‑150 (MQ=255)
                                                                                     tATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGAc                                                         >  2:1406499/1‑150 (MQ=255)
                                                                                                tgctTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATAcccg                                              <  2:1401869/150‑2 (MQ=255)
                                                                                                 gctTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAg                                             >  2:603027/1‑150 (MQ=255)
                                                                                                    tAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCAc                                          >  1:67518/1‑150 (MQ=255)
                                                                                                    tAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCAc                                          >  1:1295598/1‑150 (MQ=255)
                                                                                                         cTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCag                                     <  1:1223747/150‑1 (MQ=255)
                                                                                                         cTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCag                                     <  2:709817/150‑1 (MQ=255)
                                                                                                          ttttCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCaga                                    <  1:310778/150‑1 (MQ=255)
                                                                                                            ttCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAgaga                                  <  2:517526/150‑1 (MQ=255)
                                                                                                            ttCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAgaga                                  <  2:1372783/150‑1 (MQ=255)
                                                                                                              ccACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATc                                >  2:76322/1‑150 (MQ=255)
                                                                                                               cACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCg                               >  1:73754/1‑150 (MQ=255)
                                                                                                                aCCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGc                              <  2:598391/150‑1 (MQ=255)
                                                                                                                 ccGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGcc                             >  2:1031046/1‑150 (MQ=255)
                                                                                                                        gggACGCGTTGCCATATTTTGCAACGTGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAg                                             >  2:52949/2‑127 (MQ=39)
                                                                                                                        gggACGCGTTGCCATATTTTGCAACGTGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAg                                             <  1:52949/126‑1 (MQ=39)
                                                                                                                        gggACGCGTTGCCATATTTTGCAACGTGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCaa                      >  2:375339/2‑150 (MQ=255)
                                                                                                                         ggACGCGTTGCCATATTTTGCAACGTGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAg                     >  1:1336403/1‑150 (MQ=255)
                                                                                                                         ggACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAg                     <  1:590238/150‑1 (MQ=255)
                                                                                                                          gACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGa                    <  1:218016/150‑1 (MQ=255)
                                                                                                                            cgcGTTGCCATATTTTGCAACGTGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGaaa                  >  1:308643/3‑150 (MQ=255)
                                                                                                                             acGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCag                                     >  1:217159/1‑130 (MQ=255)
                                                                                                                             acGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCag                                     <  2:217159/130‑1 (MQ=255)
                                                                                                                              agTTGCAAAATATGGCAACGCGTCCCGCGATGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTCGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCAACGCAGAGATCGACCCGGCAAGAAATa                >  2:95296/2‑150 (MQ=255)
                                                                                                                               gTTTCAAAATATGGCAACGAGACCCGCGCAGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATAc               >  2:669417/1‑150 (MQ=255)
                                                                                                                                    aaaaTATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAag          <  1:12715/150‑1 (MQ=255)
                                                                                                                                    aaaaTATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAag          <  1:720906/150‑1 (MQ=255)
                                                                                                                                     aaaTATGGCAACGCGTCACCCGCTGTGGATACTGTTAGTCGCGATTATCATGATGCTGGTCTGGCTGGTTCGTAAACCCAGAAGAAATTACAGTATCGTGACAATAATACCCAGCACAGCAGAGATCGCCCCCGCAAGAAATACCGAaga         >  2:1406492/1‑150 (MQ=255)
                                                                                                                                       atatGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGAAGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGt       <  2:462835/150‑1 (MQ=255)
                                                                                                                                            gCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAACCa  <  1:375353/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAACCA  >  NC_000913/2438737‑2439026

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