Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,569,172 G→A pseudogene (175/921 nt) yjiP → pseudogene, transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,569,1720GA99.0% 350.1 / NA 96pseudogene (175/921 nt)yjiPpseudogene, transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (42/53);  minor base C (0/1);  total (42/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GCATGATGCATTATTTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGGCGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCATAT  >  NC_000913/4569024‑4569317
                                                                                                                                                    |                                                                                                                                                 
gccttctgCATTATTTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGacg                                                                                                                                                 >  1:1013855/7‑151 (MQ=255)
 cATGATGCATTATTTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGCTTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCAAGCTTTGTCGATGAAAAATTGCGGacg                                                                                                                                                 >  2:952010/1‑150 (MQ=255)
  atgatgCATTATTTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCt                                                                                                                                               >  1:1080703/1‑151 (MQ=255)
  atgatgCATTATTTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCt                                                                                                                                               <  1:1169315/151‑1 (MQ=255)
     atgCATTATTTAAAACCTTTCTCACGCTCCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGTTaca                                                                                                                                            <  2:847275/151‑1 (MQ=255)
       gCATTATTTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACt                                                                                                                                          >  1:1045216/1‑151 (MQ=255)
         attattTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTcc                                                                                                                                        <  1:116/151‑1 (MQ=255)
         attattTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATACACTTACCTAAAGAGTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTcc                                                                                                                                        <  2:363013/151‑1 (MQ=255)
           tattTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCg                                                                                                                                       <  2:777705/150‑1 (MQ=255)
            attTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGat                                                                                                                                     <  1:713126/151‑1 (MQ=255)
              ttAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCTATATCTTAAAACTGTAATCCGCCAGCTTTGTCGATGAAAAAGTGCGGACGCTACACTCCGata                                                                                                                                    <  2:321188/150‑1 (MQ=255)
              ttAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGatat                                                                                                                                   >  1:261748/1‑151 (MQ=255)
               tAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTCCCTAACGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATCtt                                                                                                                                  >  2:409964/1‑151 (MQ=255)
               tAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTACGATAtt                                                                                                                                  >  2:146932/1‑151 (MQ=255)
                aaaaCCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATAtt                                                                                                                                  <  1:211102/150‑1 (MQ=255)
                  aaCCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATAtttt                                                                                                                                >  2:631160/1‑150 (MQ=255)
                   aCCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTAt                                                                                                                              <  1:340385/151‑1 (MQ=255)
                   aCCTTTCCCACGCACCCTGTCACTGCACGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTATGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCGACACTCCGTTCTTTTAt                                                                                                                              >  2:191714/1‑151 (MQ=255)
                     cTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAATACTGGAATCCGCCAGCGGTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATg                                                                                                                             <  2:731773/150‑1 (MQ=255)
                       ttCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTc                                                                                                                          >  1:642088/1‑151 (MQ=255)
                        tctcACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCg                                                                                                                         >  1:999426/1‑151 (MQ=255)
                           cACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGt                                                                                                                       >  1:650146/1‑150 (MQ=255)
                              gCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTaa                                                                                                                     >  1:43073/1‑149 (MQ=255)
                                 cccTGACCCTGCGCTGGATTTTATGGAGAGTCACTTACCTAAAGATTTACGTTAACTTTGCGATCTCGTTCGCTTAAAACTGGAGTCCGCCATCATTGTCGATGAAAAAGTGCGGACGCTCCTCTCCGATATTTTATGGTCAGTTAAGAcc                                                                                                                <  2:84003/151‑1 (MQ=255)
                                           gcgcGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCg                                                                                                      <  1:463407/151‑1 (MQ=255)
                                            cgcgGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCg                                                                                                      >  1:880153/1‑150 (MQ=255)
                                             gcgGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGAt                                                                                                    <  2:143990/151‑1 (MQ=255)
                                              cgGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATg                                                                                                   <  1:482554/151‑1 (MQ=255)
                                                 gATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCt                                                                                                <  2:1013783/151‑1 (MQ=255)
                                                    tttATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCtat                                                                                              <  1:316326/150‑1 (MQ=255)
                                                     ttATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTatat                                                                                            >  1:83484/1‑151 (MQ=255)
                                                      tATGGAGATTCACTTGCCTAAAGATTTACGTGAACTGTGCGTTCTCGATAGCTTAAACCTGGTTTACGTCAGCGTTGTCCATGAAAAATTGCGGCCTCTACACTCCGATATTTTATGGTCTCTAAAGACCCGCGAAGGCGATTGCTatat                                                                                            <  2:411423/150‑1 (MQ=255)
                                                      tATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTatat                                                                                            >  2:1136554/1‑150 (MQ=255)
                                                      tATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATc                                                                                           >  1:852881/1‑151 (MQ=255)
                                                        tGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTatat                                                                                            >  2:956783/1‑148 (MQ=255)
                                                         ggAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTAt                                                                                        >  1:585201/1‑151 (MQ=255)
                                                           agaTTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTAtgt                                                                                      <  2:1138430/151‑1 (MQ=255)
                                                               tCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATgtggt                                                                                   <  2:434719/150‑1 (MQ=255)
                                                               tCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATgtggt                                                                                   <  1:487988/150‑1 (MQ=255)
                                                                  cTTACCTTAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGAtt                                                                               >  2:65513/1‑151 (MQ=255)
                                                                  cTTACCGTAAGATTTACGTGAACTGTGCGACCTCGATAGCGTAAAACTGGAACCGGCCAGGTTTGTCGATGAAAAATTGCGGACGCTACACTCGGATAGTTTATGGTCAGTAAACTCGCGCGAAGGCGATGGCTATATCTAtgatgtgcc                                                                                <  2:128876/150‑8 (MQ=255)
                                                                     aCCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTg                                                                              >  2:1016788/1‑149 (MQ=255)
                                                                      ccTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGTTGGCTATATCTATGTGGTGATTGaa                                                                            >  2:647324/1‑150 (MQ=255)
                                                                      ccTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAAc                                                                           >  1:896886/1‑151 (MQ=255)
                                                                         aaaGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACa                                                                          <  2:732621/149‑1 (MQ=255)
                                                                         aaaGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACa                                                                          <  1:333408/149‑1 (MQ=255)
                                                                         aaaGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACAt                                                                         <  2:277042/150‑1 (MQ=255)
                                                                         aaaGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACAt                                                                         <  1:498622/150‑1 (MQ=255)
                                                                         aaaGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATc                                                                        >  2:66208/1‑151 (MQ=255)
                                                                         aaaGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATc                                                                        >  1:1150019/1‑151 (MQ=255)
                                                                          aaGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTTGAGGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCTGTAAAGACGCGCGAAGGCGATGGCTATATCTATGGGTTGATTGAACATCa                                                                       <  2:1123034/151‑1 (MQ=255)
                                                                          aaGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCa                                                                       <  1:271/151‑1 (MQ=255)
                                                                           aGATTTACGTGAACTGTGCGATCTCGAGAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTCCACTCCGATATTTTATGGTCGGTAAAGACCCTCGAAGGCGATGGCTATATCTATGTGGTGAGTGAACATCag                                                                      <  2:798621/151‑1 (MQ=255)
                                                                            gATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATc                                                                        >  1:159696/1‑148 (MQ=255)
                                                                            gATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATc                                                                        >  1:491504/1‑148 (MQ=255)
                                                                            gATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCaga                                                                     <  1:599746/151‑1 (MQ=255)
                                                                                 aCGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCcg                                                                 <  2:890973/150‑1 (MQ=255)
                                                                                   gTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCcgcg                                                               >  1:656436/1‑150 (MQ=255)
                                                                                   gTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGa                                                              <  1:929663/151‑1 (MQ=255)
                                                                                     gAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAgg                                                            <  1:706651/151‑1 (MQ=255)
                                                                                      aaCTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGa                                                           <  1:969187/151‑1 (MQ=255)
                                                                                          gtgCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACa                                                         <  1:361783/149‑1 (MQ=255)
                                                                                           tgCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACAtt                                                       >  2:1173621/1‑150 (MQ=255)
                                                                                           tgCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTc                                                      >  1:406325/1‑151 (MQ=255)
                                                                                            gCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTc                                                      <  1:295213/150‑1 (MQ=255)
                                                                                            gCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTc                                                      <  1:930817/150‑1 (MQ=255)
                                                                                              gacctcGTTAGCTTAAAACTGGAATCCGCCAGATTTGTCGATGAAAAATTGCGGACGCTACACTCCGATGTTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATGTATGTGGTGATTGAACAACAGAGCCGCGAGGACATTCaga                                                   <  2:791885/148‑3 (MQ=255)
                                                                                                   cGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCatat                                                  <  2:202009/147‑1 (MQ=255)
                                                                                                    gATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCt                                             <  1:1179690/151‑1 (MQ=255)
                                                                                                     aTAGCTTAAAACTGTAATCCGCCAGCTTGGTCGATGAAAAATTGCGGACGCTACACTCCTATATTTTATGGTCGATAAAGACCCGCGAAGGCGATGGCTATATCTATTTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCt                                             <  2:915327/150‑1 (MQ=255)
                                                                                                     aTAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCtt                                            >  2:761324/1‑151 (MQ=255)
                                                                                                     aTAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCtt                                            <  2:620660/151‑1 (MQ=255)
                                                                                                        gCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCg                                         <  2:276295/151‑1 (MQ=255)
                                                                                                         cTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGATACACTCAGACATTTTATGGTCGGTAAAGACCAGCGAAGGCGATGGCTATATCTCTGTGGCGATTGAACATCAGAGCTGCGAGGACATTCATATGTCCTTTc                                          >  2:1125/1‑149 (MQ=255)
                                                                                                             aaaCTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTg                                     >  2:47482/1‑150 (MQ=255)
                                                                                                              aaCTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGAt                                   <  1:943618/151‑1 (MQ=255)
                                                                                                                     aTCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTAt                            >  1:606235/1‑151 (MQ=255)
                                                                                                                        cGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTcc                         >  1:408511/1‑151 (MQ=255)
                                                                                                                         gCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTcc                         >  2:307062/1‑150 (MQ=255)
                                                                                                                         gCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTc                          >  2:721627/1‑149 (MQ=255)
                                                                                                                         gCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCa                        <  1:583321/151‑1 (MQ=255)
                                                                                                                          ccAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCAt                       <  1:1070895/151‑1 (MQ=255)
                                                                                                                           cAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATg                      <  1:1103050/151‑1 (MQ=255)
                                                                                                                             gCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGc                    >  1:300106/1‑151 (MQ=255)
                                                                                                                             gCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGc                    <  1:929932/151‑1 (MQ=255)
                                                                                                                                ttGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCg                   <  2:896315/149‑1 (MQ=255)
                                                                                                                                  gTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGa               <  1:682344/151‑1 (MQ=255)
                                                                                                                                    cGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATg             <  2:897388/151‑1 (MQ=255)
                                                                                                                                     gATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGAt              <  2:535014/149‑1 (MQ=255)
                                                                                                                                     gATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGt            <  2:43067/151‑1 (MQ=255)
                                                                                                                                      aTGAAAAATTGCGGACGCTTCACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGAt              <  2:204310/148‑1 (MQ=255)
                                                                                                                                       tGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGt            >  1:657737/1‑149 (MQ=255)
                                                                                                                                          aaaaTTGGGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCGATGGCGGGGATGTAgcg        >  2:45596/1‑150 (MQ=255)
                                                                                                                                            aaTTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTCGAGcc      >  2:404237/1‑150 (MQ=255)
                                                                                                                                            aaTTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCa     <  2:146074/151‑1 (MQ=255)
                                                                                                                                             aTTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCat    >  2:611508/1‑151 (MQ=255)
                                                                                                                                                gCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCatat  >  1:1134563/1‑150 (MQ=255)
                                                                                                                                                 cGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCatat  >  1:401990/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCATGATGCATTATTTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGGCGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCATAT  >  NC_000913/4569024‑4569317

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