Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,825,400 T→C 7.4% V64A (GTG→GCG)  recG → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,825,4000TC 7.4% 8653.7 / inf 2748V64A (GTG→GCG) recGATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base T (1317/1222);  new base C (107/97);  total (1426/1322)
Fisher's exact test for biased strand distribution p-value = 8.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.68e-01

Only 100 of 2782 total aligned reads displayed.

ACGGGCGTTGGCGCAGCACTTAGTAACAAACTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTA  >  NC_000913/3825252‑3825548
                                                                                                                                                    |                                                                                                                                                    
  gggCGTTGGCGCAGCACTTAGTAACAAACTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTgg                                                                                                                                                    >  4:2266797/1‑149 (MQ=255)
  gggCGTTGGCGCAGCACTTAGTAACAAACTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTgg                                                                                                                                                    <  3:1465041/149‑1 (MQ=255)
  gggCGTTGGCGCAGCACTTAGTAACAAACTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTgg                                                                                                                                                    >  1:4928916/1‑149 (MQ=255)
  gggCGTTGGCGCAGCACTTAGTAACAAACTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTgg                                                                                                                                                    >  4:3370983/1‑149 (MQ=255)
    gCGTTGGCGCAGCACTTAGTAACAAACTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGaa                                                                                                                                                  <  4:6208041/149‑1 (MQ=255)
           cgcAGCACTTAGTAACAAACTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAg                                                                                                                                           >  3:4745838/1‑149 (MQ=255)
            gcagcaCTTAGTAACAAACTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGt                                                                                                                                          >  1:6031628/1‑149 (MQ=255)
             cagcaCTTAGTAACAAACTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTACCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTg                                                                                                                                         <  4:4040229/149‑1 (MQ=255)
              agcaCTTAGTAACAAACTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGc                                                                                                                                        <  3:253562/149‑1 (MQ=255)
                  cTTAGTAACAAACTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGCGGAAGGCGAAGTGCTGaa                                                                                                                                    <  3:6295070/149‑1 (MQ=255)
                     aGTAACAAACTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTg                                                                                                                                 >  3:7055803/1‑149 (MQ=255)
                      gTAACAAACTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGc                                                                                                                                <  1:723004/149‑1 (MQ=255)
                      gTAACAAACTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGc                                                                                                                                <  1:6480585/149‑1 (MQ=255)
                        aaCAAACTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCaa                                                                                                                              >  1:5168041/1‑149 (MQ=255)
                         aCAAACTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGCGGAAGGCGAAGTGCTGAACTGCAat                                                                                                                             >  1:6950672/1‑149 (MQ=255)
                           aaaCTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAatat                                                                                                                           <  1:3202217/149‑1 (MQ=255)
                            aaCTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAgg                                                                                                                                                <  3:4073458/127‑1 (MQ=255)
                            aaCTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATAtc                                                                                                                          >  2:4426193/1‑149 (MQ=255)
                               tGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTcc                                                                                                                       >  3:3839685/1‑149 (MQ=255)
                               tGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTcc                                                                                                                       <  2:5444885/149‑1 (MQ=255)
                                  cGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTtc                                                                                                                    <  1:547123/149‑1 (MQ=255)
                                   gAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTtcg                                                                                                                   >  1:6229766/1‑149 (MQ=255)
                                   gAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTACCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTtcg                                                                                                                   >  3:282936/1‑149 (MQ=255)
                                     aaaTCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTtcggc                                                                                                                 <  2:2551294/149‑1 (MQ=255)
                                         cAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTtcggcggtc                                                                                                             <  1:2992669/149‑1 (MQ=255)
                                           aCCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGACGTGCTGAACTGCAATATCTCCTtcggcggtcgg                                                                                                           >  2:1976424/1‑149 (MQ=255)
                                           aCCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTtcggcggtcgg                                                                                                           >  3:941319/1‑149 (MQ=255)
                                            ccTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTtcggcggtcggc                                                                                                          >  4:4908972/1‑149 (MQ=255)
                                             cTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTtcggcggtcggcg                                                                                                         >  4:5947906/1‑149 (MQ=255)
                                               gCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGga                                                                                                       <  3:3633582/149‑1 (MQ=255)
                                                cATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGgat                                                                                                      <  4:3143527/149‑1 (MQ=255)
                                                 aTACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGgatg                                                                                                     <  2:3537415/149‑1 (MQ=255)
                                                   aCCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGgatgat                                                                                                   <  3:3370983/149‑1 (MQ=255)
                                                    ccGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGatgatg                                                                                                  <  2:6275687/149‑1 (MQ=255)
                                                      gTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAgg                                                                                                                                                <  2:1276537/101‑1 (MQ=255)
                                                      gTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGAc                                                                                                <  1:674542/149‑1 (MQ=255)
                                                      gTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGAc                                                                                                <  4:2584807/149‑1 (MQ=255)
                                                       tGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGAcc                                                                                               >  1:42340/1‑149 (MQ=255)
                                                            gATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGgat                                                                                                      <  1:4646078/137‑1 (MQ=255)
                                                            gATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGgat                                                                                                      <  4:5680270/137‑1 (MQ=255)
                                                             aTCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAg                                                                                         <  2:2859217/149‑1 (MQ=255)
                                                             aTCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAg                                                                                         >  3:1803535/1‑149 (MQ=255)
                                                                  ctctTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTtc                                                                                                                    >  2:679727/1‑117 (MQ=255)
                                                                  ctctTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAg                                                                                    <  3:1220058/149‑1 (MQ=255)
                                                                   tctTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGc                                                                                   >  1:1014536/1‑149 (MQ=255)
                                                                     ttACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTTCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGa                                                                                 >  3:5806728/1‑149 (MQ=255)
                                                                     ttACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGaa                                                                                                                                            >  1:1292961/1‑90 (MQ=255)
                                                                     ttACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGc                                                                                                                                <  3:6781959/102‑1 (MQ=255)
                                                                     ttACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGa                                                                                 <  2:4928916/149‑1 (MQ=255)
                                                                      tACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGAc                                                                                <  1:2168906/149‑1 (MQ=255)
                                                                         acCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGCGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGt                                                                             >  1:3955739/1‑149 (MQ=255)
                                                                           cTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGcc                                                                                           <  1:3460578/133‑1 (MQ=255)
                                                                           cTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGAc                                                                                >  4:745185/1‑144 (MQ=255)
                                                                             tCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCg                                                                         >  2:701032/1‑149 (MQ=255)
                                                                              ccccTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCgg                                                                        <  3:2315864/149‑1 (MQ=255)
                                                                                ccTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCa                                                                      >  4:1286465/1‑149 (MQ=255)
                                                                                ccTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCa                                                                      >  3:6561944/1‑149 (MQ=255)
                                                                                     gcTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGATCCGGCATCCTc                                                                 <  3:533048/149‑1 (MQ=255)
                                                                                      cTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCa                                                                <  2:1526762/149‑1 (MQ=255)
                                                                                       tACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCAc                                                               >  2:99871/1‑149 (MQ=255)
                                                                                       tACGAAGATCGCACCCATCTCTACCCCATCGGACAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCAc                                                               >  3:4520316/1‑149 (MQ=255)
                                                                                          gAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGa                                                                                 >  2:3074021/1‑128 (MQ=255)
                                                                                          gAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCAt                                                            <  4:6889546/149‑1 (MQ=255)
                                                                                             gATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATAATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATgcg                                                         >  4:4695669/1‑149 (MQ=255)
                                                                                               tCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCt                                                       >  3:5137378/1‑149 (MQ=255)
                                                                                                  cACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCtttt                                                    <  4:4242338/149‑1 (MQ=255)
                                                                                                    cccATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTc                                                  <  2:6498850/149‑1 (MQ=255)
                                                                                                     ccATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCAc                                                               <  2:232278/135‑1 (MQ=255)
                                                                                                     ccATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCa                                                 >  1:2723051/1‑149 (MQ=255)
                                                                                                     ccATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCa                                                 >  4:1997896/1‑149 (MQ=255)
                                                                                                       aTCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAAt                                               >  2:1790786/1‑149 (MQ=255)
                                                                                                       aTCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGCGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAAt                                               <  2:3791949/149‑1 (MQ=255)
                                                                                                             aCCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAgc                                         <  1:5792442/149‑1 (MQ=255)
                                                                                                              ccccATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGc                                                                                            >  2:739635/1‑97 (MQ=255)
                                                                                                              ccccATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGc                                                                                   <  4:4465587/106‑1 (MQ=255)
                                                                                                               cccATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTcc                                                                                                                       >  1:3701489/1‑69 (MQ=255)
                                                                                                               cccATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAgcgc                                       >  3:4297339/1‑149 (MQ=255)
                                                                                                                 cATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCgg                                     >  3:3101928/1‑149 (MQ=255)
                                                                                                                    cGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCaa                                  >  2:6070051/1‑149 (MQ=255)
                                                                                                                     ggAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAAt                                 <  2:4660870/149‑1 (MQ=255)
                                                                                                                       agaACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGCTCAGCGACGGTTCCGGAATCCTCACCATGCGCTTTTTCAATTTCAGCGCCGCAATGa                               <  4:1441305/149‑1 (MQ=255)
                                                                                                                          actactGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCgg                                     >  3:5582868/1‑140 (MQ=255)
                                                                                                                          actactGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGaaaa                            <  4:1690754/149‑1 (MQ=255)
                                                                                                                           ctactGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGaaaaa                           <  1:1871171/149‑1 (MQ=255)
                                                                                                                            tactGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAAt                          <  4:535089/149‑1 (MQ=255)
                                                                                                                               tGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGc                       <  2:5168041/149‑1 (MQ=255)
                                                                                                                               tGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGc                       <  4:3612011/149‑1 (MQ=255)
                                                                                                                               tGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGc                       >  2:3322885/1‑149 (MQ=255)
                                                                                                                                  cGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTc                                                                           >  3:2053834/1‑94 (MQ=255)
                                                                                                                                    ggCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTggc                  <  1:2402409/149‑1 (MQ=255)
                                                                                                                                       gTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTggcggc               <  4:2855259/149‑1 (MQ=255)
                                                                                                                                        tttATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGgcggcg              <  4:6504378/149‑1 (MQ=255)
                                                                                                                                         ttATGCCACGGCGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGcggcgg             <  3:6091258/149‑1 (MQ=255)
                                                                                                                                         ttATGCCACGGCGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTACCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGcggcgg             <  4:5139516/149‑1 (MQ=255)
                                                                                                                                             gCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGgccg         >  2:5740724/1‑149 (MQ=255)
                                                                                                                                               cACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGgccgcc       >  3:4922108/1‑149 (MQ=255)
                                                                                                                                                aCGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGccgccg      <  1:1582155/149‑1 (MQ=255)
                                                                                                                                                 cGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCgt     <  1:99132/149‑1 (MQ=255)
                                                                                                                                                    tGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTa  >  2:4104800/1‑149 (MQ=255)
                                                                                                                                                    tGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGCccgccgccag  <  4:380187/149‑5 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGGGCGTTGGCGCAGCACTTAGTAACAAACTGGCGAAAATCAACCTGCATACCGTGCAGGATCTACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTA  >  NC_000913/3825252‑3825548

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