Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,115,423 C→T C23Y (TGT→TAT)  mdtG ← putative drug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,115,4230CT98.2% 419.9 / ‑1.6 109C23Y (TGT→TAT) mdtGputative drug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (66/41);  total (67/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

Only 100 of 109 total aligned reads displayed.

GTGAGGCGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAG  >  NC_000913/1115277‑1115570
                                                                                                                                                  |                                                                                                                                                   
gTGAGGCGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGcc                                                                                                                                                 >  1:660105/1‑151 (MQ=255)
  gAGGCGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCt                                                                                                                                                >  2:564174/1‑150 (MQ=255)
  gAGGCGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCt                                                                                                                                                >  1:1191457/1‑150 (MQ=255)
   aGGCGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCt                                                                                                                                                >  2:623822/1‑149 (MQ=255)
       gATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGcc                                                                                                                                            >  2:253714/1‑149 (MQ=255)
        aTGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAgg                                                                                                                                         <  1:1634588/151‑1 (MQ=255)
         tGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGc                                                                                                                                        >  1:958177/1‑151 (MQ=255)
             cGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGcgacg                                                                                                                                    >  2:1507697/1‑151 (MQ=255)
                 aaTAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATc                                                                                                                                 >  2:629471/1‑150 (MQ=255)
                  aTAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAg                                                                                                                               <  2:344049/151‑1 (MQ=255)
                    aaaaaTGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGgt                                                                                                                             <  1:499713/151‑1 (MQ=255)
                      aaaTGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGgttt                                                                                                                           >  1:1538891/1‑151 (MQ=255)
                       aaTGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGgtttc                                                                                                                          >  2:1043071/1‑151 (MQ=255)
                        aTGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTa                                                                                                                                               >  1:425663/1‑129 (MQ=255)
                        aTGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTa                                                                                                                                               <  2:425698/129‑1 (MQ=255)
                        aTGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGgtttcg                                                                                                                         <  1:1218671/151‑1 (MQ=255)
                         tGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGgtttcgt                                                                                                                        >  2:1378050/1‑151 (MQ=255)
                         tGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGgtttcgt                                                                                                                        >  1:948445/1‑151 (MQ=255)
                          gTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGgtttcgt                                                                                                                        <  1:1579368/150‑1 (MQ=255)
                            aaTGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGgtttcgtttc                                                                                                                     >  1:1401682/1‑151 (MQ=255)
                              tGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTcc                                                                                                                    >  1:1608610/1‑150 (MQ=255)
                               gCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTcc                                                                                                                    >  2:656623/1‑149 (MQ=255)
                                cTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGt                                                                                                                 >  2:1621263/1‑151 (MQ=255)
                                 tGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGtt                                                                                                                >  2:1480159/1‑151 (MQ=255)
                                  gAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGtt                                                                                                                <  1:1293321/150‑1 (MQ=255)
                                  gAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGtt                                                                                                                <  1:868798/150‑1 (MQ=255)
                                    aGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCCTTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTtat                                                                                                             <  1:1250705/151‑1 (MQ=255)
                                     gACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTtata                                                                                                            <  2:345396/151‑1 (MQ=255)
                                      aCAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTtat                                                                                                             >  2:1153531/1‑149 (MQ=255)
                                      aCAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTtat                                                                                                             >  1:697322/1‑149 (MQ=255)
                                        aaTACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAggg                                                                                                         >  1:354189/1‑151 (MQ=255)
                                         aTACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAgggg                                                                                                        >  2:722283/1‑151 (MQ=255)
                                           aCCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGgtg                                                                                                      >  2:1243459/1‑151 (MQ=255)
                                           aCCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGgtg                                                                                                      >  1:313068/1‑151 (MQ=255)
                                            ccGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGgtg                                                                                                      <  2:1228029/150‑1 (MQ=255)
                                              ggACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTc                                                                                                    >  1:1035677/1‑150 (MQ=255)
                                              ggACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTc                                                                                                    >  2:1622341/1‑150 (MQ=255)
                                               gACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCAt                                                                                                  >  1:1449676/1‑151 (MQ=255)
                                                   acaTATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTc                                                                                              <  2:660182/151‑1 (MQ=255)
                                                     atatTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTcac                                                                                            >  1:600969/1‑151 (MQ=255)
                                                     atatTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTcac                                                                                            >  1:1200826/1‑151 (MQ=255)
                                                     atatTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTcac                                                                                            >  1:1315517/1‑151 (MQ=255)
                                                      tatTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTcac                                                                                            >  1:648157/1‑150 (MQ=255)
                                                          cAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGt                                                                                       >  2:1533654/1‑151 (MQ=255)
                                                              gCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGAc                                                                                    >  2:31483/1‑150 (MQ=255)
                                                               cGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGGGACa                                                                                   >  2:1348409/1‑150 (MQ=255)
                                                                 gAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACAt                                                                                  >  2:271931/1‑149 (MQ=255)
                                                                     gACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAAt                                                                            >  2:1519582/1‑151 (MQ=255)
                                                                       ccGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATcc                                                                          >  2:238112/1‑151 (MQ=255)
                                                                        cGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCg                                                                         <  2:60792/151‑1 (MQ=255)
                                                                         tgTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGc                                                                        >  2:363814/2‑151 (MQ=255)
                                                                         ggTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGc                                                                        >  1:1693937/1‑151 (MQ=255)
                                                                          gTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCt                                                                       <  2:1292985/151‑1 (MQ=255)
                                                                                ccAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGt                                                                 <  2:1060301/151‑1 (MQ=255)
                                                                                      tgctCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGcc                                                            >  1:902244/1‑150 (MQ=255)
                                                                                          cAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTcaca                                                                                           >  1:736078/1‑115 (MQ=255)
                                                                                          cAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTcaca                                                                                           <  2:736136/115‑1 (MQ=255)
                                                                                          cAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAgg                                                         <  2:1191547/149‑1 (MQ=255)
                                                                                          cAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGgcg                                                       <  1:144543/151‑1 (MQ=255)
                                                                                           aaCGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGgcgc                                                      >  1:109229/1‑151 (MQ=255)
                                                                                           aaCGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGgcg                                                       <  1:545732/150‑1 (MQ=255)
                                                                                               tAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGt                                                   >  1:961313/1‑150 (MQ=255)
                                                                                                  agGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAAc                                               >  2:690834/1‑151 (MQ=255)
                                                                                                   ggggTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACa                                              <  1:185675/151‑1 (MQ=255)
                                                                                                   ggggTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACa                                              <  1:1172119/151‑1 (MQ=255)
                                                                                                   ggggTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACa                                              <  1:927429/151‑1 (MQ=255)
                                                                                                       tAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATa                                            >  1:170877/1‑149 (MQ=255)
                                                                                                        aagaagGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATaaga                                         <  2:1672076/151‑1 (MQ=255)
                                                                                                         agaagGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATaaga                                         >  1:1257197/1‑150 (MQ=255)
                                                                                                             gggCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGaaa                                    <  2:1260148/151‑1 (MQ=255)
                                                                                                             gggCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGaaa                                    <  1:271907/151‑1 (MQ=255)
                                                                                                               gCATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGaaa                                    >  2:447372/1‑149 (MQ=255)
                                                                                                                cATTACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGaaa                                    >  1:34094/1‑148 (MQ=255)
                                                                                                                  ttACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGAt                                >  1:679112/1‑150 (MQ=255)
                                                                                                                   tACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATa                               >  1:1531197/1‑150 (MQ=255)
                                                                                                                   tACCAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATa                               <  1:1136240/150‑1 (MQ=255)
                                                                                                                     ccAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAg                             <  2:236138/150‑1 (MQ=255)
                                                                                                                     ccAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGc                            >  2:1682090/1‑151 (MQ=255)
                                                                                                                     ccAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGc                            >  1:721307/1‑151 (MQ=255)
                                                                                                                      cAGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGca                           <  1:623766/151‑1 (MQ=255)
                                                                                                                       aGACTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGcac                          <  1:72496/151‑1 (MQ=255)
                                                                                                                          cTGAAGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCacac                        >  2:433821/1‑150 (MQ=255)
                                                                                                                             aaGGCGGCACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAAt                    >  1:1361501/1‑151 (MQ=255)
                                                                                                                                gcggcACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTAt                 <  1:1042992/151‑1 (MQ=255)
                                                                                                                                 cggcACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATg                >  1:849379/1‑151 (MQ=255)
                                                                                                                                  ggcACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATgc               <  1:1301761/151‑1 (MQ=255)
                                                                                                                                  ggcACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATgc               >  2:1524170/1‑151 (MQ=255)
                                                                                                                                    cACCGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATgcgc             >  1:1135187/1‑151 (MQ=255)
                                                                                                                                      ccGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGccc           <  1:780533/151‑1 (MQ=255)
                                                                                                                                       cGGTAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCg          >  2:784199/1‑151 (MQ=255)
                                                                                                                                         gTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGAc        <  1:253691/151‑1 (MQ=255)
                                                                                                                                          tAAGAAAATAGCCTGGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACt       >  1:1377376/1‑151 (MQ=255)
                                                                                                                                          tAAGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACt       >  2:166688/1‑151 (MQ=255)
                                                                                                                                           aaGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACtt      >  1:462181/1‑151 (MQ=255)
                                                                                                                                           aaGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACt       <  2:1225466/150‑1 (MQ=255)
                                                                                                                                            aGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTc     >  2:693925/1‑151 (MQ=255)
                                                                                                                                            aGAAAATAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTc     <  1:511054/151‑1 (MQ=255)
                                                                                                                                              aaaaTAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTcc    >  1:203333/1‑150 (MQ=255)
                                                                                                                                              aaaaTAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTcc    >  1:1017035/1‑150 (MQ=255)
                                                                                                                                               aaaTAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAg  <  2:1315612/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTGAGGCGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAG  >  NC_000913/1115277‑1115570

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