Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,130,038 C→T 5.1% A31T (GCA→ACA)  flgM ← anti‑sigma factor for FliA (sigma 28)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,130,0380CT 5.1% 260.1 / 2.9 78A31T (GCA→ACA) flgManti‑sigma factor for FliA (sigma 28)
Reads supporting (aligned to +/- strand):  ref base C (37/37);  new base T (2/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

CGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCA  >  NC_000913/1129891‑1130180
                                                                                                                                                   |                                                                                                                                              
cGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTgc                                                                                                                                                <  6:218378/148‑1 (MQ=255)
 ggTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTgcc                                                                                                                                               >  5:295386/1‑148 (MQ=255)
 ggTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTCCCGGGTTGAATCAGTTTTGCTTGCGCGTCGATTAACGTCACACTGGAGCTGGTGGAGGCGGTTGTTTTTgca                                                                                                                                               >  6:337430/1‑147 (MQ=255)
        aTTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGc                                                                                                                                        >  6:366276/1‑148 (MQ=255)
          ttttAGTTCACCGTTACGAATCGCCCGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTg                                                                                                                                      <  8:342856/148‑1 (MQ=255)
             tAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGGTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGtt                                                                                                                                    <  3:194346/147‑1 (MQ=255)
                  cACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTa                                                                                                                               <  8:335872/147‑1 (MQ=255)
                      gTTACGAATCGCCAGTTTTAACGCTTCGACACGGTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGgc                                                                                                                          <  4:212506/148‑1 (MQ=255)
                            aagcGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCAGCAGTTTTGCTTGCGCGTCGCTTAACGGCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCa                                                                                                                    <  5:337430/145‑1 (MQ=255)
                             aTCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGcgcg                                                                                                                       >  6:132130/1‑144 (MQ=255)
                             aTCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGcgcg                                                                                                                       <  5:132130/144‑1 (MQ=255)
                               cGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTg                                                                                                                 >  8:156209/1‑148 (MQ=255)
                                 ccAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGATTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTgc                                                                                                                                                >  7:236211/1‑115 (MQ=255)
                                     ttttAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTcg                                                                                                           >  8:153546/1‑148 (MQ=255)
                                      tttAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTcgc                                                                                                          >  5:111016/1‑148 (MQ=255)
                                      tttAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTcgc                                                                                                          >  8:342378/1‑148 (MQ=255)
                                           cGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGt                                                                                                     >  6:246639/1‑148 (MQ=255)
                                           cGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGt                                                                                                     <  5:366276/148‑1 (MQ=255)
                                           cGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGt                                                                                                     <  6:295386/148‑1 (MQ=255)
                                                cGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGaa                                                                                                 <  7:153546/147‑1 (MQ=255)
                                                 gACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACg                                                                                               <  7:342378/148‑1 (MQ=255)
                                                     cGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGc                                                                                           >  2:306003/1‑148 (MQ=255)
                                                      gTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCt                                                                                          >  4:343559/1‑148 (MQ=255)
                                                        tCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAg                                                                                                                   <  4:15668/121‑1 (MQ=255)
                                                        tCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAg                                                                                                                   >  3:15668/1‑121 (MQ=255)
                                                        tCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTa                                                                                        >  3:200818/1‑148 (MQ=255)
                                                            gactGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAgg                                                                                    >  7:133392/4‑148 (MQ=255)
                                                             aTTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGGGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGACCGGTGCTTACAGGc                                                                                   >  3:292254/1‑148 (MQ=255)
                                                               tGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGt                                                                                                                                 <  6:3209/100‑1 (MQ=255)
                                                               tGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGt                                                                                                                                 >  5:3209/1‑100 (MQ=255)
                                                                 atatCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTc                                                                                <  8:133392/147‑1 (MQ=255)
                                                                  tatCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCag                                                                              <  6:111016/148‑1 (MQ=255)
                                                                  tatCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCag                                                                              <  5:246639/148‑1 (MQ=255)
                                                                   atCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCTTCGCTTTACGTCTCACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTAGCGGCTCGTCAGTGGTTTCGCGCGGTTGAACGGTGc                                                                                           <  8:30181/134‑1 (MQ=255)
                                                                    tCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCtt                                                                                 <  7:75236/143‑1 (MQ=255)
                                                                    tCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCtt                                                                                 >  8:75236/1‑143 (MQ=255)
                                                                         gctgcCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGc                                                                                           <  3:83437/128‑1 (MQ=255)
                                                                         gctgcCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGc                                                                                           >  4:83437/1‑128 (MQ=255)
                                                                         gctgcCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCg                                                                         >  5:241896/1‑146 (MQ=255)
                                                                                  ttGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCa                                                              <  1:306003/148‑1 (MQ=255)
                                                                                   tGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCa                                                              >  6:320083/1‑147 (MQ=255)
                                                                                   tGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGACCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCaa                                                             >  6:123897/1‑148 (MQ=255)
                                                                                      aTCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGc                                                                                           >  6:380040/1‑115 (MQ=255)
                                                                                      aTCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGc                                                                                           <  5:380040/115‑1 (MQ=255)
                                                                                       tCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGa                                                                        >  1:184336/1‑133 (MQ=255)
                                                                                       tCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGa                                                                        <  2:184336/133‑1 (MQ=255)
                                                                                            tttGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGTCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGa                                                                                                  >  7:176841/1‑102 (MQ=255)
                                                                                            tttGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGTCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGa                                                                                                  <  8:176841/102‑1 (MQ=255)
                                                                                            tttGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGTCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGa                                                                        <  6:168765/128‑1 (MQ=255)
                                                                                            tttGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGTCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGa                                                                        >  5:168765/1‑128 (MQ=255)
                                                                                            tttGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCtt                                                                                         <  8:259109/111‑1 (MQ=255)
                                                                                            tttGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCtt                                                                                         >  7:259109/1‑111 (MQ=255)
                                                                                             ttgcttgcGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGa                                                                 >  3:13354/1‑134 (MQ=255)
                                                                                             ttgcttgcGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGa                                                                 <  4:13354/134‑1 (MQ=255)
                                                                                                cttgcGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGGTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGttt                                                 <  1:81200/147‑1 (MQ=255)
                                                                                                  tgcGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCtt                                                                                         >  5:367398/1‑105 (MQ=255)
                                                                                                  tgcGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCtt                                                                                         <  6:367398/105‑1 (MQ=255)
                                                                                                   gcgcGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTAt                                               <  3:343559/146‑1 (MQ=255)
                                                                                                           cTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGa                                     <  4:200818/148‑1 (MQ=255)
                                                                                                            ttAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTc                                                                                                                     >  7:58363/1‑67 (MQ=255)
                                                                                                            ttAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTc                                                                                                                     <  8:58363/67‑1 (MQ=255)
                                                                                                            ttAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCtt                                                                                         <  4:193992/95‑1 (MQ=255)
                                                                                                            ttAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCtt                                                                                         >  3:193992/1‑95 (MQ=255)
                                                                                                             tAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGaa                                                                       >  8:352345/1‑112 (MQ=255)
                                                                                                             tAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGaa                                                                       <  7:352345/112‑1 (MQ=255)
                                                                                                                   cacaCTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCtt                             >  8:200606/1‑148 (MQ=255)
                                                                                                                   cacaCTGGTGCTGGTGGAGGCGGTAGTTTTTGCCGCACGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGATGCGAAGTGCGATCAATACTCATCGTTTATTCCTCATTGAGGGCGCtt                             >  3:59252/1‑148 (MQ=255)
                                                                                                                        ttgtgctggtgGTGGCGGTTGTTTTTGCCGCCAGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCTGTTGAACGGTGCTTACAGGCTTCAGAGGAGAAGTTCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATc                        <  4:59252/146‑1 (MQ=255)
                                                                                                                        tggtgctggtgGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTAt                         >  6:209909/1‑147 (MQ=255)
                                                                                                                         ggtgctggtgGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCa                       >  1:244065/1‑148 (MQ=255)
                                                                                                                          gtgctggtgGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATc                        <  2:244065/146‑1 (MQ=255)
                                                                                                                                 tggAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATc                                                               >  2:374572/1‑100 (MQ=255)
                                                                                                                                 tggAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATc                                                               <  1:374572/100‑1 (MQ=255)
                                                                                                                                  ggAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCt              >  5:328275/1‑148 (MQ=255)
                                                                                                                                  ggAGGCGGATGTTTTTGCCGCCCGGCTCTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGACCGGTGCTTACACGATTCAGAGGCGAAGTGCGATCAATACTCATGGTGTATTCCGCAATGAGGGCGCATTTATCATCTGTAGCt              >  6:211730/1‑148 (MQ=255)
                                                                                                                                         gttgttTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCttatttat       <  6:328275/148‑1 (MQ=255)
                                                                                                                                          ttgttTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGt                                                   <  6:23278/103‑1 (MQ=255)
                                                                                                                                          ttgttTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGt                                                   >  5:23278/1‑103 (MQ=255)
                                                                                                                                          ttgttTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATc      >  6:112309/1‑148 (MQ=255)
                                                                                                                                           tgttTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCg     >  3:254576/1‑148 (MQ=255)
                                                                                                                                             tttttGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCgg    <  5:320083/147‑1 (MQ=255)
                                                                                                                                               tttGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCa  <  7:200606/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCA  >  NC_000913/1129891‑1130180

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