Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 465,603 A→G 5.0% intergenic (‑56/+9) queC ← / ← ybaE 7‑cyano‑7‑deazaguanine synthase/putative protein YbaE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913465,6030AG 5.0% 246.4 / 2.5 79intergenic (‑56/+9)queC/ybaE7‑cyano‑7‑deazaguanine synthase/putative protein YbaE
Reads supporting (aligned to +/- strand):  ref base A (42/33);  new base G (2/2);  total (44/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CACGCAATGGACTTCATCATATTGTTGTAATGCCTGCACCAGACAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGT  >  NC_000913/465455‑465746
                                                                                                                                                    |                                                                                                                                               
cACGCAATGGACTTCATCATATTGTTGTAATGCCTGCACCAGACAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGAc                                                                                                                                                <  1:1015402/150‑1 (MQ=255)
    cAATGGACTTCATCATATTGTTGTAATGCCTGCACCAGACAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGc                                                                                                                                            <  2:1490573/150‑1 (MQ=255)
          aCTTCATCATATTGTTGTAATGCCTGCACCAGACAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCa                                                                                                                                      >  1:884089/1‑150 (MQ=255)
              catcatATTGTTGTAATGCCTGCACCAGACAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGc                                                                                                                                  <  1:1641514/150‑1 (MQ=255)
                 catATTGTTGTAATGCCTGCACCAGACAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGAt                                                                                                                               <  2:1335868/150‑1 (MQ=255)
                      tgttgtAATGCCTGCACCAGACAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGgcgc                                                                                                                          >  1:917369/1‑150 (MQ=255)
                       gttgtAATGCCTGCACCAGACAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGGCCAGCTATTCATTGCGATGGcgcg                                                                                                                         >  2:38576/1‑150 (MQ=255)
                       gttgtAATGCCTGCACCAGACAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGGCCAGCTATTCATTGCGATGGcgcg                                                                                                                         >  2:1039984/1‑150 (MQ=255)
                         tgtAATGCCTGCACCAGACAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCggg                                                                                                                       <  1:1547026/150‑1 (MQ=255)
                            aaTGCCTGCACCAGACAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAg                                                                                                                    >  2:546698/1‑150 (MQ=255)
                            aaTGCCTGCACCAGACAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAg                                                                                                                    <  1:1253672/150‑1 (MQ=255)
                              tGCCTGCACCAGACAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGcc                                                                                                                  >  2:419465/1‑150 (MQ=255)
                                ccTGCACCAGACAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAg                                                                                                                >  2:1857203/1‑150 (MQ=255)
                                     aCCAGACAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTc                                                                                                           <  1:85229/150‑1 (MQ=255)
                                      ccAGACAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCg                                                                                                          >  2:376144/1‑150 (MQ=255)
                                           cAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTaaa                                                                                                     >  2:1955369/1‑150 (MQ=255)
                                                  tggAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAAc                                                                                              >  1:566291/1‑150 (MQ=255)
                                                   ggAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAAcc                                                                                             <  1:1857203/150‑1 (MQ=255)
                                                      aTCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGc                                                                                          <  1:389110/150‑1 (MQ=255)
                                                         cTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCcgc                                                                                       >  2:1221203/1‑150 (MQ=255)
                                                         cTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCcgc                                                                                       >  2:964011/1‑150 (MQ=255)
                                                             ccTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCgg                                                                                   >  2:1318360/1‑150 (MQ=255)
                                                             ccTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCgg                                                                                   >  2:1593373/1‑150 (MQ=255)
                                                             ccTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCgg                                                                                   >  2:1654065/1‑150 (MQ=255)
                                                              cTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCGTTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGgt                                                                                  >  2:434288/1‑150 (MQ=255)
                                                              cTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGgt                                                                                  >  2:1998301/1‑150 (MQ=255)
                                                               tCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGgtg                                                                                 >  2:1319848/1‑150 (MQ=255)
                                                                     gAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTc                                                                           >  1:1714959/1‑150 (MQ=255)
                                                                       acacAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGc                                                                         <  2:273640/150‑1 (MQ=255)
                                                                       acacAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGc                                                                         >  2:283822/1‑150 (MQ=255)
                                                                       acacAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCACGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGc                                                                         >  1:358492/1‑150 (MQ=255)
                                                                        cacaACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGca                                                                        >  2:936270/1‑150 (MQ=255)
                                                                        cacaACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGca                                                                        >  2:943197/1‑150 (MQ=255)
                                                                         acaACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGcac                                                                       <  1:1221203/150‑1 (MQ=255)
                                                                          caACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGcaca                                                                      >  1:1035225/1‑150 (MQ=255)
                                                                          caACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGcaca                                                                      >  2:1370589/1‑150 (MQ=255)
                                                                                 aGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCa                                                               <  2:136608/150‑1 (MQ=255)
                                                                                  gCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCAc                                                              >  2:1558013/1‑150 (MQ=255)
                                                                                   cACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACg                                                             >  2:497399/1‑150 (MQ=255)
                                                                                    aCGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGc                                                            >  2:1158473/1‑150 (MQ=255)
                                                                                     cGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGcc                                                           <  1:419465/150‑1 (MQ=255)
                                                                                             aGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCACACACCGTTCACGCCGGAGGGGc                                                  <  1:1000644/150‑1 (MQ=255)
                                                                                             aGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCACACACCGTTCACGCCGGAGGGGc                                                  <  1:1007502/150‑1 (MQ=255)
                                                                                                 aaTCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCGTTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACt                                               <  1:434288/150‑1 (MQ=255)
                                                                                                 aaTCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACt                                               <  1:81219/150‑1 (MQ=255)
                                                                                                 aaTCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACt                                               >  1:1514469/1‑150 (MQ=255)
                                                                                                       gCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCg                                         >  2:1944077/1‑150 (MQ=255)
                                                                                                        cTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCga                                        >  2:755682/1‑150 (MQ=255)
                                                                                                           gAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCgatag                                     >  2:1579652/1‑150 (MQ=255)
                                                                                                            aaCAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCgatagt                                    <  1:943197/150‑1 (MQ=255)
                                                                                                            aaCAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCgatagt                                    <  1:936270/150‑1 (MQ=255)
                                                                                                             aCAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCgatagtg                                   >  1:1993174/1‑150 (MQ=255)
                                                                                                              cAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCCGGAGGGGCACTAATGCgatagtga                                  >  1:744338/1‑150 (MQ=255)
                                                                                                                aTAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCgatagtgata                                <  1:1158473/150‑1 (MQ=255)
                                                                                                                 tAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCgatagtgatag                               <  2:382340/150‑1 (MQ=255)
                                                                                                                  aaaaGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCgatagtgatagt                              >  2:609033/1‑150 (MQ=255)
                                                                                                                   aaaGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGtt                             >  2:1743381/1‑150 (MQ=255)
                                                                                                                     aGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGa                           <  2:884089/150‑1 (MQ=255)
                                                                                                                     aGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGa                           <  1:1061439/150‑1 (MQ=255)
                                                                                                                      gCGTTATTGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCTCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGaa                          <  2:1496469/150‑1 (MQ=255)
                                                                                                                            tGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGATGGGCACTAATGCGATAGTGATAGTTGAACAgcgg                    <  1:1370589/150‑1 (MQ=255)
                                                                                                                            tGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAgcgg                    <  2:1714959/150‑1 (MQ=255)
                                                                                                                               tAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGt                 >  1:211873/1‑150 (MQ=255)
                                                                                                                               tAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGt                 >  1:1030828/1‑150 (MQ=255)
                                                                                                                                aaCGCCTGCGATTAACCCGGGCCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTc                <  1:1039984/150‑1 (MQ=255)
                                                                                                                                aaCGCCTGCGATTAACCCGGGCCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTc                <  1:38576/150‑1 (MQ=255)
                                                                                                                                  cGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAg              <  2:358492/150‑1 (MQ=255)
                                                                                                                                    ccTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCg            >  1:769821/1‑150 (MQ=255)
                                                                                                                                       gCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCg         <  2:565366/150‑1 (MQ=255)
                                                                                                                                       gCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCg         >  1:1707439/1‑150 (MQ=255)
                                                                                                                                       gCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCg         >  2:1946773/1‑150 (MQ=255)
                                                                                                                                         gATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCa       <  1:1319848/150‑1 (MQ=255)
                                                                                                                                          aTTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCAt      <  2:1993174/150‑1 (MQ=255)
                                                                                                                                           ttAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATc     >  1:1593533/1‑150 (MQ=255)
                                                                                                                                            tAACCCGGACCCGCTATTCATTGCGATGGCGCGGGAAGCCATGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCt    <  1:470811/150‑1 (MQ=255)
                                                                                                                                             aaCCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCGGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTg   <  1:1998301/150‑1 (MQ=255)
                                                                                                                                              aCCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTc                >  1:1305444/1‑136 (MQ=255)
                                                                                                                                              aCCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTc                <  2:1305444/136‑1 (MQ=255)
                                                                                                                                              aCCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGt  >  1:1757527/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CACGCAATGGACTTCATCATATTGTTGTAATGCCTGCACCAGACAGGTGGTGGAATCCTGACCTCCACTGAACACAACGACAGCACGTTTCATAGATAATCCTGCTTGAACAATAAAAGCGTTATGGTAACGCCTGCGATTAACCCGGACCAGCTATTCATTGCGATGGCGCGGGAAGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGT  >  NC_000913/465455‑465746

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