Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 446,786 C→T 9.9% intergenic (‑120/+29) 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,7860CT 9.9% 191.8 / 7.2 71intergenic (‑120/+29)yajR/cyoEputative major facilitator superfamily transport protein YajR/heme O synthase
Reads supporting (aligned to +/- strand):  ref base C (34/27);  new base T (5/2);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGA  >  NC_000913/446641‑446930
                                                                                                                                                 |                                                                                                                                                
ccTGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAg                                                                                                                                              >  2:151401/1‑150 (MQ=255)
  tGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGcc                                                                                                                                                  <  2:188766/144‑1 (MQ=255)
  tGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGcc                                                                                                                                                  >  1:188766/1‑144 (MQ=255)
             tattatCGTTCATTAGCAGCTCTGTTTTAAATTCGCTTCTAGTTTATAGCGGCAAACGTATGTGGTATACATGCGGTTTTTTTAGTAATTATTTCTTGAGGGATATGTTTCTAATGCTTAGAAAAATGGTGCTGTAGCACCTTTTTAATa                                                                                                                                 <  2:1866996/150‑1 (MQ=255)
             tataATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATa                                                                                                                                 >  1:1911095/1‑150 (MQ=255)
              ataATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATag                                                                                                                                <  2:385255/150‑1 (MQ=255)
              ataATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTCTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATag                                                                                                                                <  2:303624/150‑1 (MQ=255)
                aaTCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATagag                                                                                                                              <  2:968488/150‑1 (MQ=255)
                 aTCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATagaga                                                                                                                             >  1:57809/1‑150 (MQ=255)
                  tCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATagagag                                                                                                                            <  1:1143850/150‑1 (MQ=255)
                   cGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAgg                                                                                                                           >  2:276000/1‑150 (MQ=255)
                      tCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGttt                                                                                                                        <  2:1088625/150‑1 (MQ=255)
                      tCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGttt                                                                                                                        <  2:1187998/150‑1 (MQ=255)
                       cATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGtttt                                                                                                                       >  1:1317629/1‑150 (MQ=255)
                       cATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGtttt                                                                                                                       >  1:593704/1‑150 (MQ=255)
                        aTTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTg                                                                                                                      <  1:1902557/150‑1 (MQ=255)
                         ttACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGt                                                                                                                     <  2:1317629/150‑1 (MQ=255)
                          tACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGtt                                                                                                                    <  2:1822128/150‑1 (MQ=255)
                           accaccTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTa                                                                                                                   <  2:435408/150‑1 (MQ=255)
                            ccaccTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTAc                                                                                                                  >  2:378409/1‑150 (MQ=255)
                             caccTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTAcc                                                                                                                 <  1:151401/150‑1 (MQ=255)
                             caccTATGTTTTAAATTCGCGACAAGTGTAAAGCTGCAAACTAATGGGGTAAACATGCGGATTTTATAGTAAATATTTCGAGAGGTATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTAcc                                                                                                                 <  2:1882017/150‑1 (MQ=255)
                               ccTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACcac                                                                                                               <  2:552133/150‑1 (MQ=255)
                                ctctGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACcaca                                                                                                              >  1:960234/1‑150 (MQ=255)
                                    gTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcagc                                                                                                          >  1:1164911/1‑150 (MQ=255)
                                    gTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcagc                                                                                                          <  2:960234/150‑1 (MQ=255)
                                         aaaTTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcagcagcca                                                                                                     >  1:122134/1‑150 (MQ=255)
                                            ttCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAGTATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcagcagccagca                                                                                                  >  2:1125169/1‑150 (MQ=255)
                                             tCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcagcagccagcag                                                                                                 >  2:2944/1‑150 (MQ=255)
                                               gcgACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCg                                                                                               >  1:1995260/1‑150 (MQ=255)
                                                    tAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGc                                                                                          <  1:510752/150‑1 (MQ=255)
                                                     aGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCg                                                                                         >  2:1133811/1‑150 (MQ=255)
                                                       tgtAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAg                                                                                       >  1:1700025/1‑150 (MQ=255)
                                                       tgtAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAg                                                                                       <  2:248242/150‑1 (MQ=255)
                                                        gtAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGt                                                                                      >  1:230899/1‑150 (MQ=255)
                                                          aaaGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTcc                                                                                    <  2:593704/150‑1 (MQ=255)
                                                          aaaGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCAGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTcc                                                                                    >  1:1872430/1‑150 (MQ=255)
                                                             gCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGt                                                                                 >  2:307754/1‑150 (MQ=255)
                                                             gCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGt                                                                                 <  2:132890/150‑1 (MQ=255)
                                                                          ggggTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCAGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCaa                                                                    >  1:1838889/1‑150 (MQ=255)
                                                                              tAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGa                                                                <  2:122134/150‑1 (MQ=255)
                                                                                     gCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACatcatc                                                         >  1:1864624/1‑150 (MQ=255)
                                                                                       ggATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCAc                                                       >  2:1688297/1‑150 (MQ=255)
                                                                                             tttAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGagag                                                 >  1:97411/1‑150 (MQ=255)
                                                                                                 gTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCa                                             >  1:1237509/1‑150 (MQ=255)
                                                                                                   aaaTATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCAGTAGCACCTTTTTAATAGAGAGGTTTTGTTAACACACAGCAGCAAGCAGCGTATGCGAGTCCGGTACCATAAAATCAAGGACATCATCACCGAGAGGGCAGTg                                          >  1:1655998/1‑150 (MQ=255)
                                                                                                    aaTATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTg                                          >  1:465236/1‑150 (MQ=255)
                                                                                                      tatTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGAt                                        >  1:687633/1‑150 (MQ=255)
                                                                                                      tatTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGAt                                        <  1:307754/150‑1 (MQ=255)
                                                                                                         ttCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGc                                     >  2:232234/1‑150 (MQ=255)
                                                                                                         ttCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGc                                     >  2:225243/1‑150 (MQ=255)
                                                                                                           cGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCga                                   >  1:166350/1‑150 (MQ=255)
                                                                                                            gagaGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCgat                                  <  2:57809/150‑1 (MQ=255)
                                                                                                                gggATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATa                              >  2:637029/1‑150 (MQ=255)
                                                                                                                gggATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATa                              >  1:493685/1‑150 (MQ=255)
                                                                                                                   atatGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCATCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATagaa                           >  2:1703441/1‑149 (MQ=255)
                                                                                                                    tatGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAgaga                          <  2:97411/150‑1 (MQ=255)
                                                                                                                         ttCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCg                     >  1:1227398/1‑150 (MQ=255)
                                                                                                                               tGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGc               <  1:276000/150‑1 (MQ=255)
                                                                                                                                gCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCt              <  2:1700025/150‑1 (MQ=255)
                                                                                                                                 tctattaaaaaGGTGCTACGGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCtt             >  1:294088/7‑150 (MQ=255)
                                                                                                                                 cTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCtt             >  2:508888/1‑150 (MQ=255)
                                                                                                                                 cTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCtt             >  1:161054/1‑150 (MQ=255)
                                                                                                                                  tAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTg            <  2:1237509/150‑1 (MQ=255)
                                                                                                                                  ctattaaaaaGGTGCTACGGCACCTTATTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGATAGTCCGGTACCATAAAATCAACTGACATAATCACCGAGAGGCCAGTGATGGAGATGATAGAGAATCCGAACAGCTTg            >  2:1955060/6‑150 (MQ=255)
                                                                                                                                     ttaaaaaGGTGCTACGGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgc         >  1:1328807/3‑150 (MQ=255)
                                                                                                                                      taaaaaGGTGCTACGGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgcg        >  2:163719/2‑150 (MQ=255)
                                                                                                                                      aaaaaaGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgcg        >  1:1647497/1‑150 (MQ=255)
                                                                                                                                       aaaaaGGTGCTACGGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgcac       >  1:473336/1‑148 (MQ=255)
                                                                                                                                       aaaaaGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgcgc       <  2:1164911/150‑1 (MQ=255)
                                                                                                                                        aaaaGGTGCTGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGcc      <  2:1872430/150‑1 (MQ=255)
                                                                                                                                        aaaaGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGcc      <  2:1864624/150‑1 (MQ=255)
                                                                                                                                            ggTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGa  >  2:1174750/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CCTGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGA  >  NC_000913/446641‑446930

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: