Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,022,007 C→T 9.4% G21G (GGC→GGT tatA → TatABCE protein translocation system subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,022,0070CT 9.4% 340.6 / 19.1 115G21G (GGC→GGTtatATatABCE protein translocation system subunit
Reads supporting (aligned to +/- strand):  ref base C (51/53);  new base T (5/6);  total (56/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.33e-01

Only 100 of 116 total aligned reads displayed.

ACGCTGATTTTTTCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAGGCG  >  NC_000913/4021860‑4022151
                                                                                                                                                   |                                                                                                                                                
acGCTGATTTTTTCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTCTTTGGc                                                                                                                                                  <  2:486192/148‑1 (MQ=255)
 cGCTGATTTTTTCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCa                                                                                                                                                 >  2:429966/1‑148 (MQ=255)
       ttttttCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCaaaa                                                                                                                                           >  1:528169/1‑148 (MQ=255)
          tttCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                                                                       <  3:84392/149‑2 (MQ=255)
          tttCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGc                                                                                                                                        <  1:484531/148‑1 (MQ=255)
            tCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTc                                                                                                                                      <  1:697557/148‑1 (MQ=255)
            tCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCg                                                                                                                                     <  3:806572/149‑1 (MQ=255)
              aTCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGc                                                                                                                                   >  3:233721/1‑149 (MQ=255)
               tCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGc                                                                                                                                   >  1:101353/1‑148 (MQ=255)
               tCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCt                                                                                                                                  >  3:386241/1‑149 (MQ=255)
                   tCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCAc                                                                                                                                                <  2:492862/131‑1 (MQ=255)
                   tCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCAc                                                                                                                                                >  1:492862/1‑131 (MQ=255)
                   tCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCAt                                                                                                                              >  3:246194/1‑149 (MQ=255)
                    cAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGc                                                                                                                                        >  2:581522/1‑138 (MQ=255)
                    cAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGc                                                                                                                                        <  1:581522/138‑1 (MQ=255)
                    cAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCAt                                                                                                                              >  1:130819/1‑148 (MQ=255)
                     aaGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATc                                                                                                                             <  1:429966/148‑1 (MQ=255)
                      aGGCGGGCCGTGTAACGTATAATGCGGCTTTGCTTAATCATCATCTACCACAGAGGAACATCTATGGGTGGTATCAGTATTTGGCTGTTAATGATTATTGCCGTCATCGTTGCACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCgg                                                                                                                           <  3:562804/149‑1 (MQ=255)
                       ggCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCgg                                                                                                                           >  2:91149/1‑148 (MQ=255)
                          gggCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACTACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                                                                       <  3:574257/149‑1 (MQ=255)
                              cGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGAt                                                                                                                    >  2:242039/1‑148 (MQ=255)
                                  tAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTg                                                                                                                <  2:883577/148‑1 (MQ=255)
                                    aCGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTg                                                                                                             >  3:425668/1‑149 (MQ=255)
                                        ataATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGt                                                                                                          >  1:702012/1‑148 (MQ=255)
                                           aTGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGAt                                                                                                      <  4:582639/149‑1 (MQ=255)
                                             gCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATc                                                                                                     >  2:871308/1‑148 (MQ=255)
                                               ggCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaa                                                                                                   <  1:91149/148‑1 (MQ=255)
                                               ggCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTACTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaa                                                                                                   <  2:689286/148‑1 (MQ=255)
                                                gCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGACAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaaa                                                                                                  >  2:209961/1‑148 (MQ=255)
                                                  tttGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGt                                                                                                          >  2:804483/1‑138 (MQ=255)
                                                  tttGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGt                                                                                                          <  1:804483/138‑1 (MQ=255)
                                                  tttGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGc                                                                                               >  3:643841/1‑149 (MQ=255)
                                                      tttAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCTCCAAAAAGGTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAg                                                                                                 >  3:780317/1‑143 (MQ=255)
                                                      tttAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaaac                                                                                                 <  4:780317/143‑2 (MQ=255)
                                                      tttAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCttt                                                                                            >  1:878013/1‑148 (MQ=255)
                                                       ttAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATc                                                                                                                   <  2:713204/124‑1 (MQ=255)
                                                       ttAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATc                                                                                                                   >  1:713204/1‑124 (MQ=255)
                                                           tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaaaaa                                                                                      <  4:163489/149‑1 (MQ=255)
                                                             atcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCtt                                                                                                                 <  2:826714/120‑1 (MQ=255)
                                                             atcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCtt                                                                                                                 >  1:826714/1‑120 (MQ=255)
                                                              tcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCa                                                                                   >  3:281850/1‑149 (MQ=255)
                                                                  cTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGa                                                                               >  4:148422/1‑149 (MQ=255)
                                                                     ccACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCg                                                                            >  4:13047/1‑149 (MQ=255)
                                                                      cacaGAGGAACGTGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCg                                                                            <  1:355084/148‑1 (MQ=255)
                                                                      cacaGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCga                                                                           >  4:690168/1‑149 (MQ=255)
                                                                       acaGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgat                                                                          <  4:850450/149‑1 (MQ=255)
                                                                       acaGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCga                                                                           >  2:471267/1‑148 (MQ=255)
                                                                          gagGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgatgat                                                                       >  4:315162/1‑149 (MQ=255)
                                                                             gAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGAtgatga                                                                     <  1:471267/148‑1 (MQ=255)
                                                                             gAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGaa                                                                    <  4:246194/149‑1 (MQ=255)
                                                                              aaCATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGaa                                                                    >  2:578603/1‑148 (MQ=255)
                                                                              aaCATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAAc                                                                   <  3:518593/149‑1 (MQ=255)
                                                                                     aTGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGc                                                             <  2:563007/148‑1 (MQ=255)
                                                                                         gtggtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGAt                                                        >  4:210849/1‑149 (MQ=255)
                                                                                          tggtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGAt                                                        <  2:130819/148‑1 (MQ=255)
                                                                                          tggtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATa                                                       >  4:705759/1‑149 (MQ=255)
                                                                                            gtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaaa                                                                                                  >  3:824355/1‑104 (MQ=255)
                                                                                            gtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaaa                                                                                                  <  4:824355/104‑1 (MQ=255)
                                                                                             tATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATaaa                                                     <  1:242039/148‑1 (MQ=255)
                                                                                                cAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAAcc                                                  >  2:30408/1‑148 (MQ=255)
                                                                                                  gTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGt                                               >  4:294947/1‑149 (MQ=255)
                                                                                                  gTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGt                                               >  3:487368/1‑149 (MQ=255)
                                                                                                   tATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaaaaa                                                                                      <  4:894929/109‑1 (MQ=255)
                                                                                                   tATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaaaaa                                                                                      >  3:894929/1‑109 (MQ=255)
                                                                                                    aTTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTc                                              >  2:762373/1‑148 (MQ=255)
                                                                                                    aTTTGGCAGTTATTGATATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGc                                                                                                                                        >  1:802109/1‑57 (MQ=25)
                                                                                                          cAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGc                                       <  4:281850/149‑1 (MQ=255)
                                                                                                           aGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCt                                      <  4:238057/149‑1 (MQ=255)
                                                                                                              tATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCTATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGa                                    <  2:468658/148‑1 (MQ=255)
                                                                                                                 tGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGc                                                                                                                                        <  1:751312/45‑1 (MQ=37)
                                                                                                                 tGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGc                                                                                                                                        >  2:751312/1‑45 (MQ=37)
                                                                                                                 tGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGAtttt                                >  4:200963/1‑149 (MQ=255)
                                                                                                                  gATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTa                               <  4:639/149‑1 (MQ=255)
                                                                                                                   attattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTa                               <  1:397955/148‑1 (MQ=255)
                                                                                                                    ttattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCa                                                 <  2:768082/129‑1 (MQ=255)
                                                                                                                    ttattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCa                                                 >  1:768082/1‑129 (MQ=255)
                                                                                                                    ttattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTAc                              >  2:582547/1‑148 (MQ=255)
                                                                                                                    ttattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTAc                              >  1:531124/1‑148 (MQ=255)
                                                                                                                    ttattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACt                             >  4:528195/1‑149 (MQ=255)
                                                                                                                     tattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCt                                                                                                                  >  1:759002/1‑63 (MQ=255)
                                                                                                                     tattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCt                                                                                                                  <  2:759002/63‑1 (MQ=255)
                                                                                                                          ccGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGaa                        >  1:660779/1‑148 (MQ=255)
                                                                                                                           cGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGaaaa                      <  3:96470/149‑1 (MQ=255)
                                                                                                                               aTCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTa                   <  2:528169/148‑1 (MQ=255)
                                                                                                                               aTCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTAt                  <  4:487368/149‑1 (MQ=255)
                                                                                                                               aTCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTAt                  <  3:528195/149‑1 (MQ=255)
                                                                                                                                tCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTAt                  <  1:446720/148‑1 (MQ=255)
                                                                                                                                tCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTAt                  >  2:602153/1‑148 (MQ=255)
                                                                                                                                tCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATc                 <  4:37866/149‑1 (MQ=255)
                                                                                                                                  gTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGc               >  3:185044/1‑149 (MQ=255)
                                                                                                                                     gTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGa            <  4:643841/149‑1 (MQ=255)
                                                                                                                                     gTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGa            <  3:315162/149‑1 (MQ=255)
                                                                                                                                       cgagCTTATTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAACCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATa          >  4:567658/4‑149 (MQ=255)
                                                                                                                                          gCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGaaa                       >  1:488921/1‑133 (MQ=255)
                                                                                                                                          gCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGaaa                       <  2:488921/133‑1 (MQ=255)
                                                                                                                                          gCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAg        >  2:54333/1‑148 (MQ=255)
                                                                                                                                           cTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGa                         <  3:41116/130‑1 (MQ=255)
                                                                                                                                            tttttGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAg     >  3:13168/1‑149 (MQ=255)
                                                                                                                                             ttttGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAg     <  2:878013/148‑1 (MQ=255)
                                                                                                                                               ttGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAGGc   <  1:582547/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ACGCTGATTTTTTCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAGGCG  >  NC_000913/4021860‑4022151

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: