Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,309,009 C→T 6.2% intergenic (‑85/+24) nlpI ← / ← pnp lipoprotein involved in osmotic sensitivity and filamentation/polynucleotide phosphorylase/polyadenylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,309,0090CT 6.2% 379.9 / 10.2 115intergenic (‑85/+24)nlpI/pnplipoprotein involved in osmotic sensitivity and filamentation/polynucleotide phosphorylase/polyadenylase
Reads supporting (aligned to +/- strand):  ref base C (55/52);  new base T (2/5);  total (58/57)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 123 total aligned reads displayed.

ATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGG  >  NC_000913/3308862‑3309152
                                                                                                                                                   |                                                                                                                                               
aTTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCccgccg                                                                                                                                                >  4:438768/1‑149 (MQ=255)
  tACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCcgccgc                                                                                                                                               >  2:251116/1‑148 (MQ=255)
    cTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGc                                                                                                                                            >  3:945839/1‑149 (MQ=255)
      gCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCg                                                                                                                                           >  1:509252/1‑148 (MQ=255)
       cATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCgg                                                                                                                                          <  2:783524/148‑1 (MQ=255)
       cATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGa                                                                                                                                         <  3:724845/149‑1 (MQ=255)
        aTCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAg                                                                                                                                        >  3:513052/1‑149 (MQ=255)
          ccTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAGCGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAgg                                                                                                                                       <  1:536159/148‑1 (MQ=255)
          ccTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCTCCGCTGCGGCgggggg                                                                                                                                      <  4:931412/149‑1 (MQ=255)
            tGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCa                                                                                                                                    >  4:897708/1‑149 (MQ=255)
             gCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCTCCGCTGCGGCGGGGGGctt                                                                                                                                   <  3:736400/149‑3 (MQ=255)
             gCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCTCCGCTGCGGCGGGGGGct                                                                                                                                    <  1:814472/148‑2 (MQ=255)
              cAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCaaa                                                                                                                                  <  3:897708/149‑1 (MQ=255)
               aaGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCaaa                                                                                                                                  <  2:535537/148‑1 (MQ=255)
               aaGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAAt                                                                                                                                 <  4:809383/149‑1 (MQ=255)
               aaGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAAt                                                                                                                                 <  3:231266/149‑1 (MQ=255)
                 gCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCTCCGCTGCGGCGGGGGGcttttaa                                                                                                                               <  4:56741/149‑7 (MQ=255)
                   gTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCa                                                                                                                             <  4:367071/149‑1 (MQ=255)
                   gTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCa                                                                                                                             >  3:248877/1‑149 (MQ=255)
                    tAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCa                                                                                                                             <  1:908496/148‑1 (MQ=255)
                     aaGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCaa                                                                                                                            >  2:474635/1‑148 (MQ=255)
                        tGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCtt                                                                                                                        >  4:295614/1‑149 (MQ=255)
                         gCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCtt                                                                                                                        <  1:160927/148‑1 (MQ=255)
                           tGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACt                                                                                                                     >  3:89997/1‑149 (MQ=255)
                                cAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGccc                                                                                                                >  4:102726/1‑149 (MQ=255)
                                   cGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCATACTCGCCCTg                                                                                                              <  2:466267/148‑1 (MQ=255)
                                     aaaCACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTc                                                                                                           <  4:248877/149‑1 (MQ=255)
                                       acacCAGCGCAAAAAAGGCTTCATTTCCCAGTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTca                                                                                                          <  1:629426/148‑1 (MQ=255)
                                        cacCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTcag                                                                                                         >  2:81457/1‑148 (MQ=255)
                                         acCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTcagc                                                                                                        >  2:733092/1‑148 (MQ=255)
                                          ccAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTcaggag                                                                                                      >  3:635737/1‑149 (MQ=255)
                                             gcgcAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGcc                                                                                                    >  1:372989/1‑148 (MQ=255)
                                             gcgcAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCg                                                                                                   <  4:513052/149‑1 (MQ=255)
                                              cgcAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCg                                                                                                   <  1:764398/148‑1 (MQ=255)
                                               gcAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCgg                                                                                                  <  1:474635/148‑1 (MQ=255)
                                                  aaaaaGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCt                                                                                              >  3:229883/1‑149 (MQ=255)
                                                      aGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTcc                                                                                           >  2:429136/1‑148 (MQ=255)
                                                       ggCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCgg                                                                                         <  3:548565/149‑1 (MQ=255)
                                                         cTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGtg                                                                                       >  4:253211/1‑149 (MQ=255)
                                                         cTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGtg                                                                                       <  3:295614/149‑1 (MQ=255)
                                                          ttCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGtg                                                                                       >  2:466446/1‑148 (MQ=255)
                                                          ttCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGtg                                                                                       >  2:35075/1‑148 (MQ=255)
                                                           tCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGtgct                                                                                     <  3:352256/149‑1 (MQ=255)
                                                             aTTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGtgctgc                                                                                   >  3:911007/1‑149 (MQ=255)
                                                             aTTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGtgctg                                                                                    <  1:81457/148‑1 (MQ=255)
                                                                tCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTgcagc                                                                                >  3:302507/1‑149 (MQ=255)
                                                                 cccACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCt                                                                                                               >  3:175353/1‑117 (MQ=255)
                                                                  ccACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGcagcag                                                                              >  4:390422/1‑149 (MQ=255)
                                                                   cACTCGCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGGCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGcagcagc                                                                             <  4:619544/149‑2 (MQ=255)
                                                                   cACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAgg                                                                             >  3:695471/1‑149 (MQ=255)
                                                                    aCTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAgg                                                                                                                                       <  1:163738/90‑1 (MQ=255)
                                                                    aCTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAgg                                                                                                                                       >  2:163738/1‑90 (MQ=255)
                                                                     cTCCCTAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGt                                                                            >  2:552140/1‑148 (MQ=255)
                                                                      tCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTg                                                                          <  3:390422/149‑1 (MQ=255)
                                                                       cccGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAGAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTga                                                                         <  4:539653/149‑1 (MQ=255)
                                                                       cccGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTga                                                                         <  4:817/149‑1 (MQ=255)
                                                                        ccGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTga                                                                         >  2:414722/1‑148 (MQ=255)
                                                                         cGAAGACCACGGTTGCATGAACGTGCTGTACCCGGGTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCTCCCGCGGCAGCGGAGGTCAAGTGGCAACCTTACTCGCCCTGTTGAGCAGCGGGATCTTCCGGTGCTGCAGCAGGTTgaga                                                                       >  4:799268/1‑149 (MQ=255)
                                                                             gACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTg                                                                    <  1:894123/148‑1 (MQ=255)
                                                                             gACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGc                                                                   <  3:660507/149‑1 (MQ=255)
                                                                                cACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCa                                                                >  3:186524/1‑149 (MQ=255)
                                                                                    gTTGAATGAAGGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCg                                                            >  3:61888/1‑149 (MQ=255)
                                                                                     ttGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTCAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCg                                                            >  2:289098/1‑148 (MQ=255)
                                                                                       gaatgaaCGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGctt                                                         >  4:896308/1‑149 (MQ=255)
                                                                                          tgaaCGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGcttct                                                       >  2:808290/1‑148 (MQ=255)
                                                                                           gaaCGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTa                                                                                                                       <  4:242305/83‑1 (MQ=255)
                                                                                           gaaCGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTa                                                                                                                       >  3:242305/1‑83 (MQ=255)
                                                                                           gaaCGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGcttctt                                                      >  1:695590/1‑148 (MQ=255)
                                                                                             aCGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTa                                                    <  2:386926/148‑1 (MQ=255)
                                                                                                         gggTTGCTAACAAGTCGTCCTGCCCGTTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGACCTTTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCGTCTTTAATGCTCAGACGGa                                       >  4:492326/2‑149 (MQ=255)
                                                                                                         cGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACgg                                        <  2:250173/148‑1 (MQ=255)
                                                                                                               cTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACgg                                  >  1:443116/1‑148 (MQ=255)
                                                                                                                tAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGcc                                >  3:314180/1‑149 (MQ=255)
                                                                                                                tAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGcc                                >  3:186332/1‑149 (MQ=255)
                                                                                                                  aCAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCACAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCt                              >  4:934609/1‑149 (MQ=255)
                                                                                                                     aGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCa                                             <  2:296695/131‑1 (MQ=255)
                                                                                                                     aGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCa                                             >  1:296695/1‑131 (MQ=255)
                                                                                                                     aGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTgg                            <  1:869518/148‑1 (MQ=255)
                                                                                                                     aGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGc                           >  4:767902/1‑149 (MQ=255)
                                                                                                                       gCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCg                          <  1:521527/148‑1 (MQ=255)
                                                                                                                        cGTCCTGCCCGGTTAAAAGCCCCCCGCCACAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGAt                        <  3:934609/149‑1 (MQ=255)
                                                                                                                         gTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGAt                        <  1:35075/148‑1 (MQ=255)
                                                                                                                          tCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCa                      >  4:157633/1‑149 (MQ=255)
                                                                                                                           ccTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCa                      >  1:545117/1‑148 (MQ=255)
                                                                                                                                ccGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTcc                <  3:253211/149‑1 (MQ=255)
                                                                                                                                ccGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTcc                >  3:832953/1‑149 (MQ=255)
                                                                                                                                ccGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTc                 <  2:535361/148‑1 (MQ=255)
                                                                                                                                 cGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCa               >  3:92720/1‑149 (MQ=255)
                                                                                                                                  ggTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGt                                                                                        <  1:90175/75‑1 (MQ=255)
                                                                                                                                  ggTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGt                                                                                        >  2:90175/1‑75 (MQ=255)
                                                                                                                                     tAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGaa            <  1:84590/148‑1 (MQ=255)
                                                                                                                                      aaaaGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACt          >  3:36979/1‑149 (MQ=255)
                                                                                                                                      aaaaGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACt          >  3:497592/1‑149 (MQ=255)
                                                                                                                                       tttgCCCTCCGCTGCGGCGGGGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACtt         >  3:344575/4‑149 (MQ=255)
                                                                                                                                         aGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTc       <  4:229883/149‑1 (MQ=255)
                                                                                                                                          gCCCTCCGCTGCGGCGGGGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCa      >  3:332465/1‑149 (MQ=255)
                                                                                                                                           ccccccGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCa      <  2:830385/148‑1 (MQ=255)
                                                                                                                                             ctccgctgCGGCGGGGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCATTACTTTCAcc    >  2:485917/3‑148 (MQ=255)
                                                                                                                                             ccccGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCAcc    >  1:912680/1‑148 (MQ=255)
                                                                                                                                               ccgccgCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCgg  >  2:652001/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGG  >  NC_000913/3308862‑3309152

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