Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,022,008 A→G 9.8% T22A (ACC→GCC)  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,0080AG 9.8% 335.1 / 20.5 113T22A (ACC→GCC) tatATatABCE protein translocation system subunit
Reads supporting (aligned to +/- strand):  ref base A (48/53);  new base G (5/6);  total (54/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 = 9.86e-01

Only 100 of 123 total aligned reads displayed.

CGCTGATTTTTTCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAGGCGGATAC  >  NC_000913/4021861‑4022156
                                                                                                                                                   |                                                                                                                                                    
cGCTGATTTTTTCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCa                                                                                                                                                      >  2:429966/1‑148 (MQ=255)
cGCTGATTTTTTCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCa                                                                                                                                                      >  1:468658/1‑148 (MQ=255)
      ttttttCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCaaaa                                                                                                                                                >  1:528169/1‑148 (MQ=255)
         tttCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                                                                            <  3:84392/149‑2 (MQ=255)
         tttCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                                                                            <  3:161498/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                                                                                                                                                     >  1:492862/1‑131 (MQ=255)
                  tCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCAc                                                                                                                                                     <  2:492862/131‑1 (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)
                     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)
                                gtAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTg                                                                                                                     <  4:66358/149‑1 (MQ=255)
                                 tAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTg                                                                                                                     <  2:883577/148‑1 (MQ=255)
                                   aCGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTg                                                                                                                  >  3:425668/1‑149 (MQ=255)
                                       ataATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGt                                                                                                               >  1:702012/1‑148 (MQ=255)
                                            gCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATc                                                                                                          >  2:871308/1‑148 (MQ=255)
                                              ggCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaa                                                                                                        <  1:91149/148‑1 (MQ=255)
                                              ggCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaa                                                                                                        <  2:6512/148‑1 (MQ=255)
                                              ggCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTACTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaa                                                                                                        <  2:689286/148‑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                                                                                                                        >  1:713204/1‑124 (MQ=255)
                                                      ttAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATc                                                                                                                        <  2:713204/124‑1 (MQ=255)
                                                          tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaaaaa                                                                                           <  4:233721/149‑1 (MQ=255)
                                                          tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaaaaa                                                                                           <  4:163489/149‑1 (MQ=255)
                                                          tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaaaa                                                                                            <  1:871308/148‑1 (MQ=255)
                                                            atcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCtt                                                                                                                      <  2:826714/120‑1 (MQ=255)
                                                            atcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCtt                                                                                                                      >  1:826714/1‑120 (MQ=255)
                                                            atcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAg                                                                                          <  1:58027/148‑1 (MQ=255)
                                                            atcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGc                                                                                         <  3:39263/149‑1 (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)
                                                                                   tATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGc                                                                  <  3:569485/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)
                                                                                           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)
                                                                                               cAGTAATTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAAcc                                                       >  2:446720/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)
                                                                                                   aTTTGGCAGTTATTGATATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGc                                                                                                                                             <  2:802109/57‑1 (MQ=37)
                                                                                                          aGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCt                                           <  4:238057/149‑1 (MQ=255)
                                                                                                                tGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGc                                                                                                                                             <  1:751312/45‑1 (MQ=37)
                                                                                                                 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                                                                                                                       <  2:759002/63‑1 (MQ=255)
                                                                                                                    tattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCt                                                                                                                       >  1:759002/1‑63 (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)
                                                                                                                               tCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTAt                       <  1:446720/148‑1 (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                            <  2:488921/133‑1 (MQ=255)
                                                                                                                                         gCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAg             >  2:54333/1‑148 (MQ=255)
                                                                                                                                          cTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGa                              >  4:41116/1‑130 (MQ=255)
                                                                                                                                            ttttGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAg          <  2:878013/148‑1 (MQ=255)
                                                                                                                                                   aCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAg                                                                   <  3:462124/84‑1 (MQ=255)
                                                                                                                                                   aCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAGGc        <  3:460743/143‑1 (MQ=255)
                                                                                                                                                   aCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAGGCGGATa   >  2:749359/1‑148 (MQ=255)
                                                                                                                                                   aCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAGGCGGATAc  >  3:361256/1‑149 (MQ=255)
                                                                                                                                                   aCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAGGCGGATAc  >  4:53399/1‑149 (MQ=255)
                                                                                                                                                   aCCAAAAAGCTCGGCTCCATCGGTTCCGATCTGGGTGCGTCGATCAAAGCCTTTAAAAAAGCAATGAGCGATGATGAACCAAAg                                                                   >  4:462124/1‑84 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCTGATTTTTTCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAGGCGGATAC  >  NC_000913/4021861‑4022156

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