Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,363,453 C→T 5.1% A61T (GCC→ACC) ‡ yfaY ← inactive PncC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,363,4530CT 5.1% 328.4 / 5.8 97A61T (GCC→ACC) ‡yfaYinactive PncC family protein
Reads supporting (aligned to +/- strand):  ref base C (39/53);  new base T (2/3);  total (41/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

TACTCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTAACACT  >  NC_000913/2363304‑2363598
                                                                                                                                                     |                                                                                                                                                 
tACTCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGc                                                                                                                                                   <  1:1026155/150‑1 (MQ=255)
   tCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTg                                                                                                                                                >  2:1399400/1‑150 (MQ=255)
   tCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTg                                                                                                                                                >  1:994413/1‑150 (MQ=255)
   tCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTg                                                                                                                                                <  1:854454/150‑1 (MQ=255)
   tCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTg                                                                                                                                                >  1:1241798/1‑150 (MQ=255)
     gTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCt                                                                                                                                              <  1:143557/150‑1 (MQ=255)
       ccACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGa                                                                                                                                            <  2:823247/150‑1 (MQ=255)
        cACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGAc                                                                                                                                           >  1:665800/1‑150 (MQ=255)
         aCGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACGCCGATGTGCTGatc                                                                                                                                          <  1:1168460/150‑3 (MQ=255)
         aCGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACGCCGATGTGCTGatc                                                                                                                                          <  1:1168457/150‑3 (MQ=255)
         aCGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACt                                                                                                                                          <  2:860210/150‑1 (MQ=255)
             tCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACGCCGATGTGCTGatcgtta                                                                                                                                      <  2:348882/150‑7 (MQ=255)
               gTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTc                                                                                                                                    >  2:350061/1‑150 (MQ=255)
               gTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTc                                                                                                                                    >  2:582752/1‑150 (MQ=255)
                  gAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTcgcg                                                                                                                                 <  1:565586/150‑1 (MQ=255)
                  gAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTcgcg                                                                                                                                 <  1:741764/150‑1 (MQ=255)
                       agggTTCCATCTCTTTGAGGGAGGCTTCTTGCAGCACCCGGCGGTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCCTCGCTGGTCGGTCCCAGCCCGCCGTTAACGCTCTGCACATCGGCGTGGTGACTACGTTCGCgaagaa                                                                                                                            <  2:1235182/147‑1 (MQ=255)
                       aGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGGCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCgaagaa                                                                                                                            <  2:1235181/150‑1 (MQ=255)
                         cGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATg                                                                                                                          >  1:1168967/1‑150 (MQ=255)
                           ttCCATCTCTTTGTGACTGGCTTCTTGCTGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGTGTGCGCTTAAATAATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTTCTGACTACGTTCGCGAAGAATGGt                                                                                                                        <  2:1392887/150‑1 (MQ=255)
                           ttCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGt                                                                                                                        <  2:11674/150‑1 (MQ=255)
                           ttCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGt                                                                                                                        <  2:520665/150‑1 (MQ=255)
                           ttCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGt                                                                                                                        <  1:825332/150‑1 (MQ=255)
                            tCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTg                                                                                                                       >  2:731164/1‑150 (MQ=255)
                            tCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTg                                                                                                                       >  2:379870/1‑150 (MQ=255)
                             ccATCTCTTTGAGCCAGGCTTCATGCAGCACCCGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGa                                                                                                                      <  2:915593/150‑1 (MQ=255)
                             ccATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGa                                                                                                                      <  2:590028/150‑1 (MQ=255)
                             ccATCTCTTTGAGCCAGGCGTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTCAATCATCGCTGGTCGTTCCCATCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGTAGAATGGTGa                                                                                                                      <  2:775599/150‑1 (MQ=255)
                                tctcTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTa                                                                                                                   >  1:527910/1‑150 (MQ=255)
                                    tttGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGtc                                                                                                               <  2:567138/150‑1 (MQ=255)
                                    tttGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGtc                                                                                                               >  2:808440/1‑150 (MQ=255)
                                    tttGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGtc                                                                                                               >  1:720437/1‑150 (MQ=255)
                                     ttGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGtca                                                                                                              <  2:34870/150‑1 (MQ=255)
                                       gcgCCAGGCTTCCTGCCGCCCCCGGCCTGCACCTCGCGCTGGCGCAGCGGCGAGTGCGCATAACCCATCGCTGGACGGTCCCAGCCCGCCGAGACCGACCAGCACAGCGGCGTGCTGACAACGTTCGCGAAGAATGGTGACTAAGtcatc                                                                                                            <  2:1052265/148‑1 (MQ=255)
                                        aGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGtcatca                                                                                                           >  1:1367445/1‑150 (MQ=255)
                                          ccAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAg                                                                                                         <  1:855756/150‑1 (MQ=255)
                                           cAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAgg                                                                                                        >  2:448726/1‑150 (MQ=255)
                                           cAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAgg                                                                                                        >  2:104605/1‑150 (MQ=255)
                                           cAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAgg                                                                                                        >  2:1133360/1‑150 (MQ=255)
                                           cAGGCTTCATGCAGCACCACGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAgg                                                                                                        <  1:396888/150‑1 (MQ=255)
                                            aGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGt                                                                                                       >  1:442380/1‑150 (MQ=255)
                                            aGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGt                                                                                                       >  1:442416/1‑150 (MQ=255)
                                                ttCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATc                                                                                                   <  1:299907/150‑1 (MQ=255)
                                                  cATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATccc                                                                                                 <  1:135384/150‑1 (MQ=255)
                                                  cATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATccc                                                                                                 <  1:135321/150‑1 (MQ=255)
                                                     gcagcaCCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCa                                                                                              <  2:51981/150‑1 (MQ=255)
                                                             aGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTg                                                                                      >  1:938893/1‑150 (MQ=255)
                                                               gCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcg                                                                                    >  2:1189958/1‑150 (MQ=255)
                                                                ccTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcgg                                                                                   <  1:1160985/150‑1 (MQ=255)
                                                                ccTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcgg                                                                                   <  1:1160931/150‑1 (MQ=255)
                                                                 cTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggc                                                                                  >  1:1042251/1‑150 (MQ=255)
                                                                 cTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggc                                                                                  >  2:848208/1‑150 (MQ=255)
                                                                 cTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggc                                                                                  <  1:339699/150‑1 (MQ=255)
                                                                 cTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggc                                                                                  <  1:504653/150‑1 (MQ=255)
                                                                 cTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggc                                                                                  >  1:584008/1‑150 (MQ=255)
                                                                  ttaaCCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGACGTTAACGATCAGCACATCGGCGTGCTGACTCCGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgaggcg                                                                                 >  1:1266772/4‑150 (MQ=255)
                                                                  ttCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggcg                                                                                 >  1:528212/1‑150 (MQ=255)
                                                                  ttCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggcg                                                                                 >  2:225994/1‑150 (MQ=255)
                                                                  ttCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggcg                                                                                 >  1:1266775/1‑150 (MQ=255)
                                                                         tttGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAAt                                                                          <  1:797783/150‑1 (MQ=255)
                                                                               gTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAAc                                                                    <  1:641663/150‑1 (MQ=255)
                                                                               gTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAAc                                                                    <  1:641594/150‑1 (MQ=255)
                                                                               gTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAAc                                                                    <  1:641575/150‑1 (MQ=255)
                                                                                   cAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCt                                                                <  1:509872/150‑1 (MQ=255)
                                                                                   cAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCt                                                                <  1:509913/150‑1 (MQ=255)
                                                                                     gcggcgAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCtga                                                              >  1:974584/1‑150 (MQ=255)
                                                                                     gcggcgAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCtga                                                              <  2:448035/150‑1 (MQ=255)
                                                                                          gAGTGCGCTTATATCATCGCTGGTCGGTCCCAGCGCGCCGTTAACGATCAGCTCATCGGCGTGCTGACTACGTTCGCGAATAATGGTGACTAAGTCATCTAGGTTATCCCCCACCTTATTGCGGCGATATAATGGCAACCCCTGATgaaa                                                         <  2:793551/150‑1 (MQ=255)
                                                                                          gAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATgaaa                                                         <  1:574012/150‑1 (MQ=255)
                                                                                          gAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATgaaa                                                         <  1:1226747/150‑1 (MQ=255)
                                                                                          gAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATgaaa                                                         <  2:1129755/150‑1 (MQ=255)
                                                                                          gAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATgaaa                                                         <  1:874063/150‑1 (MQ=255)
                                                                                          gAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATgaaa                                                         <  1:236516/150‑1 (MQ=255)
                                                                                                   tAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATc                                                <  1:1206813/150‑1 (MQ=255)
                                                                                                    aaaTCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCg                                               <  2:759619/150‑1 (MQ=255)
                                                                                                     aaTCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCgg                                              <  2:1313329/150‑1 (MQ=255)
                                                                                                         atcGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGgccag                                          >  1:1258251/1‑150 (MQ=255)
                                                                                                            gCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGgccagcca                                       <  2:1128063/150‑1 (MQ=255)
                                                                                                                 tcggtcCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAg                                  >  2:227570/1‑150 (MQ=255)
                                                                                                                 tcggtcCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAg                                  >  2:1124956/1‑150 (MQ=255)
                                                                                                                   ggtcCCTGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTAGGTTCGCTATTATTGGTGACTAAGTCATCAAGGTTATCCCCCACCTTTTTGCGGCGATTTAATGGCAACCCCTTATGAAATAAAAAATCGGCCAGCCAGGCAGCg                                <  1:897027/150‑1 (MQ=255)
                                                                                                                      cccAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTa                             >  2:731377/1‑150 (MQ=255)
                                                                                                                       ccAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAg                            >  1:1236040/1‑150 (MQ=255)
                                                                                                                       ccAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATAGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAg                            >  1:1236050/1‑150 (MQ=255)
                                                                                                                        cAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAgt                           <  2:146587/150‑1 (MQ=255)
                                                                                                                           cccGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGACAAAATCGGCCAGCCAGGCAGCGTTAGTGTc                        >  2:729875/1‑150 (MQ=255)
                                                                                                                                  ttAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGagt                 >  2:1060865/1‑148 (MQ=255)
                                                                                                                                    aaCGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTg               >  1:748539/1‑150 (MQ=255)
                                                                                                                                      tagtCAGCACGCCGATGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGcc             >  2:597425/4‑150 (MQ=255)
                                                                                                                                        gtCAGCACGCCGATGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCg           >  1:498339/2‑150 (MQ=255)
                                                                                                                                        aTCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCTTATTGCGGCGAGATATTGGAAACCCATGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCg           <  1:416598/150‑1 (MQ=255)
                                                                                                                                         tccgCACATGGGCGTGCTGAGTCGGTTCGCGACGAAGGGTGACTAAGTCATCTCGGTTATCCCCCACGGTATTGGGGCGTGTTAATGGCAACCCCTTATTAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCgt          <  1:551194/147‑1 (MQ=255)
                                                                                                                                            gCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTa       >  1:1375087/1‑150 (MQ=255)
                                                                                                                                             cacaTCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTaa      >  1:1144033/1‑150 (MQ=255)
                                                                                                                                               caTCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTAaca    <  2:183249/150‑1 (MQ=255)
                                                                                                                                                aTCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTAacac   <  2:1291705/150‑1 (MQ=255)
                                                                                                                                                aTCGGCGTGCTGACGACGTTCGCGAAGACTGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCTACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGTTTTGCACGTGTAacac   <  2:1291684/150‑1 (MQ=255)
                                                                                                                                                 tCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTAACACt  <  2:731528/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TACTCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTAACACT  >  NC_000913/2363304‑2363598

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