Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,788,674 A→T 6.5% intergenic (+25/‑311) ygaQ → / → csiD putative uncharacterized protein YgaQ/PF08943 family protein CsiD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,788,6740AT 6.5% 474.5 / 13.4 149intergenic (+25/‑311)ygaQ/csiDputative uncharacterized protein YgaQ/PF08943 family protein CsiD
Reads supporting (aligned to +/- strand):  ref base A (69/70);  new base T (2/8);  total (71/78)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

AAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCTTA  >  NC_000913/2788528‑2788823
                                                                                                                                                  |                                                                                                                                                     
aaaTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAAc                                                                                                                                                    <  2:1259854/150‑1 (MQ=255)
aaaTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAAc                                                                                                                                                    <  2:824305/150‑1 (MQ=255)
   ttATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACggg                                                                                                                                                 <  1:620826/150‑1 (MQ=255)
   ttATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACggg                                                                                                                                                 <  2:731433/150‑1 (MQ=255)
    tATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTAATGTTATATCTGCCCCGTTTTATCGGGGc                                                                                                                                               <  1:850016/150‑1 (MQ=255)
     aTCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGc                                                                                                                                               <  1:1123113/150‑1 (MQ=255)
     aTCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGc                                                                                                                                               <  1:1279332/150‑1 (MQ=255)
       ccccccTTGTAAAGTACTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCAg                                                                                                                                             <  1:589555/150‑1 (MQ=255)
        cccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCAGa                                                                                                                                            <  2:974113/150‑1 (MQ=255)
        cccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCAGa                                                                                                                                            <  1:879855/150‑1 (MQ=255)
          cccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGata                                                                                                                                          >  2:1298995/1‑150 (MQ=255)
          cccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGGAGCGCTCCTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCAGata                                                                                                                                          <  1:380642/150‑1 (MQ=255)
            cTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCAGatata                                                                                                                                        <  1:31840/150‑3 (MQ=255)
            cTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCAGatata                                                                                                                                        <  1:26104/150‑3 (MQ=255)
            cTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGataat                                                                                                                                        >  2:1529386/1‑150 (MQ=255)
                tAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgt                                                                                                                                    >  2:510186/1‑150 (MQ=255)
                      aTTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtttagtt                                                                                                                              <  2:1069075/150‑1 (MQ=255)
                         tACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACATGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGTGCAGATAATATGTTTAGTTTAc                                                                                                                           <  2:915418/150‑1 (MQ=255)
                            gCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTaa                                                                                                                        >  1:1288233/1‑150 (MQ=255)
                             cTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAAc                                                                                                                       >  2:1308679/1‑150 (MQ=255)
                             cTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAAc                                                                                                                       >  2:812407/1‑150 (MQ=255)
                                 gagGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTc                                                                                                                   >  2:41989/1‑150 (MQ=255)
                                 gagGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTc                                                                                                                   >  2:140039/1‑150 (MQ=255)
                                   ggCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCAt                                                                                                                 <  1:551173/150‑1 (MQ=255)
                                   ggCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTAAACTaaacatatt                                                                                                                 <  1:1094627/150‑8 (MQ=255)
                                       taaCTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTAGTTTACTAACGGTCATTTTGc                                                                                                            <  2:750480/149‑1 (MQ=255)
                                         aCTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCa                                                                                                           <  2:435038/150‑1 (MQ=255)
                                         aCTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCa                                                                                                           <  1:1331679/150‑1 (MQ=255)
                                          ctctAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAg                                                                                                          <  2:1192516/150‑1 (MQ=255)
                                          ctctAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAg                                                                                                          >  1:1244308/1‑150 (MQ=255)
                                           tctAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAAACGCCGAAATATCTCTATTTATAATTTAATGTTATATCTTCAAAGATAAAACGTGGCAGATAATATGTTTAGTTTACTAAAGGTAATTTTGCAGt                                                                                                         >  2:1100041/1‑150 (MQ=255)
                                                 tatCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGcc                                                                                                   >  2:1409509/1‑150 (MQ=255)
                                                    cAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCAtt                                                                                                >  2:811008/1‑150 (MQ=255)
                                                     aaTATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAATGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCAttt                                                                                               >  2:422275/1‑150 (MQ=255)
                                                     aaTATTCCTGGGCGTAGCGCTCTTCTTTAAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCAGATAAAACGGGGCAGATAATATGTTTACTTTACTAACGGTCATTTTGCAGTGAAGCCAttt                                                                                               >  2:1002338/1‑150 (MQ=255)
                                                       tatTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTAc                                                                                             <  2:458193/150‑1 (MQ=255)
                                                           ccTGGGCGTAGCGCCCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGtt                                                                                         >  1:1038989/1‑150 (MQ=255)
                                                             tGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGtttt                                                                                       <  2:1119907/150‑1 (MQ=255)
                                                                  gTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATc                                                                                  <  2:1130879/150‑1 (MQ=255)
                                                                    aGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGa                                                                                >  1:261720/1‑150 (MQ=255)
                                                                     gcgcTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGaa                                                                               >  1:1512934/1‑149 (MQ=255)
                                                                     gcgcTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGAc                                                                               >  2:1340815/1‑150 (MQ=255)
                                                                      cgcTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAATTTGTTATATAATATTTTTAGTTTACTTAGTGTTATTTTTTAGTTAAGTCATTTACTTTTTTTTTTCGAcc                                                                              <  2:344776/150‑1 (MQ=255)
                                                                       gcTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCa                                                                             <  2:298192/150‑1 (MQ=255)
                                                                                   gggAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtc                                                                 >  1:145540/1‑150 (MQ=255)
                                                                                   gggAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtc                                                                 <  1:201240/150‑1 (MQ=255)
                                                                                    ggAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAATGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtct                                                                <  1:422275/150‑1 (MQ=255)
                                                                                      aaaaCATGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTctct                                                              <  2:145540/150‑1 (MQ=255)
                                                                                       aaaCAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTa                                                             <  2:919422/150‑1 (MQ=255)
                                                                                       aaaCAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACAAGATAATCTGTTCTCTa                                                             >  2:1071885/1‑150 (MQ=255)
                                                                                         aCAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATTTGTTCTCTAAt                                                           >  1:609990/1‑150 (MQ=255)
                                                                                         aCAGGAGAACCGCCGAACTATATCTATCTATAATTTAATGTTATATCTTCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAAt                                                           <  1:1400085/150‑1 (MQ=255)
                                                                                          cAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATg                                                          <  1:1173502/150‑1 (MQ=255)
                                                                                               gaACCGCCGACCTATCTCTATTTAAAATTTAATGTTATATCTGCCCCGATAAACCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCCGTGAAGCCCTTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAAc                                                     <  2:1501690/150‑1 (MQ=255)
                                                                                                   cgccgAACTATGTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTccc                                                 <  2:1269218/150‑1 (MQ=255)
                                                                                                   cgccgAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCTGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTccc                                                 <  1:1494037/150‑1 (MQ=255)
                                                                                                    gccgAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTcccc                                                <  1:343143/150‑1 (MQ=255)
                                                                                                    gccgAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAACCCATTTAATGTTTTTTATCGACAAGATAATATGTTATCTAATGTTAACTcccc                                                >  2:1233524/1‑150 (MQ=255)
                                                                                                     ccgAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTccccc                                               >  2:184019/1‑150 (MQ=255)
                                                                                                     ccgAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTccccc                                               <  2:159633/150‑1 (MQ=255)
                                                                                                     ccgAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTccccc                                               <  2:825808/150‑1 (MQ=255)
                                                                                                          cTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAAcc                                          <  1:761208/150‑1 (MQ=255)
                                                                                                              ctctATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTTTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGtt                                      >  2:183959/1‑150 (MQ=255)
                                                                                                               tctATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTg                                     >  2:1425120/1‑150 (MQ=255)
                                                                                                                ctATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGc                                    <  1:1060216/150‑1 (MQ=255)
                                                                                                                ctATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGc                                    >  1:164377/1‑150 (MQ=255)
                                                                                                                 tATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCt                                   >  2:370323/1‑150 (MQ=255)
                                                                                                                 tATTTATAATTTAATGTTATATATGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTAATAACGGTCATTTTGCAGTGAAGCAATTTAATGTTTTTTATCGCCAAGATAATCTGTTCTCTAATGTTAACTCCCACTAAAATGTTGCt                                   >  2:1458252/1‑150 (MQ=255)
                                                                                                                    ttATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAAAGGTCATTTTGCAGTGAAGACATTTACTGTTTTTTATCGACAAGATAATATGTTCTAAAATGTTAAATACACCTAACCTGTTGCTTTa                                >  2:1394120/1‑150 (MQ=255)
                                                                                                                     tataATTTAATGTTTTATCTGACTGTATAAAACGGTTCATATAATATGTTTATTTTACTAACGTTCATTTTGTATTGAATCCATTTACTGTTTTTTATTGACAAGATAATCTGTTCTCTAATGTTAACTACCCCTAACCTGTTGGTTTAg                               <  2:1512934/150‑1 (MQ=255)
                                                                                                                     tataATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACt                                                    <  2:1459994/129‑1 (MQ=255)
                                                                                                                     tataATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAg                               <  1:928247/150‑1 (MQ=255)
                                                                                                                      ataATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGt                              >  1:1122113/1‑150 (MQ=255)
                                                                                                                           ttAATTTTTTATCTGTCCCGATAAAATGGGTCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCTATTTACTGTTTTTTATCGACAAGTTAATCTGTTCTCTAATGTTATCTCCCCCTAACCTGTTGCTTTAGTTATTc                         <  2:1249648/150‑1 (MQ=255)
                                                                                                                           ttAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTc                         <  1:648459/150‑1 (MQ=255)
                                                                                                                           ttAATGTTATATCTGCCCCGATAAAACGGGGAAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTc                         >  2:122882/1‑150 (MQ=255)
                                                                                                                            tAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCa                        <  1:811008/150‑1 (MQ=255)
                                                                                                                                gTTATATCTGCACCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTc                    >  2:1251105/1‑150 (MQ=255)
                                                                                                                                 ttATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTcc                   >  1:155521/1‑150 (MQ=255)
                                                                                                                                  tataTCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCt                  >  1:1230881/1‑150 (MQ=255)
                                                                                                                                   ttatCTGCCCCGTTTTATCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTg                 >  1:506658/2‑150 (MQ=255)
                                                                                                                                   atatCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTg                 <  1:108959/150‑1 (MQ=255)
                                                                                                                                    tatCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGt                <  1:752054/150‑1 (MQ=255)
                                                                                                                                      tCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtct              >  1:471858/1‑150 (MQ=255)
                                                                                                                                      tCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtct              <  2:323397/150‑1 (MQ=255)
                                                                                                                                         gCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGtt                             >  2:30564/1‑132 (MQ=255)
                                                                                                                                         gCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGtt                             <  1:30564/132‑1 (MQ=255)
                                                                                                                                         gCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTAtt                          >  1:191615/1‑135 (MQ=255)
                                                                                                                                         gCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTAtt                          <  2:191615/135‑1 (MQ=255)
                                                                                                                                         gCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAAATCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCaa           >  2:1033667/1‑149 (MQ=255)
                                                                                                                                          ccccGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACt          >  2:1151636/1‑150 (MQ=255)
                                                                                                                                           cccGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACtt         <  1:1529386/150‑1 (MQ=255)
                                                                                                                                             cGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTg       <  1:184019/150‑1 (MQ=255)
                                                                                                                                              gATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTg                 <  1:758022/139‑1 (MQ=255)
                                                                                                                                              gATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTg                 >  2:758022/1‑139 (MQ=255)
                                                                                                                                              gATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTCCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGc      <  2:508422/150‑1 (MQ=255)
                                                                                                                                               aTAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGcc     <  1:506629/150‑1 (MQ=255)
                                                                                                                                                tAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCt    <  2:755223/150‑1 (MQ=255)
                                                                                                                                                  aaaCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCTTa  >  1:143645/1‑150 (MQ=255)
                                                                                                                                                  aaaCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCTTa  >  1:1213609/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCTTA  >  NC_000913/2788528‑2788823

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: