Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,786 T→A 8.2% E166D (GAA→GAT) ‡ feaR ← DNA‑binding transcriptional activator FeaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,446,7860TA 8.2% 292.4 / 12.6 97E166D (GAA→GAT) ‡feaRDNA‑binding transcriptional activator FeaR
Reads supporting (aligned to +/- strand):  ref base T (47/42);  new base A (5/3);  total (52/45)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

CGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAA  >  NC_000913/1446640‑1446934
                                                                                                                                                  |                                                                                                                                                    
cGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCa                                                                                                                                                   >  2:1448182/1‑150 (MQ=255)
    atatCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGaaa                                                                                                                                               >  1:6658/1‑150 (MQ=255)
     tatCTCTTCGCGAATATTATCGTCTATAAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAg                                                                                                                                              >  2:72877/1‑150 (MQ=255)
          ctTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACGCACCATCGCCTGTAGCGCAGCTTCACTTTCTGAAACAGAAAGTgaag                                                                                                                                         <  2:894574/150‑4 (MQ=255)
          ctTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCg                                                                                                                                         >  1:559620/1‑150 (MQ=255)
           ttCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGAAACAGAAAGTgaagc                                                                                                                                        <  2:235795/150‑5 (MQ=255)
            tCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGAAACAGAAAGTgaagct                                                                                                                                       <  2:246689/150‑6 (MQ=255)
             cgcgAATATTATCGTCTATCAACGTACCCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCTCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGat                                                                                                                                      >  1:817080/1‑150 (MQ=255)
             cgcgAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGat                                                                                                                                      >  1:583862/1‑150 (MQ=255)
             cgcgAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTCCCGGat                                                                                                                                      >  1:350853/1‑150 (MQ=255)
                gAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGattat                                                                                                                                   >  2:289271/1‑150 (MQ=255)
                 aaTATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGattatt                                                                                                                                  <  2:948724/150‑1 (MQ=255)
                  atatTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTc                                                                                                                                 <  1:318110/150‑1 (MQ=255)
                   tattatCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCa                                                                                                                                <  1:289271/150‑1 (MQ=255)
                      tatCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGc                                                                                                                             <  1:1167324/150‑1 (MQ=255)
                       atCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGct                                                                                                                            >  1:1010845/1‑150 (MQ=255)
                           tCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTcc                                                                                                                        >  2:1283565/1‑150 (MQ=255)
                           tCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTcc                                                                                                                        <  1:33905/150‑1 (MQ=255)
                            cTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCt                                                                                                                       >  2:753710/1‑150 (MQ=255)
                             tATCAACTTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCTGATTATTCATGCTCTCCTg                                                                                                                      <  1:538587/150‑1 (MQ=255)
                             tATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTg                                                                                                                      <  1:1365151/150‑1 (MQ=255)
                             tATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTg                                                                                                                      <  2:662988/150‑1 (MQ=255)
                                cAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTaa                                                                                                                   <  1:1216027/150‑1 (MQ=255)
                                   cGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAg                                                                                                                >  1:1057470/1‑150 (MQ=255)
                                    gTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAgg                                                                                                               <  1:114933/150‑1 (MQ=255)
                                      aaCCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTGACAGGCg                                                                                                             <  1:65461/150‑1 (MQ=255)
                                        ccACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGAt                                                                                                           >  2:679526/1‑150 (MQ=255)
                                        ccACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGAt                                                                                                           <  1:651610/150‑1 (MQ=255)
                                        ccACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGAt                                                                                                           <  1:1398006/150‑1 (MQ=255)
                                         cACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACAGCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATg                                                                                                          <  1:1522015/150‑1 (MQ=255)
                                          aCTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGa                                                                                                         <  1:1614506/150‑1 (MQ=255)
                                          aCTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGa                                                                                                         >  2:357764/1‑150 (MQ=255)
                                            tttttGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGat                                                                                                                                      <  2:745034/119‑1 (MQ=255)
                                            tttttGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGat                                                                                                                                      >  1:745034/1‑119 (MQ=255)
                                             ttttGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTg                                                                                                      >  1:846224/1‑150 (MQ=255)
                                                  aaaCTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTg                                                                                                 >  2:344506/1‑150 (MQ=255)
                                                   aaCTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAATGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGc                                                                                                <  1:1316311/150‑1 (MQ=255)
                                                      tGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCAcc                                                                                             >  1:912125/1‑150 (MQ=255)
                                                            tCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGt                                                                                       <  1:1344574/150‑1 (MQ=255)
                                                             cACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTa                                                                                      <  1:864422/150‑1 (MQ=255)
                                                             cACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTa                                                                                      <  1:960067/150‑1 (MQ=255)
                                                             cACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTa                                                                                      >  2:1088392/1‑150 (MQ=255)
                                                             cACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTa                                                                                      <  1:1354805/150‑1 (MQ=255)
                                                              acgacgAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTaa                                                                                     >  1:1555452/1‑150 (MQ=255)
                                                                   gaGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAg                                                                                <  2:527758/150‑1 (MQ=255)
                                                                    aGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTATGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGc                                                                               >  2:1233161/1‑150 (MQ=255)
                                                                     ggTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCg                                                                              >  1:161393/1‑150 (MQ=255)
                                                                        tGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTcc                                                                           >  1:751639/1‑150 (MQ=255)
                                                                           aCAGATTCCCGCTGATGAAGTACCGGGGGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGt                                                                        >  1:445681/1‑150 (MQ=255)
                                                                             aGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCt                                                                      >  1:222686/1‑150 (MQ=255)
                                                                                   ccGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTgc                                                                >  2:1380281/1‑150 (MQ=255)
                                                                                    cGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTgcg                                                               >  1:578252/1‑150 (MQ=255)
                                                                                    cGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTgcg                                                               <  2:1172428/150‑1 (MQ=255)
                                                                                     gCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTgcgc                                                              >  2:1394414/1‑150 (MQ=255)
                                                                                      cTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCa                                                             <  1:406372/150‑1 (MQ=255)
                                                                                      cTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCa                                                             <  1:1448182/150‑1 (MQ=255)
                                                                                           gaAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAg                                                        >  2:188083/1‑150 (MQ=255)
                                                                                                aCCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGtttt                                                   >  2:438013/1‑150 (MQ=255)
                                                                                                 ccGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCTTCCAGTCTTTCTGCGCAGATAGGttttt                                                  <  1:570819/150‑1 (MQ=255)
                                                                                                 ccGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGttttt                                                  >  2:538332/1‑150 (MQ=255)
                                                                                                  cGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTg                                                 >  1:488476/1‑150 (MQ=255)
                                                                                                     gcgcAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATg                                              <  2:578252/150‑1 (MQ=255)
                                                                                                      cgcAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATgg                                             >  1:599719/1‑150 (MQ=255)
                                                                                                           cagACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGaaa                                        >  2:678663/1‑150 (MQ=255)
                                                                                                               cacacCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAatat                                    <  2:559620/150‑1 (MQ=255)
                                                                                                                  acCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGt                                 <  1:538332/150‑1 (MQ=255)
                                                                                                                  acCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGt                                 >  1:127843/1‑150 (MQ=255)
                                                                                                                  acCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGt                                 <  2:1010845/150‑1 (MQ=255)
                                                                                                                   ccATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGtt                                >  2:1185123/1‑150 (MQ=255)
                                                                                                                   ccATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGtt                                <  1:1283565/150‑1 (MQ=255)
                                                                                                                    cATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTc                               >  1:751324/1‑150 (MQ=255)
                                                                                                                       cGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCcag                            <  2:912125/150‑1 (MQ=255)
                                                                                                                       cGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCcag                            >  2:923618/1‑150 (MQ=255)
                                                                                                                        gCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCcagc                           <  2:817080/150‑1 (MQ=255)
                                                                                                                        gCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCcagc                           <  1:357764/150‑1 (MQ=255)
                                                                                                                           tGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCaga                        >  1:1323032/1‑150 (MQ=255)
                                                                                                                            gTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGACTTTCTGCGCAGATAGGTTTTTGATCGGGAAAATATTGTTCCAGCagag                       >  2:1254276/1‑150 (MQ=255)
                                                                                                                              aGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTg                     >  1:831160/1‑150 (MQ=255)
                                                                                                                                cgcAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCg                   <  2:599719/150‑1 (MQ=255)
                                                                                                                                  cAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTg                 <  2:751639/150‑1 (MQ=255)
                                                                                                                                  cAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTg                 <  1:996414/150‑1 (MQ=255)
                                                                                                                                   aGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTgg                <  1:753710/150‑1 (MQ=255)
                                                                                                                                    gCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGc               <  2:161393/150‑1 (MQ=255)
                                                                                                                                    ccggcACTTTCTGAAACAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGc               >  1:204819/5‑150 (MQ=255)
                                                                                                                                    ccggcACTTTCTGAAACAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGc               >  2:54561/5‑150 (MQ=255)
                                                                                                                                     cggcACTTTCTGAAACAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCa              >  2:299025/4‑150 (MQ=255)
                                                                                                                                      ggcACTTTCTGAAACAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGTTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCaa             >  1:712264/3‑150 (MQ=255)
                                                                                                                                      ggcACTTTCTGAAACAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCaa             >  1:222309/3‑150 (MQ=255)
                                                                                                                                       tCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCaaa            >  1:1242912/1‑150 (MQ=255)
                                                                                                                                       tCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACAATGGGTAAGTCAGCGTCCAGTCTTTCTGCTCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCaaa            <  2:1442834/150‑1 (MQ=255)
                                                                                                                                         aCTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAg          >  1:889937/1‑150 (MQ=255)
                                                                                                                                         aCTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGATGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAg          <  2:846224/150‑1 (MQ=255)
                                                                                                                                          cTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGt         >  2:1158318/1‑150 (MQ=255)
                                                                                                                                          cTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGt         >  2:111132/1‑150 (MQ=255)
                                                                                                                                           tttCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTa        >  2:20815/1‑150 (MQ=255)
                                                                                                                                                gTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGaaaa                                       <  2:1098013/114‑1 (MQ=255)
                                                                                                                                                gTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGaaaa                                       >  1:1098013/1‑114 (MQ=255)
                                                                                                                                                 tttCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGgaa  <  2:6658/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAA  >  NC_000913/1446640‑1446934

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: