Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,272,546 A→T 6.8% T14S (ACG→TCG)  chaC → cation transport regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,272,5460AT 6.8% 238.6 / 5.3 83T14S (ACG→TCG) chaCcation transport regulator
Reads supporting (aligned to +/- strand):  ref base A (36/41);  new base T (3/3);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

GTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTT  >  NC_000913/1272398‑1272694
                                                                                                                                                    |                                                                                                                                                       
gTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTaaaa                                                                                                                                                         <  1:856247/149‑1 (MQ=255)
        aTTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGTTTTACAATCGGCAttc                                                                                                                                                 <  1:324708/149‑2 (MQ=255)
          tGCATATTTTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGTAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGTTTTACAATCGGCAttcat                                                                                                                                               <  2:959156/149‑4 (MQ=255)
          tGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTgg                                                                                                                                               >  2:460978/1‑149 (MQ=255)
           gCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGGAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCAttcatc                                                                                                                                              <  2:167332/149‑5 (MQ=255)
             atatTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGc                                                                                                                                            >  2:436723/1‑149 (MQ=255)
               atTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGGATGCCGATTGTAAAACGGCATTTGGTGCCa                                                                                                                                          <  1:1163044/149‑1 (MQ=255)
                ttGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGCAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCAt                                                                                                                                         <  2:1308717/149‑1 (MQ=255)
                   tCCCGTTAGTGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTgaa                                                                                                                                     >  1:1382924/1‑149 (MQ=255)
                    cccGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTgaa                                                                                                                                     <  1:635307/149‑1 (MQ=255)
                    cccGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTgaa                                                                                                                                     <  1:828145/149‑1 (MQ=255)
                     ccGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTgaag                                                                                                                                    >  1:256745/1‑149 (MQ=255)
                       gTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCGGGGATGTATGCAGTGGCGGAAGTGTATGGTGATATCGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTtaagaa                                                                                                                                  >  2:163072/1‑149 (MQ=255)
                       gTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTgaagaa                                                                                                                                  <  2:561190/149‑1 (MQ=255)
                       gTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTgaagaa                                                                                                                                  <  1:1325648/149‑1 (MQ=255)
                           gTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAATACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCAc                                                                                                                              <  1:850446/149‑1 (MQ=255)
                            tGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCAct                                                                                                                             >  2:541307/1‑149 (MQ=255)
                            tGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCAct                                                                                                                             >  2:420034/1‑149 (MQ=255)
                               tttCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCtt                                                                                                                          >  2:1104645/1‑149 (MQ=255)
                                ttCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTa                                                                                                                         <  2:1473930/149‑1 (MQ=255)
                                ttCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTa                                                                                                                         >  1:1136547/1‑149 (MQ=255)
                                 tCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTAt                                                                                                                        >  2:1145421/1‑149 (MQ=255)
                                 tCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTAt                                                                                                                        <  1:586309/149‑1 (MQ=255)
                                   aaaaTGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATg                                                                                                                       >  2:1296724/1‑148 (MQ=255)
                                   aaaaTGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATg                                                                                                                       <  1:1296724/148‑1 (MQ=255)
                                    aaaTGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGt                                                                                                                     <  2:1560191/149‑1 (MQ=255)
                                       ttaGCATTAAAAATGTACGGAAGACACCTAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTcag                                                                                                                  <  2:687723/147‑1 (MQ=255)
                                       tGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTcag                                                                                                                  <  2:1564037/149‑1 (MQ=255)
                                             gTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGaac                                                                                                            >  2:1417415/1‑149 (MQ=255)
                                                aaaaTGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGaacaac                                                                                                         >  2:327420/1‑149 (MQ=255)
                                                      tCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGggcgg                                                                                                   <  2:750407/149‑1 (MQ=255)
                                                       ccGGAAGACACCAAAAAGTTGTCGCAGGTAAGTATGCAGTGGCTGAAGTGTAAGTTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGggcggc                                                                                                  <  1:1451224/149‑1 (MQ=255)
                                                       ccGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGAAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGggcggc                                                                                                  <  2:669769/149‑1 (MQ=255)
                                                        cGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCt                                                                                                 >  2:854887/1‑149 (MQ=255)
                                                          gAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGGGGCTTc                                                                                               >  2:1495021/1‑149 (MQ=255)
                                                             gACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCt                                                                                            >  2:904342/1‑149 (MQ=255)
                                                              acacCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑g                                                                                        >  1:1265564/1‑149 (MQ=255)
                                                                  cAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑ggcgg                                                                                    >  1:1368648/1‑149 (MQ=255)
                                                                  cAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑ggcgg                                                                                    >  1:1400074/1‑149 (MQ=255)
                                                                        gTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGCCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGc                                                                              >  2:863776/1‑149 (MQ=255)
                                                                        gTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTAGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGc                                                                              >  2:1375631/1‑149 (MQ=255)
                                                                         ttGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCt                                                                             >  2:1210497/1‑149 (MQ=255)
                                                                           gTCGCAGGGAAGTATGCAGGGGCGGAAGGGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTgg                                                                           <  2:147465/149‑1 (MQ=255)
                                                                            tCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGCAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGc                                                                          <  1:1535861/149‑1 (MQ=255)
                                                                            tCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGc                                                                          <  1:904342/149‑1 (MQ=255)
                                                                                aGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGt                                                                      <  1:541307/149‑1 (MQ=255)
                                                                                aGGGAAGTATGCACGGGCGGAAGTGAAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGTGGCCATTGAAGAATCACACTTATGGTCCGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGt                                                                      <  2:768494/149‑1 (MQ=255)
                                                                                 gggAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTc                                                                     <  1:1539932/149‑1 (MQ=255)
                                                                                 gggAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTc                                                                     <  2:1282392/149‑1 (MQ=255)
                                                                                 gggAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTc                                                                     <  1:331109/149‑1 (MQ=255)
                                                                                 gggAAGTATGCAGTGGCGGAAGTGTAAGGGGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTc                                                                     <  1:1457574/149‑1 (MQ=255)
                                                                                 cggAAGTATGCAGTGGCGGAAGCGTACGGTGCTAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTc                                                                     <  1:209645/148‑1 (MQ=255)
                                                                                  ggAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCg                                                                    >  2:661911/1‑149 (MQ=255)
                                                                                     aGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGAcc                                                                 <  2:256745/149‑1 (MQ=255)
                                                                                         tGCAGTGGCGGAAGTGTAAGGTGATAACTCGTGATTTATTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGAGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCtgat                                                             <  1:873406/149‑1 (MQ=255)
                                                                                           cAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCtgatga                                                           >  2:1188975/1‑149 (MQ=255)
                                                                                            aGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAg                                                          >  1:653754/1‑149 (MQ=255)
                                                                                            aGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAg                                                          >  1:498912/1‑149 (MQ=255)
                                                                                              tGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCTGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAggg                                                        >  2:833637/1‑149 (MQ=255)
                                                                                                    aaGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGCCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGgt                                                  <  1:863776/149‑1 (MQ=255)
                                                                                                      gtgtAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGgtgt                                                <  2:1151574/149‑1 (MQ=255)
                                                                                                          acggTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAATCGGCATTTGGTGCCATTGAAGACTCACTCTTATGGTCAGCAGAACAACGGGCGGCGTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGAt                                            <  2:1222524/147‑1 (MQ=255)
                                                                                                           aGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATc                                           >  1:1010934/1‑149 (MQ=255)
                                                                                                            ggTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCt                                          >  2:1502568/1‑149 (MQ=255)
                                                                                                              tGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTc                                        <  1:1293885/149‑1 (MQ=255)
                                                                                                              tGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTc                                        >  1:122335/1‑149 (MQ=255)
                                                                                                                     gcgTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATg                                 >  1:1574406/1‑149 (MQ=255)
                                                                                                                      cgTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTAGATCTACGTCTaagc                                >  2:513349/1‑146 (MQ=255)
                                                                                                                        tGATTTATTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGACCACCGGGCGGCTTCGCT‑‑‑GGCGGCGCCGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGGTATGGAt                              >  2:151023/1‑149 (MQ=255)
                                                                                                                           tttCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGt                           <  1:420034/149‑1 (MQ=255)
                                                                                                                            ttCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGtt                          <  2:1265564/149‑1 (MQ=255)
                                                                                                                            ttCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGtt                          <  2:327819/149‑1 (MQ=255)
                                                                                                                               ttgttgATTGCCGATTGTAAAACGGCTTTTGGTGCCATTGAAGAATCACTCTTATGGTCAGGAGAAGAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAt                       <  1:608214/149‑1 (MQ=255)
                                                                                                                                caccaaaTGCCGTTTTACAATCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAtg                      >  1:1566646/6‑149 (MQ=255)
                                                                                                                                 gatgaATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAtgt                     <  1:1417415/149‑1 (MQ=255)
                                                                                                                                 gatgaATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAtgt                     <  1:327420/149‑1 (MQ=255)
                                                                                                                                 accaaaTGCCGTTTTACAATCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAtgt                     >  2:345613/5‑149 (MQ=255)
                                                                                                                                   caaaTGCCGTTTTACAATCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTTTGGATCGTTGATGTgg                   >  1:1138026/3‑149 (MQ=255)
                                                                                                                                         ccGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCg             <  1:661911/149‑1 (MQ=255)
                                                                                                                                           gATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGc           >  2:682419/1‑149 (MQ=255)
                                                                                                                                             ttGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGCATCGTTGATGTTGAATCCGGCAc         >  2:60645/1‑149 (MQ=255)
                                                                                                                                             ttGTAAAACGGCATTTGGAGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCAc         <  1:1188975/149‑1 (MQ=255)
                                                                                                                                                tAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTgg      <  2:1382924/149‑1 (MQ=255)
                                                                                                                                                tAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGACACTgg      >  1:109129/1‑149 (MQ=255)
                                                                                                                                                 aaaaCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGa     >  1:374426/1‑149 (MQ=255)
                                                                                                                                                  aaaCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTg       >  2:593974/1‑149 (MQ=255)
                                                                                                                                                  aaaCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAg    <  1:1104645/149‑1 (MQ=255)
                                                                                                                                                    aCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGtt  <  2:498912/149‑1 (MQ=255)
                                                                                                                                                    aCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCAGGCACTGGAGtt  >  2:621614/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
GTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCT‑‑‑GGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTT  >  NC_000913/1272398‑1272694

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: