Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,272,541 G→T 6.1% C12F (TGT→TTT) ‡ chaC → cation transport regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,272,5410GT 6.1% 251.7 / 6.0 82C12F (TGT→TTT) ‡chaCcation transport regulator
Reads supporting (aligned to +/- strand):  ref base G (35/42);  new base T (3/2);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

GTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCAC  >  NC_000913/1272398‑1272687
                                                                                                                                               |                                                                                                                                                  
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                                                                                                                        <  1:1325648/149‑1 (MQ=255)
                       gTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTgaagaa                                                                                                                        <  2:561190/149‑1 (MQ=255)
                           gTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAATACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCAc                                                                                                                    <  1:850446/149‑1 (MQ=255)
                            tGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCAct                                                                                                                   >  2:420034/1‑149 (MQ=255)
                            tGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCAct                                                                                                                   >  2:541307/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                                                                                                              <  1:586309/149‑1 (MQ=255)
                                 tCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTAt                                                                                                              >  2:1145421/1‑149 (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)
                                                              acacCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTg                                                                                 >  1:1265564/1‑149 (MQ=255)
                                                                  cAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTggcgg                                                                             >  1:1368648/1‑149 (MQ=255)
                                                                  cAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTggcgg                                                                             >  1:1400074/1‑149 (MQ=255)
                                                                        gTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGCCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGc                                                                       >  2:863776/1‑149 (MQ=255)
                                                                        gTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTAGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGc                                                                       >  2:1375631/1‑149 (MQ=255)
                                                                         ttGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCt                                                                      >  2:1210497/1‑149 (MQ=255)
                                                                           gTCGCAGGGAAGTATGCAGGGGCGGAAGGGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTgg                                                                    <  2:147465/149‑1 (MQ=255)
                                                                            tCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGCAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGc                                                                   <  1:1535861/149‑1 (MQ=255)
                                                                            tCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGc                                                                   <  1:904342/149‑1 (MQ=255)
                                                                                aGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGt                                                               <  1:541307/149‑1 (MQ=255)
                                                                                aGGGAAGTATGCACGGGCGGAAGTGAAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGTGGCCATTGAAGAATCACACTTATGGTCCGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGt                                                               <  2:768494/149‑1 (MQ=255)
                                                                                 gggAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTc                                                              <  1:1539932/149‑1 (MQ=255)
                                                                                 gggAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTc                                                              <  2:1282392/149‑1 (MQ=255)
                                                                                 gggAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTc                                                              <  1:331109/149‑1 (MQ=255)
                                                                                 gggAAGTATGCAGTGGCGGAAGTGTAAGGGGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTc                                                              <  1:1457574/149‑1 (MQ=255)
                                                                                 cggAAGTATGCAGTGGCGGAAGCGTACGGTGCTAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTc                                                              <  1:209645/148‑1 (MQ=255)
                                                                                  ggAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCg                                                             >  2:661911/1‑149 (MQ=255)
                                                                                     aGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGAcc                                                          <  2:256745/149‑1 (MQ=255)
                                                                                         tGCAGTGGCGGAAGTGTAAGGTGATAACTCGTGATTTATTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGAGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCtgat                                                      <  1:873406/149‑1 (MQ=255)
                                                                                           cAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCtgatga                                                    >  2:1188975/1‑149 (MQ=255)
                                                                                            aGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAg                                                   >  1:653754/1‑149 (MQ=255)
                                                                                            aGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAg                                                   >  1:498912/1‑149 (MQ=255)
                                                                                              tGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCTGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAggg                                                 >  2:833637/1‑149 (MQ=255)
                                                                                                    aaGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGCCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGgt                                           <  1:863776/149‑1 (MQ=255)
                                                                                                      gtgtAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGgtgt                                         <  2:1151574/149‑1 (MQ=255)
                                                                                                          acggTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAATCGGCATTTGGTGCCATTGAAGACTCACTCTTATGGTCAGCAGAACAACGGGCGGCGTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGAt                                     <  2:1222524/147‑1 (MQ=255)
                                                                                                           aGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATc                                    >  1:1010934/1‑149 (MQ=255)
                                                                                                            ggTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCt                                   >  2:1502568/1‑149 (MQ=255)
                                                                                                              tGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTc                                 <  1:1293885/149‑1 (MQ=255)
                                                                                                              tGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTc                                 >  1:122335/1‑149 (MQ=255)
                                                                                                                     gcgTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATg                          >  1:1574406/1‑149 (MQ=255)
                                                                                                                      cgTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTAGATCTACGTCTaagc                         >  2:513349/1‑146 (MQ=255)
                                                                                                                        tGATTTATTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGACCACCGGGCGGCTTCGCTGGCGGCGCCGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGGTATGGAt                       >  2:151023/1‑149 (MQ=255)
                                                                                                                           tttCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGt                    <  1:420034/149‑1 (MQ=255)
                                                                                                                            ttCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGtt                   <  2:327819/149‑1 (MQ=255)
                                                                                                                            ttCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGtt                   <  2:1265564/149‑1 (MQ=255)
                                                                                                                               ttgttgATTGCCGATTGTAAAACGGCTTTTGGTGCCATTGAAGAATCACTCTTATGGTCAGGAGAAGAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAt                <  1:608214/149‑1 (MQ=255)
                                                                                                                                caccaaaTGCCGTTTTACAATCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAtg               >  1:1566646/6‑149 (MQ=255)
                                                                                                                                 gatgaATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAtgt              <  1:1417415/149‑1 (MQ=255)
                                                                                                                                 gatgaATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAtgt              <  1:327420/149‑1 (MQ=255)
                                                                                                                                 accaaaTGCCGTTTTACAATCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAtgt              >  2:345613/5‑149 (MQ=255)
                                                                                                                                   caaaTGCCGTTTTACAATCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTTTGGATCGTTGATGTgg            >  1:1138026/3‑149 (MQ=255)
                                                                                                                                         ccGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCg      <  1:661911/149‑1 (MQ=255)
                                                                                                                                           gATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGc    >  2:682419/1‑149 (MQ=255)
                                                                                                                                             ttGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGCATCGTTGATGTTGAATCCGGCAc  >  2:60645/1‑149 (MQ=255)
                                                                                                                                             ttGTAAAACGGCATTTGGAGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCAc  <  1:1188975/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCAC  >  NC_000913/1272398‑1272687

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: