Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,272,378 C→A 5.7% A498A (GCG→GCT uvrA ← ATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,272,3780CA 5.7% 552.2 / 15.8 157A498A (GCG→GCTuvrAATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC
Reads supporting (aligned to +/- strand):  ref base C (81/67);  new base A (3/6);  total (84/73)
Fisher's exact test for biased strand distribution p-value = 3.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 157 total aligned reads displayed.

CGTGCTCCACCACAATCACGGTATTACCGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTCTGCAATCTTCGCCCGCTGACCT  >  NC_000913/4272233‑4272525
                                                                                                                                                 |                                                                                                                                                   
cgtGCTCCACCACAATCACGGTATTACCGGGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGcc                                                                                                                                                  <  4:1250439/149‑1 (MQ=255)
     tCCACCACAATCACGGTATTACCGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACg                                                                                                                                             >  2:993336/1‑149 (MQ=255)
     tCCACCACAATCACGGTATTACCGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCGAGCCAGatt                                                                                                                                             <  3:833433/149‑3 (MQ=255)
     tCCACCACAATCACGGTATTACCGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCGAGCCAGatt                                                                                                                                             <  2:413889/149‑3 (MQ=255)
      ccaccaCAATCACGGTATTACCGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCGAGCCAGattg                                                                                                                                            <  1:1198910/149‑4 (MQ=255)
      ccaccaCAATCACGGTATTACCGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCGAGCCAGattg                                                                                                                                            <  1:44573/149‑4 (MQ=255)
       caccacAATCACGGTATTACCGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGa                                                                                                                                           >  1:687556/1‑149 (MQ=255)
         ccacAATCACGGTATTACCGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATa                                                                                                                                         <  3:1500053/149‑1 (MQ=255)
             aaTCACGGTATTACCGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCt                                                                                                                                     <  1:926936/149‑1 (MQ=255)
             aaTCACGGTATTACCGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCt                                                                                                                                     <  2:1053549/149‑1 (MQ=255)
               tCACGGTATTACCGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCtgt                                                                                                                                   >  3:45402/1‑149 (MQ=255)
                 aCGGTATTACCGAGATCGCGCAGATGGATAAGCGTACTCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGc                                                                                                                                 <  4:1369989/149‑1 (MQ=255)
                 aCGGTATTACCGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGc                                                                                                                                 >  2:921037/1‑149 (MQ=255)
                 aCGGTATTACCGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGc                                                                                                                                 <  3:1133606/149‑1 (MQ=255)
                     tATTACCGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCGAGCCAGatt                                                                                                                                             <  2:137029/133‑3 (MQ=255)
                      aTTACCGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCAc                                                                                                                            <  4:714045/149‑1 (MQ=255)
                         aCCGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGc                                                                                                                         >  3:753537/1‑149 (MQ=255)
                          ccGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGcc                                                                                                                        >  3:1061264/1‑149 (MQ=255)
                                 cgcgCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGc                                                                                                                 <  4:235060/149‑1 (MQ=255)
                                 cgcgCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGc                                                                                                                 <  2:529468/149‑1 (MQ=255)
                                   cgcAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGt                                                                                                               <  1:1236239/149‑1 (MQ=255)
                                      aGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTc                                                                                                            >  4:527522/1‑149 (MQ=255)
                                        aTGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGAt                                                                                                                                          <  2:207697/117‑1 (MQ=255)
                                        aTGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGAt                                                                                                                                          >  1:207697/1‑117 (MQ=255)
                                                gCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCATGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCgg                                                                                                  >  4:1485154/1‑149 (MQ=255)
                                                gCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGt                                                                                                               >  1:1188245/1‑136 (MQ=255)
                                                gCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGt                                                                                                               <  2:1188245/136‑1 (MQ=255)
                                                gCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCgg                                                                                                  >  4:634650/1‑149 (MQ=255)
                                                gCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCgg                                                                                                  >  1:422924/1‑149 (MQ=255)
                                                   tACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGaa                                                                                               >  1:1063624/1‑149 (MQ=255)
                                                    aCCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGgaa                                                                                              <  3:1485154/149‑1 (MQ=255)
                                                      ccAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGc                                                                                            >  2:1002927/1‑149 (MQ=255)
                                                       cAACGGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGCATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCg                                                                                           >  4:407053/1‑149 (MQ=255)
                                                       cAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCa                                                                                                                       <  4:437535/121‑1 (MQ=255)
                                                       cAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCg                                                                                           >  1:1585264/1‑149 (MQ=255)
                                                       cAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCg                                                                                           >  3:507697/1‑149 (MQ=255)
                                                       cAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAACCGTTTCTGCCg                                                                                                       >  4:591262/1‑137 (MQ=255)
                                                        aaCAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGt                                                                                          >  4:1045608/1‑149 (MQ=255)
                                                           aGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGc                                                                                                                                                   <  3:425527/89‑1 (MQ=255)
                                                           aGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGc                                                                                                                                                   >  4:425527/1‑89 (MQ=255)
                                                              cgcTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTa                                                                                    >  4:196973/1‑149 (MQ=255)
                                                               gcTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTaa                                                                                   >  4:1560348/1‑149 (MQ=255)
                                                               gcTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTaa                                                                                   <  2:226843/149‑1 (MQ=255)
                                                                 tCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAAtt                                                                                 >  3:272101/1‑149 (MQ=255)
                                                                 tCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAAtt                                                                                 >  1:1008688/1‑149 (MQ=255)
                                                                   gTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCa                                                                               <  4:843374/149‑1 (MQ=255)
                                                                    ttATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAg                                                                              >  4:1211598/1‑149 (MQ=255)
                                                                       tCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGc                                                                                            >  2:1151882/1‑132 (MQ=255)
                                                                         aCGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGa                                                                         <  3:989588/149‑1 (MQ=255)
                                                                               gTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTaa                                                                   <  4:465409/149‑1 (MQ=255)
                                                                                  cAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTAGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGa                                                                <  1:909408/149‑1 (MQ=255)
                                                                                     gCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACCCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGa                                                             >  4:905726/1‑149 (MQ=255)
                                                                                      ccGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGaa                                                            <  2:1585264/149‑1 (MQ=255)
                                                                                      ccGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGaa                                                            <  1:1415576/149‑1 (MQ=255)
                                                                                           aGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCa                                                       <  4:791217/149‑1 (MQ=255)
                                                                                             aCGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGa                                                     <  4:185086/149‑1 (MQ=255)
                                                                                               ggCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACg                                                   >  3:1017348/1‑149 (MQ=255)
                                                                                                 cTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGAt                                                 <  1:387992/149‑1 (MQ=255)
                                                                                                  tCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATc                                                <  1:993336/149‑1 (MQ=255)
                                                                                                     tCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGcc                                             >  3:1212874/1‑149 (MQ=255)
                                                                                                     tCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGcc                                             >  1:712933/1‑149 (MQ=255)
                                                                                                       cAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGa                                           <  3:196765/149‑1 (MQ=255)
                                                                                                       cAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGa                                           >  1:525184/1‑149 (MQ=255)
                                                                                                       cAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGa                                           <  2:620397/149‑1 (MQ=255)
                                                                                                           aCGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGAtctc                                       <  4:350013/149‑1 (MQ=255)
                                                                                                           aCGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGAtctc                                       >  1:897986/1‑149 (MQ=255)
                                                                                                               aCATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTa                                   >  3:1136390/1‑149 (MQ=255)
                                                                                                                 aTAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAg                                 <  3:367704/149‑1 (MQ=255)
                                                                                                                  tAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAgg                                >  3:714053/1‑149 (MQ=255)
                                                                                                                   aaCGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGaaa                                                                                              <  4:169545/86‑1 (MQ=255)
                                                                                                                   aaCGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGaaa                                                                                              >  3:169545/1‑86 (MQ=255)
                                                                                                                      gCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGAttt                            >  3:280580/1‑149 (MQ=255)
                                                                                                                      gCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGAttt                            >  2:764525/1‑149 (MQ=255)
                                                                                                                       ccAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGAtttt                           <  1:1350134/149‑1 (MQ=255)
                                                                                                                       ccAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGAtttt                           <  3:181734/149‑1 (MQ=255)
                                                                                                                         aaCCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGAtttttt                         >  1:561961/1‑149 (MQ=255)
                                                                                                                         aaCCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGAtttttt                         >  2:832848/1‑149 (MQ=255)
                                                                                                                          aCCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGGCGATCGCCGATCTCTTTAAGGATTTTTTc                        <  3:1032817/149‑1 (MQ=255)
                                                                                                                          aCCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTc                        >  4:1264172/1‑149 (MQ=255)
                                                                                                                          aCCAGGCCCGCACCAATCCGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTc                        <  2:268013/149‑1 (MQ=255)
                                                                                                                           ccAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTCt                       >  2:1595092/1‑149 (MQ=255)
                                                                                                                            cAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTCTg                      >  4:736026/1‑149 (MQ=255)
                                                                                                                            cAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTCTg                      <  2:668436/149‑1 (MQ=255)
                                                                                                                             aGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTCTGc                     >  4:468477/1‑149 (MQ=255)
                                                                                                                                cccGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAAttt                                                         >  1:24965/1‑110 (MQ=255)
                                                                                                                                  cGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTCTGCAATCt                >  2:553439/1‑149 (MQ=255)
                                                                                                                                   gCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTCTGCAATCtt               >  1:32752/1‑149 (MQ=255)
                                                                                                                                    cACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTCTGCAATCTTc              >  3:881947/1‑149 (MQ=255)
                                                                                                                                    cACCAATCTGGCTCGCCAGACGGATACCCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTCTGCAATCTTc              <  3:905726/149‑1 (MQ=255)
                                                                                                                                     aCCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGCAAGCGTTTCTGACGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGaa                                                            >  4:342797/1‑102 (MQ=255)
                                                                                                                                     aCCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTCTGCAATCTTCg             <  1:1144659/149‑1 (MQ=255)
                                                                                                                                      ccAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTCTGCAATCTTCGc            <  1:1006872/149‑1 (MQ=255)
                                                                                                                                          tCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGa                                                                >  4:522888/1‑93 (MQ=255)
                                                                                                                                          tCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGa                                                                <  3:522888/93‑1 (MQ=255)
                                                                                                                                          tCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTCTGCAATCTTCGCCCGc        <  2:1063624/149‑1 (MQ=255)
                                                                                                                                          tCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTCTGCAATCTTCGCCCGc        >  4:746677/1‑149 (MQ=255)
                                                                                                                                            tGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGaaa                                                                                              <  2:332767/61‑1 (MQ=255)
                                                                                                                                             ggCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTCTGCAATCTTCGCCCGCTGa     >  3:1388493/1‑149 (MQ=255)
                                                                                                                                                tCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTCTGCAATCTTCGCCCGCTGACCt  >  4:1206219/1‑149 (MQ=255)
                                                                                                                                                tCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTCTGCAATCTTCGCCCGCTGACCt  >  2:906008/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGTGCTCCACCACAATCACGGTATTACCGAGATCGCGCAGATGGATAAGCGTACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGCGTCAGGTAATTCAGGCCGACGTTAACGAGGAATTTCAGACGATCGCCGATCTCTTTAAGGATTTTTTCTGCAATCTTCGCCCGCTGACCT  >  NC_000913/4272233‑4272525

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: