Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,329,209 A→G W197R (TGG→CGG)  yfaQ ← tandem DUF2300 domain protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,329,2090AG97.8% 1522.6 / 8.1 416W197R (TGG→CGG) yfaQtandem DUF2300 domain protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base A (3/6);  new base G (178/227);  total (183/233)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 419 total aligned reads displayed.

GGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGT  >  NC_000913/2329061‑2329356
                                                                                                                                                    |                                                                                                                                                   
ggCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCggg                                                                                                                                                   <  5:411511/151‑1 (MQ=255)
ggCAAAACCCGATTAAGCACCGGCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCggg                                                                                                                                                   <  1:1099993/151‑1 (MQ=255)
       cccGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCAc                                                                                                                                            <  6:1448159/151‑1 (MQ=255)
        ccGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCAcc                                                                                                                                           >  2:596504/1‑151 (MQ=255)
         cGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCAccc                                                                                                                                          >  5:83913/1‑151 (MQ=255)
           aTTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGt                                                                                                                                        <  8:311888/151‑1 (MQ=255)
             tAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTa                                                                                                                                       <  5:1041658/150‑1 (MQ=255)
             tAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTa                                                                                                                                       <  1:1489058/150‑1 (MQ=255)
              aaGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTAcc                                                                                                                                     <  2:166078/151‑1 (MQ=255)
                  aCCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGAc                                                                                                                                 <  1:188338/151‑1 (MQ=255)
                   ccGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGGGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACt                                                                                                                                <  2:1342329/151‑1 (MQ=255)
                    cGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACt                                                                                                                                >  5:61096/1‑150 (MQ=255)
                    cGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTa                                                                                                                               <  6:695287/151‑1 (MQ=255)
                     gTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTAc                                                                                                                              >  1:793937/1‑151 (MQ=255)
                      tCGCATAACGGCTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGCTTTTCACCCGTACCCGACTAcc                                                                                                                             >  1:293789/1‑151 (MQ=255)
                      gcGCATAACGGGTTAGTACCGCTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTAcc                                                                                                                             <  8:1288997/150‑1 (MQ=255)
                       cGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCTGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCg                                                                                                                            <  6:554495/151‑1 (MQ=255)
                       cGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCg                                                                                                                            <  3:308669/151‑1 (MQ=255)
                          aTAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGc                                                                                                                         >  4:811631/1‑151 (MQ=255)
                            aaCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCAt                                                                                                                       >  6:261123/1‑151 (MQ=255)
                             aCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATg                                                                                                                      >  8:604740/1‑151 (MQ=255)
                                ggTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCCCCCGTACCCGACTACCGAGCATGGAt                                                                                                                   <  3:692630/151‑1 (MQ=255)
                                 gTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATg                                                                                                                  <  8:1427884/151‑1 (MQ=255)
                                   tAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGc                                                                                                                 <  3:614713/150‑1 (MQ=255)
                                    aGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCg                                                                                                                >  2:1131304/1‑150 (MQ=255)
                                    aGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCg                                                                                                                >  3:1153512/1‑150 (MQ=255)
                                     gTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGATTTCACCCGTACCCGACTACCGAGCATGGATGCGac                                                                                                              >  4:244304/1‑151 (MQ=255)
                                       aCCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGacac                                                                                                            >  8:199300/1‑151 (MQ=255)
                                        ccGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACAcc                                                                                                           >  6:1507111/1‑151 (MQ=255)
                                          gTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCg                                                                                                          >  7:1437162/1‑150 (MQ=255)
                                          gTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGa                                                                                                         <  1:1200709/151‑1 (MQ=255)
                                               gCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGTCTACCGAGCATGGATGCGACACCGAtttt                                                                                                     <  1:1270533/150‑1 (MQ=255)
                                               gCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTc                                                                                                    >  5:1193160/1‑151 (MQ=255)
                                               gCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTc                                                                                                    >  5:1494041/1‑151 (MQ=255)
                                                  ggTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGc                                                                                                  <  4:1443088/150‑1 (MQ=255)
                                                   gTGCCACTGCCAGTAACCGACAGCGGTGTACCTTCGCTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCt                                                                                                <  4:1071063/151‑1 (MQ=255)
                                                   gTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGGTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCt                                                                                                <  8:1519170/151‑1 (MQ=255)
                                                     gacACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGATTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTc                                                                                               <  1:1628254/148‑1 (MQ=255)
                                                     gCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTCACCCCTACCCGACGACCGACCATGGGAGCGACACCGATTTTACCCTCGt                                                                                             >  2:980871/1‑151 (MQ=255)
                                                           gCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGc                                                                                        >  3:1663854/1‑151 (MQ=255)
                                                             cAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGcac                                                                                      <  6:971178/151‑1 (MQ=255)
                                                                tAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGcaccac                                                                                   >  5:1301345/1‑151 (MQ=255)
                                                                 aaCCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGcaccac                                                                                   >  1:535981/1‑150 (MQ=255)
                                                                    ccACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGGCACCGATTTTCGCCTCGTCCAGCACCACATcc                                                                               <  7:923819/151‑1 (MQ=255)
                                                                    ccACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATcc                                                                               <  5:952122/151‑1 (MQ=255)
                                                                       caGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGt                                                                            <  8:1050832/151‑1 (MQ=255)
                                                                       caGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGt                                                                            >  2:490820/1‑151 (MQ=255)
                                                                         gCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTaa                                                                          <  7:1053987/151‑1 (MQ=255)
                                                                         gCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTa                                                                           >  2:1427685/1‑150 (MQ=255)
                                                                                aCCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTcc                                                                   <  8:13141/151‑1 (MQ=255)
                                                                                  cATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTg                                                                 <  8:510666/151‑1 (MQ=255)
                                                                                  cATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTg                                                                 <  2:352684/151‑1 (MQ=255)
                                                                                  cATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTg                                                                 >  7:446900/1‑151 (MQ=255)
                                                                                   aTCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCCCATCCAGTAAGACTTCCTg                                                                 >  2:1256752/1‑150 (MQ=255)
                                                                                   aTCGGGTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGa                                                                <  7:577928/151‑1 (MQ=255)
                                                                                     cGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGc                                                              <  4:67009/151‑1 (MQ=255)
                                                                                     cGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGc                                                              <  7:73213/151‑1 (MQ=255)
                                                                                      ggTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCt                                                             >  5:326132/1‑151 (MQ=255)
                                                                                      ggTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCt                                                             >  4:1368871/1‑151 (MQ=255)
                                                                                       gTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCtt                                                            <  6:801462/151‑1 (MQ=255)
                                                                                          aaCCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGt                                                         >  4:396950/1‑151 (MQ=255)
                                                                                            ccagccagCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGttt                                                       >  3:798658/1‑151 (MQ=255)
                                                                                             cagccagcAAAACCACCCTGACTGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGtttt                                                      <  6:1250765/151‑1 (MQ=255)
                                                                                             cagccagcAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGtttt                                                      <  4:483587/151‑1 (MQ=255)
                                                                                                ccagcAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTgcg                                                   >  5:693536/1‑151 (MQ=255)
                                                                                                   gcAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCgg                                                <  2:767988/151‑1 (MQ=255)
                                                                                                   gcAAAACCACCCTCACGGATTTCCTGTTTATCATCGGCAACCCATCTCCGGCTTATCACCCGTCCCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCgg                                                <  1:689232/151‑1 (MQ=255)
                                                                                                    cAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCggg                                               >  7:856816/1‑151 (MQ=255)
                                                                                                      aaaCCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTACACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTa                                             <  2:1386870/151‑1 (MQ=255)
                                                                                                      aaaCCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGGTTTCACCCGTACCCGACAACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTa                                             <  3:1273879/151‑1 (MQ=255)
                                                                                                         ccaccCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACg                                          >  5:294301/1‑151 (MQ=255)
                                                                                                            cccTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCg                                       <  6:1332594/151‑1 (MQ=255)
                                                                                                             ccTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGc                                      <  3:894256/151‑1 (MQ=255)
                                                                                                                gACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTa                                    >  3:1412494/1‑150 (MQ=255)
                                                                                                                gACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTa                                    >  6:618846/1‑150 (MQ=255)
                                                                                                                 ccGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAAt                                  <  8:1234909/150‑1 (MQ=255)
                                                                                                                 aCGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAAt                                  >  2:1018910/1‑151 (MQ=255)
                                                                                                                 aCGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAAt                                  <  7:1509560/151‑1 (MQ=255)
                                                                                                                 aCGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACACGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCATGAGCTTTGTTTTGCGCGGGTAACGTCGCTAAt                                  >  5:1225283/1‑151 (MQ=255)
                                                                                                                  cGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCCACCATGGATGCGACACCGATAATCGCCTCGTCCACCACCACATCCAGTAAGACTTCATGAGCTTTGTATTGAGCGGGTAACGTCGCTAAtt                                 >  5:1615121/1‑150 (MQ=255)
                                                                                                                   ggATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATg                                 >  5:1649004/1‑150 (MQ=255)
                                                                                                                   ggATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATg                                 >  1:880401/1‑150 (MQ=255)
                                                                                                                   ggATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGa                                <  8:920378/151‑1 (MQ=255)
                                                                                                                    gATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGaa                               <  1:1370493/151‑1 (MQ=255)
                                                                                                                    gATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGa                                >  4:1320609/1‑150 (MQ=255)
                                                                                                                      tttCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATc                             >  4:1618965/1‑151 (MQ=255)
                                                                                                                       ttCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATc                             <  6:453029/150‑1 (MQ=255)
                                                                                                                        tCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCa                            <  3:478374/150‑1 (MQ=255)
                                                                                                                        tCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCa                            >  5:2073/1‑150 (MQ=255)
                                                                                                                         ccTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACGCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTATTGAATCAag                          <  8:772675/151‑1 (MQ=255)
                                                                                                                             tttATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAagcagc                      >  7:362487/1‑151 (MQ=255)
                                                                                                                                   atcGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTa                <  2:631563/151‑1 (MQ=255)
                                                                                                                                      ggCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAAcc             <  3:505986/151‑1 (MQ=255)
                                                                                                                                        cAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGt           <  6:378227/151‑1 (MQ=255)
                                                                                                                                        cAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGt           >  2:905079/1‑151 (MQ=255)
                                                                                                                                        cAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGt           <  7:256009/151‑1 (MQ=255)
                                                                                                                                          aaCCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGtt          >  5:54470/1‑150 (MQ=255)
                                                                                                                                            ccAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTAcc       <  4:289718/151‑1 (MQ=255)
                                                                                                                                               gCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAg    <  1:1014862/151‑1 (MQ=255)
                                                                                                                                                 tCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGt  >  3:1533004/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGT  >  NC_000913/2329061‑2329356

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