Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,123,257 C→T 5.1% D422N (GAT→AAT)  cpsB ← mannose‑1‑phosphate guanyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,123,2570CT 5.1% 462.5 / 11.2 136D422N (GAT→AAT) cpsBmannose‑1‑phosphate guanyltransferase
Reads supporting (aligned to +/- strand):  ref base C (67/62);  new base T (3/4);  total (70/66)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

Only 100 of 138 total aligned reads displayed.

GCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATAGAGTCATATTTGCCCC  >  NC_000913/2123113‑2123403
                                                                                                                                                |                                                                                                                                                  
gCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAAt                                                                                                                                                 <  3:716318/148‑1 (MQ=255)
    cacaTCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTg                                                                                                                                             <  1:13363/148‑1 (MQ=255)
       atcatcCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATTGATGGTGata                                                                                                                                          <  4:311495/148‑3 (MQ=255)
               ctTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCaa                                                                                                                                                  <  1:466684/132‑1 (MQ=255)
                     gaTAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCcgc                                                                                                                            <  1:204010/148‑1 (MQ=255)
                       tAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGa                                                                                                                          >  1:494065/1‑148 (MQ=255)
                                   cgcACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCttt                                                                                                                                 >  1:756733/1‑129 (MQ=255)
                                        ttCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTcc                                                                                                         >  2:1078236/1‑148 (MQ=255)
                                         tCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATgg                                                                                                                                               >  2:715495/1‑109 (MQ=255)
                                         tCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATgg                                                                                                                                               <  1:715495/109‑1 (MQ=255)
                                         tCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTg                                                                                                                                      <  2:907715/118‑1 (MQ=255)
                                         tCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTg                                                                                                                                      >  1:907715/1‑118 (MQ=255)
                                             ttAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAacca                                                                                                                    <  4:401466/132‑1 (MQ=255)
                                                 aTCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGa                                                                                                <  1:706872/148‑1 (MQ=255)
                                                  tCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGAt                                                                                               <  1:1078236/148‑1 (MQ=255)
                                                   cGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATg                                                                                              <  1:635728/148‑1 (MQ=255)
                                                     aGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCaa                                                                                                                                                  >  3:408537/1‑94 (MQ=255)
                                                     aGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCaa                                                                                                                                                  <  4:408537/94‑1 (MQ=255)
                                                     aGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATgg                                                                                                                                               >  3:370990/1‑97 (MQ=255)
                                                                  ccGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGa                                                                               >  2:114818/1‑148 (MQ=255)
                                                                   cGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGa                                                                                                                          >  1:527783/1‑104 (MQ=255)
                                                                   cGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGa                                                                                                                          <  2:527783/104‑1 (MQ=255)
                                                                     ggTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACaa                                                                            >  3:1013243/1‑148 (MQ=255)
                                                                             aGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTcc                                                                                                                               >  2:891148/1‑89 (MQ=255)
                                                                             aGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTcc                                                                                                                               <  1:891148/89‑1 (MQ=255)
                                                                             aGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTc                                                                                                                                <  4:204845/88‑1 (MQ=255)
                                                                             aGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTc                                                                                                                                <  2:493490/88‑1 (MQ=255)
                                                                             aGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTc                                                                                                                                >  3:204845/1‑88 (MQ=255)
                                                                              ggCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTa                                                                                   >  2:539318/1‑132 (MQ=255)
                                                                              ggCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTa                                                                                   <  1:539318/132‑1 (MQ=255)
                                                                                cAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACttt                                                                                                                                        >  4:343742/1‑77 (MQ=255)
                                                                                   tGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGtt                                                              >  4:100203/1‑148 (MQ=255)
                                                                                            cccAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACttt                                                                                                                                        >  1:1006891/1‑65 (MQ=255)
                                                                                              cAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATg                                                   >  1:121583/1‑148 (MQ=255)
                                                                                               aGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGc                                                  >  4:90630/1‑148 (MQ=255)
                                                                                                    aaTATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTc                                             >  2:74872/1‑148 (MQ=255)
                                                                                                    aaTATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGATGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTg                                                      <  2:955390/139‑1 (MQ=255)
                                                                                                          aaTGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTa                                                                                   <  3:135134/104‑1 (MQ=255)
                                                                                                          aaTGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTa                                                                                   >  4:135134/1‑104 (MQ=255)
                                                                                                          aaTGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAg                                                                           >  2:872638/1‑112 (MQ=255)
                                                                                                          aaTGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAg                                                                           <  1:872638/112‑1 (MQ=255)
                                                                                                          aaTGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATg                                                   <  3:385032/136‑1 (MQ=255)
                                                                                                          aaTGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATg                                                   >  4:385032/1‑136 (MQ=255)
                                                                                                           aTGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTg                                                                                                 >  1:1089307/1‑89 (MQ=255)
                                                                                                            tGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTa                                     >  2:1001770/1‑148 (MQ=255)
                                                                                                              gACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGc                                   <  4:1013243/148‑1 (MQ=255)
                                                                                                               aCTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCg                                                 >  1:20401/1‑133 (MQ=255)
                                                                                                               aCTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCg                                                 <  2:20401/133‑1 (MQ=255)
                                                                                                                cTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCg                                                 <  3:131109/132‑1 (MQ=255)
                                                                                                                cTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCg                                                 >  4:131109/1‑132 (MQ=255)
                                                                                                                 tCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCAccc                                                                                                                 >  3:569575/1‑67 (MQ=255)
                                                                                                                 tCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCAccc                                                                                                                 <  4:569575/67‑1 (MQ=255)
                                                                                                                  cGTTTTCCCCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGc                                                                          <  2:1052247/105‑1 (MQ=255)
                                                                                                                    ttttCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCcgcg                                                                                                      <  4:56534/75‑1 (MQ=255)
                                                                                                                    ttttCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGc                             >  1:800520/1‑148 (MQ=255)
                                                                                                                     tttCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTgg                                       <  2:1091487/137‑1 (MQ=255)
                                                                                                                     tttCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTgg                                       >  1:1091487/1‑137 (MQ=255)
                                                                                                                     tttCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGcc                            >  1:74963/1‑148 (MQ=255)
                                                                                                                       tCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCcg                          >  1:903102/1‑148 (MQ=255)
                                                                                                                          ccAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGccc                                                                        >  1:1000738/1‑99 (MQ=255)
                                                                                                                          ccAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGccc                                                                        <  2:1000738/99‑1 (MQ=255)
                                                                                                                            aaGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGa                                                     >  3:419317/1‑116 (MQ=255)
                                                                                                                             aGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCa                                            <  3:934004/124‑1 (MQ=255)
                                                                                                                              gCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTc                                             >  1:31144/1‑122 (MQ=255)
                                                                                                                              gCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTc                                             <  2:31144/122‑1 (MQ=255)
                                                                                                                                 gTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATaga                <  2:903102/148‑1 (MQ=255)
                                                                                                                                 gTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATaga                <  3:100203/148‑1 (MQ=255)
                                                                                                                                   ttGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATTGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACg                                                         <  4:943803/105‑1 (MQ=255)
                                                                                                                                   ttGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGa                                                                               <  3:805250/83‑1 (MQ=255)
                                                                                                                                   ttGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGa                                                                               >  4:805250/1‑83 (MQ=255)
                                                                                                                                   ttGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTc                                                            <  3:1143795/102‑1 (MQ=255)
                                                                                                                                   ttGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTc                                                            >  4:1143795/1‑102 (MQ=255)
                                                                                                                                   ttGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACg                                                         >  3:943803/1‑105 (MQ=255)
                                                                                                                                   ttGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATAGAGt              >  2:397261/1‑148 (MQ=255)
                                                                                                                                   ttGATATCACCATCAATGGTGACTTTTGCCGTCCCCGCGTCAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCGGGTTTCACg                                                         <  4:738704/105‑1 (MQ=255)
                                                                                                                                      atatCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCa                                            <  2:217608/115‑1 (MQ=255)
                                                                                                                                      atatCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCa                                            >  1:217608/1‑115 (MQ=255)
                                                                                                                                       tatCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATAGAGTCata          >  4:953922/1‑148 (MQ=255)
                                                                                                                                       agtCACCATTGATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATagag               >  2:151850/3‑143 (MQ=255)
                                                                                                                                         tCACCATTGATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATaga                <  1:488077/140‑1 (MQ=255)
                                                                                                                                         tCACCATTGATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATaga                >  2:488077/1‑140 (MQ=255)
                                                                                                                                         tCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCAcc                                          <  4:201295/114‑1 (MQ=255)
                                                                                                                                         tCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCAcc                                          >  3:201295/1‑114 (MQ=255)
                                                                                                                                         tCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCAc                                           <  3:930422/113‑1 (MQ=255)
                                                                                                                                         tCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCAc                                           >  4:930422/1‑113 (MQ=255)
                                                                                                                                          cACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGc                                   <  1:1040894/120‑1 (MQ=255)
                                                                                                                                          cACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGc                                   >  2:1040894/1‑120 (MQ=255)
                                                                                                                                          cACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATag                 >  3:479718/1‑138 (MQ=255)
                                                                                                                                          cACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATag                 <  4:479718/138‑1 (MQ=255)
                                                                                                                                          cACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATAGAGTCata          <  3:890007/145‑1 (MQ=255)
                                                                                                                                          cACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATAGAGTCata          >  4:890007/1‑145 (MQ=255)
                                                                                                                                          cACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATAGAGTCATAttt       >  2:762830/1‑148 (MQ=255)
                                                                                                                                          cACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATAGAGTCATAttt       >  3:612346/1‑148 (MQ=255)
                                                                                                                                           aCCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATAGAGTCATATTTg      <  1:995108/148‑1 (MQ=255)
                                                                                                                                            ccATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATaga                <  2:59505/137‑1 (MQ=255)
                                                                                                                                            ccATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATAGAGTCATATTTGc     >  1:279997/1‑148 (MQ=255)
                                                                                                                                             cATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTg                                        >  2:380714/1‑112 (MQ=255)
                                                                                                                                             cATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTg                                        <  1:380714/112‑1 (MQ=255)
                                                                                                                                                cAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCgg                                                                >  4:228943/1‑85 (MQ=255)
                                                                                                                                                cAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCgg                                                                <  3:228943/85‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACCCAGTGTTCCGCGCGGTGATGGTGCATCTGTACCGACAAGCCCTCGCCCGGTTTCACGGTGATGCGTTTCACCTGGTAGCGGTCGCCCGCGTCGATAGAGTCATATTTGCCCC  >  NC_000913/2123113‑2123403

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