Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,569,910 C→A 7.7% pseudogene (912/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,9100CA 7.7% 165.6 / 5.1 61pseudogene (912/921 nt)yjiPpseudogene, transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base C (32/24);  new base A (3/2);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-01

ATAAAAGGGGAGCAGCGCATTCTTCGATTGTTGTTGCAGAATGGGGCGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTACATAACCCATTAAATATATTGAACTTTAAAGATTTTTGTAGACCTGGTCAGGCGTTCACATGGCATCCGGCAT  >  NC_000913/4569767‑4570047
                                                                                                                                               |                                                                                                                                         
ataAAAGGGGAGCAGCGCATTCTTCGATTGTTGTTGCAGAATGGGGCGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAGGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAagag                                                                                                                                      >  2:637693/1‑149 (MQ=255)
ataAAAGGGGAGCAGCGCATTCTTCGATTGTTGTTGCAGAATGGGGCGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAagag                                                                                                                                      <  2:358698/149‑1 (MQ=255)
   aaaGGGGAGCAGCGCATTCTTCGATTGTTGTTGCAGAATGGGGCGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTa                                                                                                                                   >  1:792406/1‑149 (MQ=255)
         gAGCAGCGCATTCTTCGATTGTTGTTGCAGAATGGGGCGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCaga                                                                                                                             >  1:73622/1‑149 (MQ=255)
          agcagcGCATTCTTCGATTGTTGTTGCAGAATGGGGCGGGTCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATCGCTCAAGAGCTAATCagag                                                                                                                            >  2:721300/1‑149 (MQ=255)
           gcagcGCATTCTTCGATTGTTGTTGCAGAATGGGGCGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTTAGCGCTATTCATGGCTCAAGAGCTAATCAgaga                                                                                                                           <  2:554545/149‑1 (MQ=255)
           gcagcGCATTCTTCGATTGTTGTTGCAGAATGGGGCGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGAGAACGCTATTCATGGCTCAAGAGCTAATCaaaaa                                                                                                                           >  1:332560/1‑145 (MQ=255)
               cgcATTCTTCGATTGTTGTTGCAGAATGGGGCGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGa                                                                                                                       >  1:135300/1‑149 (MQ=255)
                  aTTCTTCGATTGTTGTTGCAGAATGGGGCGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGa                                                                                                                    <  1:349082/149‑1 (MQ=255)
                      ttcGATTGTTGTTGCAGAATGGGGCGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGAGaaa                                                                                                                >  1:641799/1‑149 (MQ=255)
                          aTTGTTGTTGCAGAATGGGGCGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCa                                                                                                            >  1:592510/1‑149 (MQ=255)
                            tgttgttgCAGAATGGGGCGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGGGACGGATGACAAACGCaaa                                                                                                          >  1:58438/1‑149 (MQ=255)
                             gttgttgCAGAATGGGGCGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGGCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAag                                                                                                         <  2:624771/149‑1 (MQ=255)
                                 ttgCAGAATGGGGCGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAagcagc                                                                                                     >  2:729103/1‑149 (MQ=255)
                                   gCAGAATGGGGCGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCt                                                                                                   <  1:156836/149‑1 (MQ=255)
                                        aTGGTGCCGATCCTAGATCGATACAAAAGATTACCCGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATgc                                                                                              <  1:626731/149‑1 (MQ=255)
                                           gggCGGATCCTGACTGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCt                                                                                           >  2:492025/1‑149 (MQ=255)
                                            ggCGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTa                                                                                          <  2:683625/149‑1 (MQ=255)
                                            ggCGGATCCTGAATGGATACAAAAGATTACCACACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCGATTCATGGCTCAAGAGCTAATCAGAGACGGAAGCCAAACGCACAGCAGCGTAAAGCGCTa                                                                                          >  2:371114/1‑149 (MQ=255)
                                              cGGATCCTGAATGGATACAAAAGGTTACCGGACTATCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGCGAGCGCTATCCATGGCTCAAGAGCTAATCAGAGACGGCTCACACTGTCAAAGCAGCCTGATGCGCTATg                                                                                        <  2:405652/149‑1 (MQ=255)
                                              cGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATg                                                                                        <  1:359012/149‑1 (MQ=255)
                                                 aTCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGGCTTATTCGCCATgtat                                                                                     >  2:198554/1‑146 (MQ=255)
                                                   ccTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTAt                                                                                   <  1:637693/149‑1 (MQ=255)
                                                   ccTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAACGTGAGCGCTATTCATGGCTAAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTAt                                                                                   >  2:75044/1‑149 (MQ=255)
                                                   ccTGAAAGGATACAAGAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTAAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTAt                                                                                   >  2:470152/1‑149 (MQ=255)
                                                     tGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCa                                                                                 <  1:245369/149‑1 (MQ=255)
                                                      gAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAg                                                                                <  1:89470/149‑1 (MQ=255)
                                                        aTGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGTAGCCCGTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGt                                                                              <  2:483067/149‑1 (MQ=255)
                                                                aaaaGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGACTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGACTACATaa                                                                      >  1:259813/1‑149 (MQ=255)
                                                                 aaaGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTTTCAGGGCTACATAAc                                                                     >  1:45278/1‑149 (MQ=255)
                                                                  aaGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGCCTACATAAcc                                                                    >  2:594185/1‑149 (MQ=255)
                                                                  aaGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGACTACATAAcc                                                                    >  1:737912/1‑149 (MQ=255)
                                                                  aaGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCGCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTACATAAcc                                                                    <  2:694773/149‑1 (MQ=255)
                                                                       tACCGGACTGTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGACGCGCTATGCTTATCAGGTCTACATAACCCATTa                                                               <  2:568823/149‑1 (MQ=255)
                                                                           ggACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAACCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGGCCACATAACCCATTAAata                                                           >  2:498234/1‑149 (MQ=255)
                                                                               tttCGGCAGAGCAAATGCAGGCATTAAGACAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGACTACGTAACCCATTAAATATAtt                                                       >  1:427436/1‑149 (MQ=255)
                                                                                ttCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTCCATAACCCATTAAATATATTg                                                      >  1:342576/1‑149 (MQ=255)
                                                                                ttCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGGCTACATAACCCATTAAATATATTg                                                      >  1:742181/1‑149 (MQ=255)
                                                                                   ggCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTCATCAGGTCTCCATAACCCATTAAATATATTGAAc                                                   <  1:402482/149‑1 (MQ=255)
                                                                                    gCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATAGCTCAAGAGCTAATAAGAGACGGATGACACCCAGGAAACATCGGTATGCGCCATGCCTATCATGGCTACATAACCCACCAAATGGCTTGAACt                                                  >  1:382341/1‑149 (MQ=255)
                                                                                        agCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTACATAACCCATTAAATATATTGAACTTTaa                                              >  1:527855/1‑149 (MQ=255)
                                                                                        agCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGAAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTACATAACCCATTAAATATATTGAACTTTaa                                              <  2:73622/149‑1 (MQ=255)
                                                                                        agCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGATATGCATGGCTAAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGGCTACCTAACCCATTAGATATATTGAACTTTaa                                              >  1:387086/1‑149 (MQ=255)
                                                                                         gCAAATGCCGGCATTAAGGCCGCCCTTGCCTGAGCGTCAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTACATAACCCATTAAATATATTGAACTTTaaa                                             <  2:479239/149‑1 (MQ=255)
                                                                                             aTGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTACATAACCCATAAAATATATTGAACGTTAAAGAtt                                         >  1:697565/1‑149 (MQ=255)
                                                                                                      ttAAGGCAGCCCTTGCCTGGGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTACATAACCCATTAAATATATTGAACTTTAAAGATTTTTGTAGAc                                <  2:429483/149‑1 (MQ=255)
                                                                                                          ggCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTACATAACCCATTAAATATATTGAACTTTAAAGATTTTTGTAGACCTgg                            >  2:744277/1‑149 (MQ=255)
                                                                                                           gCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAGACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTACATAACCCATTAAATATATTGAACTTTAAAGATTTTTGTAGGCCTGGt                           >  1:538455/1‑149 (MQ=255)
                                                                                                           gCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGGCTACATAACCCATTAAATATATTGAACTTTAAAGCCTTTTGTAGACCTGGt                           >  1:301783/1‑149 (MQ=255)
                                                                                                                ccTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTACATAACCCATTAAATATATTGAACTTTAAAGATTTTTGTAGACCTGGTCAGGc                      <  1:570818/149‑1 (MQ=255)
                                                                                                                    gCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTACATAACCCATTAAATATATTGAACTTTAAAGATTTTTGTAGACCTGGTCAGGCGTTc                  <  2:135300/149‑1 (MQ=255)
                                                                                                                     ccTGAGCGTGAGCGCTATTCATGGCTAAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTACATAACCCATTAAATATATTGAACTTTAAAGATTTTTGTAGACCTGGTCAGGCGTTca                 <  1:470152/149‑1 (MQ=255)
                                                                                                                     ccTGAACGTGAGCGCTATTCATGGCTAAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTACATAACCCATTAAATATATTGAACTTTAAAGATTTTTGTAGACCTGGTCAGGCGTTca                 <  1:75044/149‑1 (MQ=255)
                                                                                                                        gcgcgtgagcgCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTACATAACCCATTAAATATATTGAACTTTAAAGATTTTTGTAGACCTGGTCAGGCGTTCACAt              <  2:710622/147‑1 (MQ=255)
                                                                                                                        gagcgtgagcgCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGGCTACATAACCCATTAAATATATTGACCTTTAAAGATTTTTGTAGACCTGGTCAGGCGTTCACAt              >  2:259015/1‑149 (MQ=255)
                                                                                                                           cgtgagcgCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTACATAACCCATTAAATATATTGAACTTTAAAGATTTTTGTAGACCTGGTCAGGCGTTCACATGGc           <  2:58438/149‑1 (MQ=255)
                                                                                                                               agcgCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGGCTACATAACCCATTAAATATATTGAACTTTAAAGATTTTTGTAGACCTGGTCAGGCGTTCACATGGCATcc       >  1:682103/1‑149 (MQ=255)
                                                                                                                               agcgCTATTCATGGCTCAAGAGCGAATCAGAGACTGATGACAAACGCAAAGCAACCTGATGCGCTATGCTTATCAGGTCTACATAACCCATTAAATATATTGAACTTTAAAGATTTTTGTAGACCTGGTCAGGCGTTCAAATGGCATcc       >  1:439063/1‑149 (MQ=255)
                                                                                                                                   cTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTACATAACCCATTAAATATATTGAACTTTAAAGATTTTTGTAGACCTGGTCAGGCGTTCACATGGCATCCGGCa   <  1:717995/149‑1 (MQ=255)
                                                                                                                                    tATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGGCTACATAACCCATTAAATATATTGAACTTTAAAGATTTTTGTATACCTGGTCAGGCGTTCACATGGCATCCGGCAt  >  1:290453/1‑149 (MQ=255)
                                                                                                                                    tATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGGCTACATAACCCATTAAATATATTGAACTTTAAAGATTTTTGTAGACCTGGTCAGGCGTTCAGATGGCATCCGGCAt  >  1:463881/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                         
ATAAAAGGGGAGCAGCGCATTCTTCGATTGTTGTTGCAGAATGGGGCGGATCCTGAATGGATACAAAAGATTACCGGACTTTCGGCAGAGCAAATGCAGGCATTAAGGCAGCCCTTGCCTGAGCGTGAGCGCTATTCATGGCTCAAGAGCTAATCAGAGACGGATGACAAACGCAAAGCAGCCTGATGCGCTATGCTTATCAGGTCTACATAACCCATTAAATATATTGAACTTTAAAGATTTTTGTAGACCTGGTCAGGCGTTCACATGGCATCCGGCAT  >  NC_000913/4569767‑4570047

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑