Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,835,243 T→C 5.2% Q229R (CAG→CGG)  yicJ ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,835,2430TC 5.2% 234.7 / 2.6 85Q229R (CAG→CGG) yicJputative transporter
Reads supporting (aligned to +/- strand):  ref base T (42/36);  new base C (3/2);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-01

ATCAGGTTACCCACGCAATAAGTGGTGAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGAATGCCACCACGGAAAGGACCGCGATACCGCCCTGGAAACCGA  >  NC_000913/3835096‑3835390
                                                                                                                                                   |                                                                                                                                                   
aTCAGGCTACCCACGCAATAAGTGGTGAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTg                                                                                                                                                    <  2:499424/149‑1 (MQ=255)
 tCAGGTTACCCACGCAATAAGTGGTGAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTgg                                                                                                                                                   >  1:340560/1‑149 (MQ=255)
     gTTACCCACGCAATAAGTGGTGAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGt                                                                                                                                               >  2:597364/1‑149 (MQ=255)
       tACTCACGCAATAAGTGGTGAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGAGATTGAAAATGGTTAGTATACCGACAATCCGCCACTGGTCGttt                                                                                                                                             <  1:579572/149‑1 (MQ=255)
        aCCCACGCAATAAGTTGTGAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGtttt                                                                                                                                            <  2:384116/149‑1 (MQ=255)
        aCCCACGCAATAAGTGGTGAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGtttt                                                                                                                                            >  1:139909/1‑149 (MQ=255)
        aCCCACGCAATAAGTGGTGAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGtttt                                                                                                                                            <  2:900411/149‑1 (MQ=255)
         cccACGCAATAAGTGGTGAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTg                                                                                                                                           >  2:950100/1‑149 (MQ=255)
          ccACGCAATAAGTGGTGAGAAAAGCGACAAACACTTCCGGCGTGCCCAGCATCGATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGc                                                                                                                                          >  1:512726/1‑149 (MQ=255)
          ccACGCAATAAGTGGTGAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCCCACCGCTGGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGc                                                                                                                                          >  1:32896/1‑149 (MQ=255)
            aggCAATAAGTGGTGAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCCCCGCGGACGCACACCGCAAGGATGTGAAAAATGGTTAGTAAACCGACAAACCACCAATGGGCGGATTGcaa                                                                                                                                        >  1:270714/3‑147 (MQ=255)
            aCGCAATAAGTGGTGAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTTACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCCGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCa                                                                                                                                        >  1:512494/1‑149 (MQ=255)
            aCGCAATAAGTGGTGAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACACTCCGCCACTGCTCGTTTTGCCa                                                                                                                                        >  2:170610/1‑149 (MQ=255)
                aaTAAGTGGTGAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTCCCCGata                                                                                                                                    >  1:406171/1‑149 (MQ=255)
                        gtgAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTaa                                                                                                                            <  1:359834/149‑1 (MQ=255)
                        gtgAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACCGGTCGTTTTGCCAGATATCACg                                                                                                                               >  2:683687/1‑146 (MQ=255)
                        gtgAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACCGGTCGTTTTGCCAGATATCACg                                                                                                                               <  1:683698/146‑1 (MQ=255)
                        gtgAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACCGGTCGTTTTGCCAGATATCACg                                                                                                                               <  1:683687/146‑1 (MQ=255)
                        gtgAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACCGGTCGTTTTGCCAGATATCACg                                                                                                                               >  2:683698/1‑146 (MQ=255)
                              aaaGCCACAAACACTTCCCGGGTGCCCAAAATCCATGTGACGTAATACATCAGCGCCCCACCGCGTACGCAAACCGCTAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTc                                                                                                                      <  2:561631/149‑1 (MQ=255)
                                aGCGACAAACCCTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATAGTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTccc                                                                                                                    <  2:1133844/149‑1 (MQ=255)
                                      aaaCACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGCTCGTTTTGCCAGATATCACGTAAATCGTCCCGCACAg                                                                                                              >  2:329546/1‑149 (MQ=255)
                                      aaaCACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCCCCACCGGTCGTTTTGCCAGAAATCACGTAAATCTTCCCGCATAg                                                                                                              >  2:25525/1‑149 (MQ=255)
                                       aaCACTTCCGGCGTGCCCAAAATCCTTGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGa                                                                                                             >  1:579386/1‑149 (MQ=255)
                                           cTTCCGGCGAGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACgtt                                                                                                         <  2:874825/149‑1 (MQ=255)
                                            ttCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATAACACGTAAATCTTCCCGCATAGACgttg                                                                                                        >  2:153924/1‑149 (MQ=255)
                                              ccGGCTTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCTCCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACgttgtt                                                                                                      <  1:90440/149‑1 (MQ=255)
                                                  cGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAg                                                                                                  <  1:286755/149‑1 (MQ=255)
                                                     gCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAggtg                                                                                               <  2:689080/149‑1 (MQ=255)
                                                      cccAAAATCCATGTGACGTAATACATCATCGCCCCGCCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTGCCCGCATAGACGTTGTTGTAggtgg                                                                                              <  1:643090/149‑1 (MQ=255)
                                                       ccAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGGCGTTGTTGTAggtggt                                                                                             >  1:609121/1‑149 (MQ=255)
                                                       ccAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGGAggtggt                                                                                             >  1:857084/1‑149 (MQ=255)
                                                        cAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGgtggtg                                                                                            >  1:881644/1‑149 (MQ=255)
                                                        cAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGgtggtg                                                                                            >  2:70076/1‑149 (MQ=255)
                                                        cAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGAAGTTGTTGgtggtggtg                                                                                            >  1:719770/1‑149 (MQ=255)
                                                         aaaaTCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACCCCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGc                                                                                           >  1:1089771/1‑149 (MQ=255)
                                                         aaaaTCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGc                                                                                           >  2:888576/1‑149 (MQ=255)
                                                         aaaaTCCAGGTGACGCAATCCCTCATCGCCCCACCGCGTACCCCCGCCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGc                                                                                           <  1:1117456/149‑1 (MQ=255)
                                                             tCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCa                                                                                       >  1:28712/1‑149 (MQ=255)
                                                               cttgtgACGTAATACATCATCGCTCCCCCGCGTACGCACACCGCCAGGAGATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAAc                                                                                     <  1:457982/147‑1 (MQ=255)
                                                               aaTGTTACGTAATAAATCATCGCCCCACCGCGTTAGCACACCGCCGGCATATTCAAAATGGTTAGTAAACCGACAATCCGGCACTGGTCGGTTTGCCAGATATCACGTAAAGCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAAc                                                                                     <  1:322125/148‑1 (MQ=255)
                                                                aTGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGCTATCCCGTAAATCTTCCCGCGTAGACGTAGTTGTAGTTGACACTTCAAcg                                                                                    >  1:80588/1‑149 (MQ=255)
                                                                    gACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGtt                                                                                >  1:609246/1‑149 (MQ=255)
                                                                    gACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGtt                                                                                >  1:769339/1‑149 (MQ=255)
                                                                     aCGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCACCCCGTTc                                                                               >  2:543128/1‑149 (MQ=255)
                                                                       gTAATACGTCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCtt                                                                             <  2:150692/149‑1 (MQ=255)
                                                                             catcatCGCCCCACCGCGTACGCACACCGCCAGGATATTCAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAgtgg                                                                       <  2:86406/149‑1 (MQ=255)
                                                                                 atcGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTgaa                                                                   <  2:153463/149‑1 (MQ=255)
                                                                                   cGCCACACCGCGTACGAACACCGACAGGATATTGAAAATGGTTAGTAAACCGACAAGCCGCCACTGGTCGTTTTGCCAGAGATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTgaaga                                                                 <  1:254110/149‑1 (MQ=255)
                                                                                     ccccACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGCCAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCAGCATAGCCGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGaaa                                                               >  1:286995/1‑149 (MQ=255)
                                                                                     ccccACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAGGaaa                                                               >  2:87841/1‑149 (MQ=255)
                                                                                       cccgcgcgTACGCACACCGCTAGGATATTGAAAATGGTTCGTAAACCGACAATCGTCCACTGGGCTTTTTGCCTGATATCAAGTAGATCTTCACGCCTAGACGTGGTTGTAGTCGGTGCTTCAACGCGTTCTTTGGAGGTGGAGGAACa                                                             <  2:555035/145‑1 (MQ=255)
                                                                                        cACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGGGAAGAAACaa                                                            >  2:457073/1‑149 (MQ=255)
                                                                                          ccGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACaaaa                                                          <  1:597364/149‑1 (MQ=255)
                                                                                            gctgactcacaCCGCCAGGATATTGAAAATGGTTAGTAAACCGACAACCCACCTCTGGTCGCTTTGCCAGATATCACTTAAATCTTCCCGCATAGTCGTTGTTGTAGGTAGTGCTTCAACGCGTTCTTTACTGGTGAAGAAACAAAATg                                                        <  2:270714/142‑1 (MQ=255)
                                                                                            gcgTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCAAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGATCTTTAGTGGGGACGAAACaaaaag                                                        >  1:706563/1‑147 (MQ=255)
                                                                                             cggaCGCACACCGCCAGGAGATTGAAAATGGTTAGTAACCCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGc                                                       <  1:19133/146‑1 (MQ=255)
                                                                                             cgTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGc                                                       <  1:525033/149‑1 (MQ=255)
                                                                                                   cacaCCGCCAGGATATTGAAAATGGTTAGCAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGcat                                                 <  2:139909/149‑1 (MQ=255)
                                                                                                    acacCGCCAGGATATTGAAAATGGTTAGTAAACCGACCATCCGCCACTGGTCGATTTGGCAGATATCGCGCACATCTTCCAACATAGACGTTGTCGTAGGTGGTGCTTCAACGCGTTCTCTAGTGGGGAAGAAACAAAATGCCAGcatc                                                >  2:478087/1‑149 (MQ=255)
                                                                                                         gCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGaa                                           <  2:28712/149‑1 (MQ=255)
                                                                                                           cAGGATATTTAAAATGGAAAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTACTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATACCAGCATCATGAAta                                         >  2:1114286/1‑148 (MQ=255)
                                                                                                              gagagtGAAAATGGTTACTAAACCGACAATCCGCCACGGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGGCGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGAATGcca                                      <  2:286995/144‑1 (MQ=255)
                                                                                                              gATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGTCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCCTCAACGCGTTCTTTAGTGGGGAAGAAACAAAATGCCAGCATCAGGAATGccc                                      >  2:953603/1‑148 (MQ=255)
                                                                                                              gATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGAATGcca                                      <  1:950100/149‑1 (MQ=255)
                                                                                                                   tGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGAATGCCACCACg                                 >  2:497263/1‑149 (MQ=255)
                                                                                                                     aaaaTGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGAATGCCACCACGGa                               >  1:877156/1‑149 (MQ=255)
                                                                                                                      aaaTGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCAAGCATCATGAATGCCACCACGGaa                              <  1:856009/149‑1 (MQ=255)
                                                                                                                       aaTGGTTGGTAAACCGACAATCCGCCTCTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGAATGCCACCACGGaaa                             <  1:335556/149‑1 (MQ=255)
                                                                                                                          ggTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAAAGCCAGCATCATGAATGCCACCACGGAAAGGa                          >  1:69398/1‑149 (MQ=255)
                                                                                                                            ttAGTAAACCGACAATCCGCCACTGGTTGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGAATGCCACCACGGAAAGGAcc                        >  2:459167/1‑149 (MQ=255)
                                                                                                                              aGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTGCTTTAGTGCTGACGAATCAAAATGCCAGCATCATGAATGCCACTACGTAAAGGACcgc                      >  2:1036560/1‑149 (MQ=255)
                                                                                                                                ttaaCCCACACTCCGGCACTGGTTGATTTATCAGATAGCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGTCAGCATCTTGAATGGCACCACGCAAAGGACCGCGa                    <  2:525631/147‑1 (MQ=255)
                                                                                                                                tAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGAATGCCACCACGGAAAGGACCGCGa                    <  2:579386/149‑1 (MQ=255)
                                                                                                                                tAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGAATGCCACCACGGAAAGGACCGCGa                    >  1:778062/1‑149 (MQ=255)
                                                                                                                                tAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGAATGCCACCACGGAAAGGAACGCGa                    >  2:660527/1‑149 (MQ=255)
                                                                                                                                tAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATAATGAATGCCACCACGGAAAGGACCGCGa                    >  2:298305/1‑149 (MQ=255)
                                                                                                                                 aaaCCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGGTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGCAATCCACCACGGAAAGGACCGGGAt                   >  2:231688/1‑149 (MQ=255)
                                                                                                                                   aCCGACAGTCCGCTCTGGGTCCTTTTACCAGATATCACGTAAATCTTCCCGCTTAGACGTTGTTGTAGGTGGTGCTGCAACGCGTTCTTTAGTTGTGAAGAAACAAAATGCCCGCATTATTTTCGCCACCACGGAAAGGACCGCGATAc                 <  1:280016/149‑1 (MQ=255)
                                                                                                                                      gACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGAATGCCACCACGGAAAGGACCGCGATACCGc              <  2:769339/149‑1 (MQ=255)
                                                                                                                                         aaTCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGAATGCCACCACGGAAAGGACCGCGATACCGCCCt           <  2:609246/149‑1 (MQ=255)
                                                                                                                                              gCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGAATGCCACCACGGAAAGGACCGCGATACCGCCCTGGaaa      >  2:395878/1‑149 (MQ=255)
                                                                                                                                              gCCACTGGTCGTTGTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGAATGCCACCACGGAAAGGACCGCGATACCGCCCTGGaaa      <  1:170610/149‑1 (MQ=255)
                                                                                                                                              gCCACTCGGCGATTTGCCAGTTATCAGGTATATCTTCCCGCAGAGACGTGGTTGTAGGGGGTGCTTCAACGCGTCCTTTAGTGGTGAAGAAACAAAATGCCAGCAGCATGAATGCCACCACGGAAAGGACCGCGATACCGCCCTGGaaa      <  2:32896/149‑1 (MQ=255)
                                                                                                                                                  cTGGTCGTTTTGCCAGATATCATGTACTTGTTCCCGCATAGACGTTGTTGCAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGAATGCCACCACGGAAAGGACCGCGATACCGCCCTGGAAACCGa  <  1:153924/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATCAGGTTACCCACGCAATAAGTGGTGAGAAAAGCGACAAACACTTCCGGCGTGCCCAAAATCCATGTGACGTAATACATCATCGCCCCACCGCGTACGCACACCGCCAGGATATTGAAAATGGTTAGTAAACCGACAATCCGCCACTGGTCGTTTTGCCAGATATCACGTAAATCTTCCCGCATAGACGTTGTTGTAGGTGGTGCTTCAACGCGTTCTTTAGTGGTGAAGAAACAAAATGCCAGCATCATGAATGCCACCACGGAAAGGACCGCGATACCGCCCTGGAAACCGA  >  NC_000913/3835096‑3835390

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