Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,785 T→A 8.0% T167S (ACA→TCA)  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,7850TA 8.0% 289.5 / 10.9 97T167S (ACA→TCA) 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.73e-01

CGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAA  >  NC_000913/1446636‑1446934
                                                                                                                                                     |                                                                                                                                                     
cGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGt                                                                                                                                                       <  2:416320/150‑1 (MQ=255)
    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:1398006/150‑1 (MQ=255)
                                            ccACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGAt                                                                                                           <  1:651610/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:1354805/150‑1 (MQ=255)
                                                                 cACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTa                                                                                      >  2:1088392/1‑150 (MQ=255)
                                                                 cACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTa                                                                                      <  1:960067/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:127843/1‑150 (MQ=255)
                                                                                                                      acCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGt                                 <  1:538332/150‑1 (MQ=255)
                                                                                                                      acCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGt                                 <  2:1010845/150‑1 (MQ=255)
                                                                                                                       ccATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGtt                                <  1:1283565/150‑1 (MQ=255)
                                                                                                                       ccATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGtt                                >  2:1185123/1‑150 (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               >  2:54561/5‑150 (MQ=255)
                                                                                                                                        ccggcACTTTCTGAAACAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGc               >  1:204819/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)
                                                                                                                                                     |                                                                                                                                                     
CGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAA  >  NC_000913/1446636‑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: