Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,890,943 T→C 100% T109A (ACC→GCC)  tsaB ← tRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein; binding partner and protease for TsaD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,890,9430TC99.6% 412.4 / ‑6.5 112T109A (ACC→GCC) tsaBtRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein; binding partner and protease for TsaD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (57/54);  minor base A (0/1);  total (57/55)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 112 total aligned reads displayed.

GGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAA  >  NC_000913/1890797‑1891088
                                                                                                                                                  |                                                                                                                                                 
ggAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTg                                                                                                                                                 >  4:1392823/1‑149 (MQ=255)
 gAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTgc                                                                                                                                                <  2:1006791/149‑1 (MQ=255)
  aaaGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTgcg                                                                                                                                               >  3:462697/1‑149 (MQ=255)
  aaaGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTgc                                                                                                                                                <  2:526471/148‑1 (MQ=255)
  aaaGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTgc                                                                                                                                                >  1:526471/1‑148 (MQ=255)
  aaaGCTGTTGCATTCGATCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTgcg                                                                                                                                               <  2:934488/149‑1 (MQ=255)
      cTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGt                                                                                                                                           <  2:525071/149‑1 (MQ=255)
             aTTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTgcg                                                                                                                                    >  1:1863684/1‑149 (MQ=255)
              ttCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTgc                                                                                                                                                <  4:191626/136‑1 (MQ=255)
              ttCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTgc                                                                                                                                                >  3:191626/1‑136 (MQ=255)
                cGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCa                                                                                                                                 <  3:643036/149‑1 (MQ=255)
                 gTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCAc                                                                                                                                >  2:1760294/1‑149 (MQ=255)
                 gTTCATGGACGATTTCGGGCTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCAc                                                                                                                                >  2:727522/1‑149 (MQ=255)
                  ttCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACg                                                                                                                               >  3:348965/1‑149 (MQ=255)
                      tGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTgc                                                                                                                                     >  3:1836721/1‑139 (MQ=255)
                      tGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTgc                                                                                                                                     <  4:1836721/139‑1 (MQ=255)
                         aCGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTgc                                                                                                                                                >  3:191929/1‑125 (MQ=255)
                         aCGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCAc                                                                                                                                <  4:798499/141‑1 (MQ=255)
                         aCGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTg                                                                                                                        >  1:1534192/1‑149 (MQ=255)
                         aCGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGGAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTgc                                                                                                                                                <  4:191929/125‑1 (MQ=255)
                             tttCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCGGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCAc                                                                                                                                <  2:1812144/137‑1 (MQ=255)
                             tttCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTgcg                                                                                                                                               >  3:448171/1‑122 (MQ=255)
                             tttCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTgcg                                                                                                                                               <  4:448171/122‑1 (MQ=255)
                             tttCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCAc                                                                                                                                >  1:1812144/1‑137 (MQ=255)
                             tttCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTgcg                                                                                                                      <  4:1531782/147‑1 (MQ=255)
                             tttCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTgcg                                                                                                                      >  3:1531782/1‑147 (MQ=255)
                                    tttGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCa                                                                                                             >  4:567399/1‑149 (MQ=255)
                                      tGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCAtt                                                                                                           >  3:121201/1‑149 (MQ=255)
                                       gAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTa                                                                                                          <  1:727522/149‑1 (MQ=255)
                                       gAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTa                                                                                                          <  4:462697/149‑1 (MQ=255)
                                         gTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGc                                                                                                        <  2:1534192/149‑1 (MQ=255)
                                          tACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCg                                                                                                       >  3:1596420/1‑149 (MQ=255)
                                           aCGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGt                                                                                                      <  4:121201/149‑1 (MQ=255)
                                             ggCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTgg                                                                                                    <  1:1130144/149‑1 (MQ=255)
                                              gCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGcccc                                                                                                                           <  3:1178378/125‑1 (MQ=255)
                                              gCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGcccc                                                                                                                           >  4:1178378/1‑125 (MQ=255)
                                                ttCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGgaga                                                                                                 >  2:1250468/1‑149 (MQ=255)
                                                  cGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTgcgc                                                                                                                                              >  3:1383704/1‑102 (MQ=255)
                                                  cGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTgcgc                                                                                                                                              <  4:1383704/102‑1 (MQ=255)
                                                    gTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGcc                                                                                             >  2:14996/1‑149 (MQ=255)
                                                       tcttcACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATTGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAat                                                                                          <  2:1222004/149‑1 (MQ=255)
                                                        cttcACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAatc                                                                                         >  1:1102784/1‑149 (MQ=255)
                                                        cttcACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAatc                                                                                         <  2:1102784/149‑1 (MQ=255)
                                                            aCCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCg                                                                                     <  1:14996/149‑1 (MQ=255)
                                                                tGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTaa                                                                                 <  3:1662790/149‑1 (MQ=255)
                                                                  ccAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAAtt                                                                               <  1:1250468/149‑1 (MQ=255)
                                                                        aCCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGGGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCCTGGAGACGCCAATCATCGGTACTTCCgcgc                                                                         >  2:1373579/1‑149 (MQ=255)
                                                                        aCCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCgcgc                                                                         >  1:1780021/1‑149 (MQ=255)
                                                                          cGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCa                                                                       >  2:114459/1‑149 (MQ=255)
                                                                             tttCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAgc                                                                    <  1:1373579/149‑1 (MQ=255)
                                                                                        tGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCg                                                                                     >  4:1175287/1‑121 (MQ=255)
                                                                                        tGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCg                                                                                     <  1:1872045/121‑1 (MQ=255)
                                                                                        tGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCg                                                                                     <  3:1175287/121‑1 (MQ=255)
                                                                                        tGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTg                                                         >  2:840561/1‑149 (MQ=255)
                                                                                              tCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGAt                                                   <  1:1760294/149‑1 (MQ=255)
                                                                                                ggCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGc                                                                                              >  1:758355/1‑104 (MQ=255)
                                                                                                ggCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGc                                                                                              <  2:758355/104‑1 (MQ=255)
                                                                                                  cccAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCa                                               >  2:1112061/1‑149 (MQ=255)
                                                                                                    cAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCaa                                              >  4:352119/1‑148 (MQ=255)
                                                                                                    cAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCaa                                              <  3:352119/148‑1 (MQ=255)
                                                                                                    cAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAAt                                             >  4:538817/1‑149 (MQ=255)
                                                                                                    cAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAAt                                             >  2:167634/1‑149 (MQ=255)
                                                                                                     aGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATa                                            >  4:1715801/1‑149 (MQ=255)
                                                                                                      gTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATAc                                           >  2:1224933/1‑149 (MQ=255)
                                                                                                       tAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTa                                                                                                          >  3:835206/1‑85 (MQ=255)
                                                                                                          aCTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAAt                                       <  1:758837/149‑1 (MQ=255)
                                                                                                          aCTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAAt                                       <  2:1863684/149‑1 (MQ=255)
                                                                                                          aCTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAAt                                       >  4:1082871/1‑149 (MQ=255)
                                                                                                            ttCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATgc                                     >  3:1371310/1‑149 (MQ=255)
                                                                                                                cccATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAAt                                       >  4:614315/1‑143 (MQ=255)
                                                                                                                cccATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAAt                                       <  3:614315/143‑1 (MQ=255)
                                                                                                                 ccATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGgaga                                                                                                 <  2:1384155/84‑1 (MQ=255)
                                                                                                                 ccATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGgaga                                                                                                 >  1:1384155/1‑84 (MQ=255)
                                                                                                                 ccATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGAt                                                   <  3:1391281/130‑1 (MQ=255)
                                                                                                                 ccATTCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGAt                                                   >  4:1391281/1‑130 (MQ=255)
                                                                                                                    ttCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCa                                                                >  4:1856078/1‑114 (MQ=255)
                                                                                                                    ttCGCGCGTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCa                                                                <  3:1856078/114‑1 (MQ=255)
                                                                                                                       gcgcgTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTa                          >  1:1801569/1‑149 (MQ=255)
                                                                                                                         gcgTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTaaa                        <  2:1780021/149‑1 (MQ=255)
                                                                                                                          cgTCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAg                       <  3:538817/149‑1 (MQ=255)
                                                                                                                            tCAATGTCTGCCCGCACCCTATCTGCGCCGAATTTGCCCCCCGCACTTTGTGCCAGCACCTCTAGCGTGGAGATGTCTCTCAACAGCAACTCCGCGCCAAGCGCCAGCCCTTGCGCGGTGCTAATACCTATGCGCTCGCCAAGAAAGCt                     >  4:732620/1‑149 (MQ=255)
                                                                                                                            tCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAgcg                                                                   >  1:862912/1‑103 (MQ=255)
                                                                                                                            tCAATGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAgcg                                                                   <  2:862912/103‑1 (MQ=255)
                                                                                                                              aaTGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGTCCTTGCGCGATGCCAATACCAAt                                       <  2:265378/129‑1 (MQ=255)
                                                                                                                              aaTGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAAt                                       >  1:265378/1‑129 (MQ=255)
                                                                                                                              aaTGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTAc                   <  1:1224933/149‑1 (MQ=255)
                                                                                                                               aTGGCTGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTaaa                        >  2:205560/1‑143 (MQ=255)
                                                                                                                                    tGCCAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGc             >  1:824978/1‑149 (MQ=255)
                                                                                                                                    atccAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTaaa                        <  4:201947/136‑1 (MQ=255)
                                                                                                                                      ccAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTaa                                                                                 >  2:1749051/1‑79 (MQ=255)
                                                                                                                                      ccAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTaa                                                                                 <  1:1749051/79‑1 (MQ=255)
                                                                                                                                       cAGCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCcgc          >  4:270214/1‑149 (MQ=255)
                                                                                                                                         gCACGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCgcgc        >  2:1873950/1‑149 (MQ=255)
                                                                                                                                           aCGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGa                                                    <  3:1861670/103‑1 (MQ=255)
                                                                                                                                           aCGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGa                                                    >  4:1861670/1‑103 (MQ=255)
                                                                                                                                           aCGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGccc      >  1:1280759/1‑149 (MQ=255)
                                                                                                                                           aCGCGGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGccc      <  1:1873950/149‑1 (MQ=255)
                                                                                                                                            cgcgGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCgg                                                                                    >  2:317149/1‑70 (MQ=255)
                                                                                                                                            cgcgGGCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCgg                                                                                    <  1:317149/70‑1 (MQ=255)
                                                                                                                                               gggCTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTaa  <  3:567399/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAA  >  NC_000913/1890797‑1891088

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