Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 886,196 C→T 6.1% A382T (GCT→ACT) ‡ ybjJ ← inner membrane protein YbjJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913886,1960CT 6.1% 312.4 / 4.7 98A382T (GCT→ACT) ‡ybjJinner membrane protein YbjJ
Reads supporting (aligned to +/- strand):  ref base C (44/48);  new base T (3/3);  total (47/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCATAAACCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTGCATCGGGGCCG  >  NC_000913/886048‑886344
                                                                                                                                                    |                                                                                                                                                    
aGCCATAAACCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGcc                                                                                                                                                     <  1:1262040/150‑1 (MQ=255)
 gCCATAAACCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCcg                                                                                                                                                    >  2:471250/1‑150 (MQ=255)
     tAAACCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGGCTATTGTCGCgaa                                                                                                                                                <  1:1228847/150‑2 (MQ=255)
     tAAACCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCgag                                                                                                                                                <  2:593529/150‑1 (MQ=255)
       aaCCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGGCTATTGTCGCgaaag                                                                                                                                              <  1:297733/150‑4 (MQ=255)
       aaCCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGaa                                                                                                                                              >  2:1536665/1‑150 (MQ=255)
        aCCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGGCTATTGTCGCgaaagc                                                                                                                                             <  1:1503028/150‑5 (MQ=255)
         ccGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATa                                                                                                                                            >  1:42525/1‑150 (MQ=255)
         ccGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATa                                                                                                                                            >  1:36439/1‑150 (MQ=255)
            cTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAAcc                                                                                                                                         >  1:1542138/1‑150 (MQ=255)
              cACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATACATGTATTCCGCAATTAATTTAATGCTCATCAACTATTCTCGATAGCCGTCTGAGTTTTTGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAg                                                                                                                                       >  2:559194/1‑150 (MQ=255)
               aCGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAgc                                                                                                                                      <  2:110283/150‑1 (MQ=255)
                  gTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGcc                                                                                                                                   >  1:703915/1‑150 (MQ=255)
                     ataGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGt                                                                                                                                <  2:234407/150‑1 (MQ=255)
                     ataGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGt                                                                                                                                <  1:706697/150‑1 (MQ=255)
                        ggTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACa                                                                                                                             >  1:1240880/1‑150 (MQ=255)
                          ttttttGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAAc                                                                                                                           >  1:696987/1‑150 (MQ=255)
                          ttttttGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAAc                                                                                                                           <  1:1536665/150‑1 (MQ=255)
                           tttttGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAAcc                                                                                                                          <  2:696987/150‑1 (MQ=255)
                              ttGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATGACCGGCGCCAGTGCAGccgcc                                                                                                                       >  1:687980/1‑147 (MQ=255)
                              ttGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGc                                                                                                                       >  1:1605912/1‑150 (MQ=255)
                                gATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCAt                                                                                                                     <  1:1319617/150‑1 (MQ=255)
                                  tCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTg                                                                                                                   <  1:91629/150‑1 (MQ=255)
                                    gCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCa                                                                                                                 >  2:24967/1‑150 (MQ=255)
                                    gCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCa                                                                                                                 >  1:235783/1‑150 (MQ=255)
                                       tAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTa                                                                                                              >  1:1570911/1‑150 (MQ=255)
                                       tAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTa                                                                                                              >  2:1174752/1‑150 (MQ=255)
                                         aGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACg                                                                                                            >  2:1047170/1‑150 (MQ=255)
                                         aGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACg                                                                                                            <  2:808283/150‑1 (MQ=255)
                                              gTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATc                                                                                                       >  2:302109/1‑150 (MQ=255)
                                              gTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATc                                                                                                       <  2:806789/150‑1 (MQ=255)
                                               tACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATcc                                                                                                      <  2:456547/150‑1 (MQ=255)
                                                   atccatGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATaa                                                                                                  <  1:24967/150‑1 (MQ=255)
                                                    tccatGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTAc                                                                                                             <  2:680812/138‑1 (MQ=255)
                                                    tccatGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTAc                                                                                                             >  1:680812/1‑138 (MQ=255)
                                                     ccatGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGATATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATg                                                                                                >  1:575090/1‑150 (MQ=255)
                                                      catGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGt                                                                                               >  1:1479988/1‑150 (MQ=255)
                                                       atGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGtt                                                                                              <  1:1143274/150‑1 (MQ=255)
                                                          tCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGc                                                                                           >  1:578201/1‑150 (MQ=255)
                                                           cTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGcc                                                                                          <  1:1107402/150‑1 (MQ=255)
                                                           cTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGcc                                                                                          <  2:1605912/150‑1 (MQ=255)
                                                            tACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCg                                                                                         >  1:502453/1‑150 (MQ=255)
                                                            tACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCg                                                                                         <  2:1249538/150‑1 (MQ=255)
                                                             aCCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGTGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCga                                                                                        <  1:1472723/150‑1 (MQ=255)
                                                             aCCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCga                                                                                        <  1:48720/150‑1 (MQ=255)
                                                             aCCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCga                                                                                        <  1:836051/150‑1 (MQ=255)
                                                                 cAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATa                                                                                    <  2:920263/150‑1 (MQ=255)
                                                                  aattaattTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATac                                                                                   >  2:85495/1‑149 (MQ=255)
                                                                      aatttaatGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCca                                                                               <  2:1542138/150‑1 (MQ=255)
                                                                       atttaatGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCcag                                                                              >  1:518331/1‑150 (MQ=255)
                                                                       atttaatGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCcag                                                                              >  2:977620/1‑150 (MQ=255)
                                                                       atttaatGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCcag                                                                              >  2:1056153/1‑150 (MQ=255)
                                                                          taatGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCcagcag                                                                           >  1:482197/1‑150 (MQ=255)
                                                                            atGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAgcg                                                                         <  2:104075/150‑1 (MQ=255)
                                                                               cTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAgcggcg                                                                      >  1:460355/1‑150 (MQ=255)
                                                                               cTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAgcggcg                                                                      <  1:543627/150‑1 (MQ=255)
                                                                                 catcaACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCggcggc                                                                    <  1:1340077/150‑1 (MQ=255)
                                                                                 catcaACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCggcggc                                                                    >  2:1069984/1‑150 (MQ=255)
                                                                                  atcaACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGcc                                                                   <  2:1286524/150‑1 (MQ=255)
                                                                                    caACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCcg                                                                 >  2:243355/1‑150 (MQ=255)
                                                                                    caACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCcg                                                                 >  2:858822/1‑150 (MQ=255)
                                                                                       cTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCcgacg                                                              <  1:243355/150‑1 (MQ=255)
                                                                                       cTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCcgacg                                                              >  2:1471531/1‑150 (MQ=255)
                                                                                          ttCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAgg                                                           <  2:1284750/150‑1 (MQ=255)
                                                                                              ccATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAg                                                       <  1:1366387/150‑1 (MQ=255)
                                                                                                 tCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCa                                                    >  2:124481/1‑150 (MQ=255)
                                                                                                     cGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATa                                                <  2:719085/150‑1 (MQ=255)
                                                                                                              tttGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGt                                       <  2:1570911/150‑1 (MQ=255)
                                                                                                                  gTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCt                                   <  2:1240880/150‑1 (MQ=255)
                                                                                                                   tATCGGGTTTGGCGACGGTTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTa                                  >  2:933008/1‑150 (MQ=255)
                                                                                                                   tATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTa                                  >  1:1539494/1‑150 (MQ=255)
                                                                                                                    aTCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTAc                                 >  1:445317/1‑150 (MQ=255)
                                                                                                                    aTCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTAc                                 <  2:578201/150‑1 (MQ=255)
                                                                                                                     tCGTGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTAcc                                <  1:350029/150‑1 (MQ=255)
                                                                                                                     tCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTAcc                                >  1:188001/1‑150 (MQ=255)
                                                                                                                        ggTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACcaca                             <  1:1047170/150‑1 (MQ=255)
                                                                                                                             ggCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGAc                        >  2:27149/1‑150 (MQ=255)
                                                                                                                              gCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGAcg                       >  1:1194747/1‑150 (MQ=255)
                                                                                                                               cGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGAcgc                      >  2:1180618/1‑150 (MQ=255)
                                                                                                                               cGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGAcgc                      >  2:566504/1‑150 (MQ=255)
                                                                                                                                  ttattctCGCGGCTATTGTCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCggg                   >  2:1067170/7‑150 (MQ=255)
                                                                                                                                  cGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCggg                   <  2:36439/150‑1 (MQ=255)
                                                                                                                                  cGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCggg                   <  2:42525/150‑1 (MQ=255)
                                                                                                                                   tattctCGCGGCTATTGTCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGt                  >  1:432519/6‑150 (MQ=255)
                                                                                                                                    gCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGtt                 <  1:1570900/150‑1 (MQ=255)
                                                                                                                                      tttCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTgg               <  2:1187543/150‑1 (MQ=255)
                                                                                                                                       ttCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGt              <  2:460355/150‑1 (MQ=255)
                                                                                                                                        tCGCGGCTATTGTCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTg             >  2:895846/1‑150 (MQ=255)
                                                                                                                                        tCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTg             <  2:575090/150‑1 (MQ=255)
                                                                                                                                        tCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTg             <  1:85495/150‑1 (MQ=255)
                                                                                                                                        tCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTg             <  1:27149/150‑1 (MQ=255)
                                                                                                                                        tCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTg             <  1:1174752/150‑1 (MQ=255)
                                                                                                                                         cgcgACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTGc            >  2:1524192/1‑150 (MQ=255)
                                                                                                                                         cgcgACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTGc            >  1:541235/1‑150 (MQ=255)
                                                                                                                                            gACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTGCATc         <  1:471250/150‑1 (MQ=255)
                                                                                                                                             aCAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTGCATCg        <  1:1056153/150‑1 (MQ=255)
                                                                                                                                               aaTAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTGCATCggg      >  2:1294524/1‑150 (MQ=255)
                                                                                                                                                  aGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTGCATCGGGGcc   <  2:518331/150‑1 (MQ=255)
                                                                                                                                                   gCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTGCATCGGGGCCg  <  1:858822/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCCATAAACCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTGCATCGGGGCCG  >  NC_000913/886048‑886344

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: