Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 446,787 G→A 6.0% intergenic (‑121/+28) yajR ← / ← cyoE putative major facilitator superfamily transport protein YajR/heme O synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913446,7870GA 6.0% 317.4 / 5.6 100intergenic (‑121/+28)yajR/cyoEputative major facilitator superfamily transport protein YajR/heme O synthase
Reads supporting (aligned to +/- strand):  ref base G (47/47);  new base A (4/2);  total (51/49)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 101 total aligned reads displayed.

CACCTGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGT  >  NC_000913/446639‑446936
                                                                                                                                                    |                                                                                                                                                     
cACCTGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGt                                                                                                                                                      >  2:220805/1‑150 (MQ=255)
  ccTGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAg                                                                                                                                                    <  2:597126/150‑1 (MQ=255)
  ccTGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTAGCGACTAGTGTAAAGCTGCAAACGAATGGGTTAAACATGAGGATTTTTTATTAAATATTTCGAGAGTGATATGTTTCTAATGCTAAGAAAAAAGGTGCAGTAg                                                                                                                                                    <  2:85089/150‑1 (MQ=255)
   attgCGTTATTTTATTATCGTTAATTAAGACCTCTTTTTTTAATTCTCGTATTGTGTAAATGGGCATAGGTATGGGGTAAATATGCTGATTTTTTAGTAAATATTTCTATATTGATATGTTTGTAATTCTAATTAAAATGGTGCCGTAgt                                                                                                                                                   <  2:1254859/147‑2 (MQ=255)
    tGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCa                                                                                                                                                  <  2:798272/150‑1 (MQ=255)
    tGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCa                                                                                                                                                  <  1:1198720/150‑1 (MQ=255)
    tGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCa                                                                                                                                                  >  2:1090551/1‑150 (MQ=255)
     ggCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCAc                                                                                                                                                 >  2:1163650/1‑150 (MQ=255)
     ggCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTAGCGACAAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTATTGCTAAGAAAAAAGGTGCCGTAGCAc                                                                                                                                                 >  2:883495/1‑150 (MQ=255)
       cGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCt                                                                                                                                               >  2:924908/1‑150 (MQ=255)
       cGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCt                                                                                                                                               >  1:1174871/1‑150 (MQ=255)
         tccttttATAAACCTTCATTACCACCTCTGTCTTAAATTCGCCACTAGTGTAAAGCGGCAAACGACTGGGTTAAACAAGCGGATTTTTTAGTAAATAATTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCttt                                                                                                                                             <  2:1098372/147‑1 (MQ=255)
         tCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCttt                                                                                                                                             >  2:1235997/1‑150 (MQ=255)
          cATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCtttt                                                                                                                                            >  1:1299066/1‑150 (MQ=255)
          cATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCtttt                                                                                                                                            <  2:1006660/150‑1 (MQ=255)
           aTTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCTACGGCACCttttt                                                                                                                                           <  1:1433885/150‑1 (MQ=255)
           aTTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCttttt                                                                                                                                           <  2:1671649/150‑1 (MQ=255)
            ttttATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTa                                                                                                                                          >  1:648569/1‑150 (MQ=255)
            ttttATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTa                                                                                                                                          <  2:396940/150‑1 (MQ=255)
             tttATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCTACGGCACCttttttc                                                                                                                                         <  2:525996/150‑3 (MQ=255)
              ttATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAAt                                                                                                                                        >  2:1611322/1‑150 (MQ=255)
               tataATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATa                                                                                                                                       <  1:1090551/150‑1 (MQ=255)
                ataATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATag                                                                                                                                      >  1:1098687/1‑150 (MQ=255)
                  aaTCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATagag                                                                                                                                    <  1:1341602/150‑1 (MQ=255)
                     cGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAgg                                                                                                                                 >  2:102911/1‑150 (MQ=255)
                         cATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAGGTGCAGTAGCACCTTTTTAATAGAGAGGTTTTg                                                                                                                            >  2:832982/1‑150 (MQ=255)
                         cATTACCACATCTGTTTTAAATTCGCGACTACTCTAAACCGGCAAACGAATGGGGTAAACATCCGCATATTTTAGTAAATATTTCGAGAGGTATATGATTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGACAGGtttt                                                                                                                             >  2:582938/1‑150 (MQ=255)
                            tACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGtt                                                                                                                          <  1:220805/150‑1 (MQ=255)
                              ccaccTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTAc                                                                                                                        <  2:130555/150‑1 (MQ=255)
                               caccTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTAcc                                                                                                                       <  1:1335048/150‑1 (MQ=255)
                               caccTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTAcc                                                                                                                       <  1:1037726/150‑1 (MQ=255)
                               caccTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCCGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTAcc                                                                                                                       <  1:1329275/150‑1 (MQ=255)
                                   tctGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACcacac                                                                                                                   >  2:1344720/1‑150 (MQ=255)
                                   tctGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACcacac                                                                                                                   <  2:822304/150‑1 (MQ=255)
                                   tctGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTGTTGTTACcacac                                                                                                                   <  2:1280452/150‑1 (MQ=255)
                                     tGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcag                                                                                                                 <  1:595669/150‑1 (MQ=255)
                                        tttAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcagcag                                                                                                              <  2:440927/150‑1 (MQ=255)
                                         ttAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCAGCTTTTTAATAGAGAGGTTTTGTTACCACAcagcagc                                                                                                             >  2:105536/1‑150 (MQ=255)
                                             aTTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcagcagccagc                                                                                                         >  1:712613/1‑150 (MQ=255)
                                               tCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcagcagccagcag                                                                                                       >  1:110207/1‑150 (MQ=255)
                                                 gcgACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCg                                                                                                     <  1:1029777/150‑1 (MQ=255)
                                                   gACTAGTGTAAATCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcagcag                                                                                                              <  2:402852/139‑1 (MQ=255)
                                                   gACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTa                                                                                                   <  1:102911/150‑1 (MQ=255)
                                                   gACTAGTGTAAAACGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcagcag                                                                                                              >  1:402852/1‑139 (MQ=255)
                                                       aGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCg                                                                                               <  1:467958/150‑1 (MQ=255)
                                                       aGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCg                                                                                               <  1:1235997/150‑1 (MQ=255)
                                                         tgtAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAg                                                                                             <  2:1461223/150‑1 (MQ=255)
                                                          gtAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTg                                                                                                                            >  2:12310/1‑118 (MQ=255)
                                                          gtAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTg                                                                                                                            <  1:12310/118‑1 (MQ=255)
                                                          gtAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGt                                                                                            <  2:1174871/150‑1 (MQ=255)
                                                              aGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCgg                                                                                        <  1:1611322/150‑1 (MQ=255)
                                                               gCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCg                                                                                               <  2:935888/142‑1 (MQ=255)
                                                               gCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCg                                                                                               >  1:935888/1‑142 (MQ=255)
                                                                cGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTa                                                                                      <  1:455752/150‑1 (MQ=255)
                                                                 ggCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACAACACAGCAGCCAGCAGCGTATGCGAGTCCGGTAc                                                                                     >  1:580640/1‑150 (MQ=255)
                                                                   cAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATTTTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTCCCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCa                                                                                   >  2:1287436/1‑150 (MQ=255)
                                                                   cAAACGAATGGGGTAAACATGCGGATTTTTTAGAAAATATTTCCAGAGCGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTATTTAATAGAGAGGATTTGTTACCACACAGCAGCCAGCAGCCTATGCGAGTCAGGTACCa                                                                                   >  2:1222091/1‑150 (MQ=255)
                                                                     aaCGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGATAAGAAAAAAGGTGCCGTAGCCCGTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATa                                                                                 >  2:1179626/1‑150 (MQ=255)
                                                                          aTGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATc                                                                            <  1:541535/150‑1 (MQ=255)
                                                                             gggTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAAc                                                                         <  2:752624/150‑1 (MQ=255)
                                                                               gTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACgg                                                                       <  1:430166/150‑1 (MQ=255)
                                                                                tAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGa                                                                      >  1:59297/1‑150 (MQ=255)
                                                                                           aTTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCg                                                           >  1:1257681/1‑150 (MQ=255)
                                                                                             tttttAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAATAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTAGCATAAAATCAACGGACATCATCACCgag                                                         <  2:1303251/150‑1 (MQ=255)
                                                                                               tttAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGagag                                                       >  2:1668992/1‑150 (MQ=255)
                                                                                                     aaaTATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGt                                                 <  2:1098687/150‑1 (MQ=255)
                                                                                                        tatTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGAt                                              >  2:546710/1‑150 (MQ=255)
                                                                                                            tCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGGGGGCAGTGATGGCg                                          <  1:1344720/150‑1 (MQ=255)
                                                                                                             cGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACTGACATCATCACCGAGAGGGCATTGATgtaga                                         <  2:1600450/150‑5 (MQ=255)
                                                                                                                   ggATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATag                                   <  2:712613/150‑1 (MQ=255)
                                                                                                                    gATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATaga                                  >  1:999751/1‑150 (MQ=255)
                                                                                                                     atatGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATagag                                 <  1:546710/150‑1 (MQ=255)
                                                                                                                      tatGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAgaga                                >  2:272954/1‑150 (MQ=255)
                                                                                                                        tGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTAcc                                                                                    <  1:211714/96‑1 (MQ=255)
                                                                                                                        tGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTAcc                                                                                    >  2:211714/1‑96 (MQ=255)
                                                                                                                        tGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAg                              >  2:27483/1‑150 (MQ=255)
                                                                                                                           ttCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCg                           >  2:197865/1‑150 (MQ=255)
                                                                                                                            tCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTAGCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGa                          >  2:1366886/1‑150 (MQ=255)
                                                                                                                              tAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAAc                        >  1:131649/1‑150 (MQ=255)
                                                                                                                               aaTGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGagag                                                       <  1:883548/118‑1 (MQ=255)
                                                                                                                               aaTGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGagag                                                       <  1:1372465/118‑1 (MQ=255)
                                                                                                                               aaTGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGagag                                                       >  2:883548/1‑118 (MQ=255)
                                                                                                                               aaTGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGagag                                                       >  2:1372465/1‑118 (MQ=255)
                                                                                                                                aTGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAg                      >  1:319785/1‑150 (MQ=255)
                                                                                                                                  ctctattaaaaaGGTGCTACGGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCt                    >  1:1605172/8‑150 (MQ=255)
                                                                                                                                   cTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCtt                   >  2:233948/1‑150 (MQ=255)
                                                                                                                                       gAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgc               >  1:1035000/1‑150 (MQ=255)
                                                                                                                                       gAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgc               >  1:604227/1‑150 (MQ=255)
                                                                                                                                         aaaaaGGTGCTACGGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGCCATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgcgc             >  1:895071/1‑150 (MQ=255)
                                                                                                                                         aaaaaGGTGCTACGGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgcgc             >  1:175676/1‑150 (MQ=255)
                                                                                                                                          aaaaGGTGCTACGGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGcc            >  1:1371181/1‑150 (MQ=255)
                                                                                                                                          aaaaGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCACCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACCGCTTGAGCGcc            >  2:1264850/1‑150 (MQ=255)
                                                                                                                                             aGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAg         <  1:1163650/150‑1 (MQ=255)
                                                                                                                                                tGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCgatg                                       >  2:1135305/1‑117 (MQ=255)
                                                                                                                                                tGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCgatg                                       <  1:1135305/117‑1 (MQ=255)
                                                                                                                                                 gCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTc     >  2:1645915/1‑150 (MQ=255)
                                                                                                                                                  ccGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCt    >  1:1558640/1‑150 (MQ=255)
                                                                                                                                                  ccGTAGCACCTATTTAATAGAGAGGTTTTCTTACCCCACAGCAGCCAGCAGCGTATGCGAGTCCGCTACCATAAAATCAACGGACATCATCAACGAGAGGGCAGAGATGGCGATGATAGAGAAGACGAACAGCTTGCGCGCCCAGATTCt    >  2:1232919/1‑150 (MQ=255)
                                                                                                                                                   cGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTg   <  1:321944/150‑1 (MQ=255)
                                                                                                                                                    gTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTAACATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAACCAGAACAGCTTGCGCGCCCAGATTCTGt  >  1:591912/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CACCTGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGT  >  NC_000913/446639‑446936

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