Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,082,702 A→C 6.0% intergenic (+31/+175) loiP → / ← speB metalloprotease LoiP/agmatinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,082,7020AC 6.0% 251.0 / 7.7 83intergenic (+31/+175)loiP/speBmetalloprotease LoiP/agmatinase
Reads supporting (aligned to +/- strand):  ref base A (39/39);  new base C (2/3);  total (41/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 = 6.40e-01

TTGAAAAACTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGT  >  NC_000913/3082553‑3082849
                                                                                                                                                     |                                                                                                                                                   
ttGAAAAACTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATaa                                                                                                                                                     >  1:781739/1‑150 (MQ=255)
 tGAAAAACTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAg                                                                                                                                                    >  1:445859/1‑150 (MQ=255)
  gAAAAACTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAgg                                                                                                                                                   >  1:548724/1‑150 (MQ=255)
                 cTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATc                                                                                                                                    <  1:1492215/150‑1 (MQ=255)
                 cTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATc                                                                                                                                    >  1:1423482/1‑150 (MQ=255)
                   ggAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATGCGGCGTGAACGCCTTATCTg                                                                                                                                  <  1:12875/150‑1 (MQ=255)
                   ggAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTg                                                                                                                                  <  1:1132239/150‑1 (MQ=255)
                   ggAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTg                                                                                                                                  >  1:1164984/1‑150 (MQ=255)
                    gaagaaGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGa                                                                                                                                 <  2:769617/150‑1 (MQ=255)
                     aagaagGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATGCGGCGTGAACGCCTTATCTGAc                                                                                                                                <  2:617341/150‑1 (MQ=255)
                     aagaagGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGAc                                                                                                                                <  2:1249120/150‑1 (MQ=255)
                     aagaagGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGAc                                                                                                                                <  2:1224936/150‑1 (MQ=255)
                         agGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATGCGGCGTGAACGCCTTATCTGAcctac                                                                                                                            <  1:1363843/150‑5 (MQ=255)
                         agGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGa                                                                                                                            >  1:1294653/1‑150 (MQ=255)
                         agGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGa                                                                                                                            >  1:421253/1‑150 (MQ=255)
                          ggTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGaa                                                                                                                           <  1:974157/150‑1 (MQ=255)
                           gTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAAGGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAACGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAAt                                                                                                                          <  1:1130554/150‑1 (MQ=255)
                             cGCCTAAGCTCAATGTTTGACGACCATCCTGCATCGGCCGAACGCGCCCACCATATTCGCGATCTCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATgg                                                                                                                        <  2:981517/150‑1 (MQ=255)
                              gCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGc                                                                                                                       <  1:410574/150‑1 (MQ=255)
                               ccAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGca                                                                                                                      >  1:789831/1‑150 (MQ=255)
                                 aaaGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGcaca                                                                                                                    <  1:54173/150‑1 (MQ=255)
                                    gCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACAttt                                                                                                                 <  2:445859/150‑1 (MQ=255)
                                     cTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTg                                                                                                                >  2:201061/1‑150 (MQ=255)
                                       cAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGATCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATTGCACATTTGtt                                                                                                              <  2:1511983/150‑1 (MQ=255)
                                       cAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGtt                                                                                                              >  2:573999/1‑150 (MQ=255)
                                         aTGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTAc                                                                                                            <  2:781739/150‑1 (MQ=255)
                                          tGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACGTTTGTTAcc                                                                                                           >  1:765471/1‑150 (MQ=255)
                                          tGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTAcc                                                                                                           >  2:530646/1‑150 (MQ=255)
                                           gTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCt                                                                                                          >  1:592822/1‑150 (MQ=255)
                                             ttGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTtg                                                                                                        <  1:855326/150‑1 (MQ=255)
                                              tGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTtgt                                                                                                       >  2:505705/1‑150 (MQ=255)
                                               gacgacCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTtgtg                                                                                                      >  1:845775/1‑150 (MQ=255)
                                                   acCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCa                                                                                                  >  1:1327296/1‑150 (MQ=255)
                                                         cTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCg                                                                                            <  1:416683/150‑1 (MQ=255)
                                                           gCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGga                                                                                          >  1:739752/1‑150 (MQ=255)
                                                           gCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGga                                                                                          <  1:1407978/150‑1 (MQ=255)
                                                              tCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatgc                                                                                       >  1:498391/1‑150 (MQ=255)
                                                                cgccgAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatgcga                                                                                     >  2:24383/1‑150 (MQ=255)
                                                                   cgAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGAGGGGAATAAGTAAAGCCGGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatgcgatgc                                                                                  <  2:322996/150‑1 (MQ=255)
                                                                      aCGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTgg                                                                               <  1:847489/150‑1 (MQ=255)
                                                                        gcgcCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGTCACATTTGTTACCTTGTGCGCATTGCCGGATGCGTTGCTGgcg                                                                             <  2:1494963/150‑1 (MQ=255)
                                                                            ccAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATc                                                                         >  1:328947/1‑150 (MQ=255)
                                                                               gCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTa                                                                      >  2:766010/1‑150 (MQ=255)
                                                                                  tatTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATcc                                                                   >  1:245003/1‑150 (MQ=255)
                                                                                  tatTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATcc                                                                   <  2:1423482/150‑1 (MQ=255)
                                                                                  tatTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATcc                                                                   <  1:1258519/150‑1 (MQ=255)
                                                                                      cgcgATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGcc                                                               >  2:1187364/1‑150 (MQ=255)
                                                                                      cgcgATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGcc                                                               <  2:421253/150‑1 (MQ=255)
                                                                                      cgcgATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGcc                                                               <  2:328947/150‑1 (MQ=255)
                                                                                       gcgATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCt                                                              >  1:1185875/1‑150 (MQ=255)
                                                                                            cGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACggg                                                         <  2:568232/150‑1 (MQ=255)
                                                                                            cGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACggg                                                         >  1:568232/1‑150 (MQ=255)
                                                                                                tGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGTACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGcc                                                     <  2:1294653/150‑1 (MQ=255)
                                                                                                tGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGcc                                                     >  2:1213883/1‑150 (MQ=255)
                                                                                                tGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGcc                                                     <  1:201061/150‑1 (MQ=255)
                                                                                                tGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGcc                                                     <  1:639290/150‑1 (MQ=255)
                                                                                                tGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGATGGCGCATCTTATTCGGCCTACGGGTGcc                                                     <  2:1442864/150‑1 (MQ=255)
                                                                                                 gAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCg                                                    <  2:450700/150‑1 (MQ=255)
                                                                                                      cgGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGt                                               <  2:932085/150‑1 (MQ=255)
                                                                                                       ggATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTa                                              <  2:245003/150‑1 (MQ=255)
                                                                                                         aTGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAgg                                            >  1:1031164/1‑150 (MQ=255)
                                                                                                                      aaGCATGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGt                               >  2:932222/1‑150 (MQ=255)
                                                                                                                       aGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGtt                              <  2:356622/150‑1 (MQ=255)
                                                                                                                        gCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTc                             <  1:649575/150‑1 (MQ=255)
                                                                                                                          gttgtgttgTCGAACTTAGGTCATTTAAGGTGTTCTCGGCGCATCTGACATGATTTGCACTTTTGTTACCTTGTGCTAATTGCAGGATGAGATGCTTGCGTATCTTTTCCGGGCTACGTGTGGCGAACGTAGGTCGGATAAGGCGTTCAc                           <  2:765471/147‑1 (MQ=255)
                                                                                                                           cggccgtgcCGACCGCAGGTCAGATAACCCGTTCCCGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCAcg                          <  2:1475428/145‑1 (MQ=38)
                                                                                                                            ggtggtGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCAcgc                         <  2:592822/150‑1 (MQ=255)
                                                                                                                             gtggtgTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCAcgcc                        >  2:81223/1‑150 (MQ=255)
                                                                                                                             gtggtgTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCAcgcc                        >  1:928290/1‑150 (MQ=255)
                                                                                                                              tggtgTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCAcgccg                       >  2:1298992/1‑150 (MQ=255)
                                                                                                                              tggtgTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCAcgccgc                      <  2:461583/150‑1 (MQ=255)
                                                                                                                                 tgtCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTTGCGAATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCAGATAAGGCGTTCACGCCGCAt                    >  2:661289/1‑150 (MQ=255)
                                                                                                                                    cGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCg                 >  2:287291/1‑150 (MQ=255)
                                                                                                                                    cGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCg                 >  1:570524/1‑150 (MQ=255)
                                                                                                                                      aaCGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGAc               <  2:1327296/150‑1 (MQ=255)
                                                                                                                                        cGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAt             <  2:1331363/150‑1 (MQ=255)
                                                                                                                                          catgTCAGATGCGGCGTGAACGCCTTATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTa           >  2:1455836/4‑150 (MQ=38)
                                                                                                                                           atgTCAGATGCGGCGTGAACGCCTTATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTaa          >  1:1033756/3‑150 (MQ=38)
                                                                                                                                            ggTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGTCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAAt         >  2:471462/1‑150 (MQ=255)
                                                                                                                                            ggTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAAt         >  1:470773/1‑150 (MQ=255)
                                                                                                                                            ggTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAAt         >  2:76375/1‑150 (MQ=255)
                                                                                                                                              tCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACTCCGCATCCGACATTAATgg       >  1:937423/1‑150 (MQ=255)
                                                                                                                                                   tAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGt  >  2:807169/1‑150 (MQ=255)
                                                                                                                                                   tAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGt  >  2:117303/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTGAAAAACTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGT  >  NC_000913/3082553‑3082849

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