Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 103,574 G→T 6.7% A140A (GCG→GCT ftsQ → cell division protein FtsQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913103,5740GT 6.7% 198.0 / 3.6 60A140A (GCG→GCTftsQcell division protein FtsQ
Reads supporting (aligned to +/- strand):  ref base G (29/27);  new base T (2/2);  total (31/29)
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.69e-01

TGTCAACATCATCCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAA  >  NC_000913/103427‑103720
                                                                                                                                                   |                                                                                                                                                  
tGTCAACATCATCCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCgga                                                                                                                                                  <  1:1420034/150‑1 (MQ=255)
  tCAACATCATCCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCggaag                                                                                                                                                >  2:310934/1‑150 (MQ=255)
   cAACATCATCCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCggaagg                                                                                                                                               >  1:1368923/1‑150 (MQ=255)
    aaCATCATCCAGACGCAAATAGATCAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCggaagga                                                                                                                                              <  2:997994/150‑1 (MQ=255)
    aaCATCATCCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCggaagga                                                                                                                                              >  1:1079982/1‑150 (MQ=255)
     aCATCATCCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCggaaggaa                                                                                                                                             >  1:1315011/1‑150 (MQ=255)
       atcatcCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAAt                                                                                                                                           >  1:445568/1‑150 (MQ=255)
        tcatcCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATa                                                                                                                                          <  2:5900/150‑1 (MQ=255)
         catcCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATAc                                                                                                                                         <  2:1315011/150‑1 (MQ=255)
          atcCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATAcc                                                                                                                                        <  1:552943/150‑1 (MQ=255)
           tcCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCt                                                                                                                                       >  1:629526/1‑150 (MQ=255)
             cAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTc                                                                                                                                     >  1:1291386/1‑150 (MQ=255)
                 cGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCTGAAGGTATTTCCTTCcgcg                                                                                                                                 <  2:528739/150‑1 (MQ=255)
                     aaTAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCTGAAGGTATTTCCTTCCGCGtcta                                                                                                                             <  1:62461/150‑4 (MQ=255)
                      aTAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTgccg                                                                                                                            <  2:445568/150‑1 (MQ=255)
                              acGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACg                                                                                                                    <  1:293138/150‑1 (MQ=255)
                                 ccTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCAc                                                                                                                 >  2:257981/1‑150 (MQ=255)
                                      cGTGGATTAAGCAGGTGATCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCagca                                                                                                            <  2:396394/150‑1 (MQ=255)
                                          gATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCagcaagca                                                                                                        >  2:32550/1‑150 (MQ=255)
                                          gATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCagcaagca                                                                                                        <  1:438448/150‑1 (MQ=255)
                                          gATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCagcaagca                                                                                                        >  1:753512/1‑150 (MQ=255)
                                          gATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCagcaagca                                                                                                        <  2:254565/150‑1 (MQ=255)
                                            ttAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAgg                                                                                                      >  2:214333/1‑150 (MQ=255)
                                              aaGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTg                                                                                                    <  2:1291386/150‑1 (MQ=255)
                                                gCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCt                                                                                                  <  2:629526/150‑1 (MQ=255)
                                                   ggTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTcc                                                                                               <  2:1175797/150‑1 (MQ=255)
                                                     tGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCaa                                                                                             >  2:533521/1‑150 (MQ=255)
                                                      gAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAAt                                                                                            <  2:862877/150‑1 (MQ=255)
                                                        gCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGc                                                                                          <  2:393271/150‑1 (MQ=255)
                                                         cGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCt                                                                                         >  1:1034274/1‑150 (MQ=255)
                                                          gTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGAAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTg                                                                                        <  2:17860/150‑1 (MQ=255)
                                                           tCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGt                                                                                       >  2:711002/1‑150 (MQ=255)
                                                            cAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCTTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTa                                                                                      >  2:901004/1‑150 (MQ=255)
                                                                       ggCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAg                                                                           >  2:161978/1‑150 (MQ=255)
                                                                        gCCTTATGAATTAAAGATTAATCTGGTTGAATATGTGAAGATTGCGCGGTGGAATGATCAACATATGGTGGAATCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAgg                                                                          >  2:838738/1‑150 (MQ=255)
                                                                           tgatgaATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAg                                                                       >  1:1132967/1‑150 (MQ=255)
                                                                            gatgaATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAAAATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCAGAAGGCAgc                                                                      >  2:1041798/1‑150 (MQ=255)
                                                                               gaATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGcc                                                                   <  2:1327797/150‑1 (MQ=255)
                                                                               gaATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGcc                                                                   <  1:533521/150‑1 (MQ=255)
                                                                                 aTTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCaa                                                                 >  2:683394/1‑150 (MQ=255)
                                                                                      aGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAg                                                            <  2:346390/150‑1 (MQ=255)
                                                                                       gATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgt                                                           <  2:906022/150‑1 (MQ=255)
                                                                                       gATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgt                                                           <  1:683394/150‑1 (MQ=255)
                                                                                              ctcgTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAACGAAATACCTTCAGCGTGCCGCCACAACGCCCAAGCAAGCAGGTGCTTCCAATCCTGTATGGCCCGGAAGGCACCGCCAATCAAGTGTTGCAg                                                    >  2:112772/4‑150 (MQ=255)
                                                                                                ggTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCgg                                                                              >  1:500839/1‑122 (MQ=255)
                                                                                                ggTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCgg                                                                              <  2:500839/122‑1 (MQ=255)
                                                                                                   tGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTa                                               <  1:711002/150‑1 (MQ=255)
                                                                                                     aaTATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATc                                             >  1:71289/1‑150 (MQ=255)
                                                                                                            gCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAAt                                      >  1:713678/1‑150 (MQ=255)
                                                                                                            gCCGATTGCGCGGTGGAATGATCAACATATCGTAGACGCGGAAGTAAATACATAAAGCGTGACGCCATAACGCACCAGCAAGCAGGTGCTTACAATGCTGTATTGCCCGGAAGCAAGCGCAAATGAAGTGTTGCAGGGCTATAGCGAAAt                                      >  2:1394561/1‑150 (MQ=255)
                                                                                                             ccGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATg                                     <  2:1368923/150‑1 (MQ=255)
                                                                                                              cGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATgg                                    <  1:214333/150‑1 (MQ=255)
                                                                                                                 ttGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGc                                 <  2:753512/150‑1 (MQ=255)
                                                                                                                               gATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAg                   >  2:1402160/1‑150 (MQ=255)
                                                                                                                                aTCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGc                                 >  1:1051473/1‑135 (MQ=255)
                                                                                                                                aTCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGc                                 <  2:1051473/135‑1 (MQ=255)
                                                                                                                                aTCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAgg                  >  2:1181573/1‑150 (MQ=255)
                                                                                                                                aTCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAgg                  >  1:212747/1‑150 (MQ=255)
                                                                                                                                aTCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAgg                  >  1:1268056/1‑150 (MQ=255)
                                                                                                                                   aaCATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACa               >  1:1316080/1‑150 (MQ=255)
                                                                                                                                     cATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAgg                  <  2:445023/145‑1 (MQ=255)
                                                                                                                                     cATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAgg                  >  1:445023/1‑145 (MQ=255)
                                                                                                                                         ggcggcaCGCTGAAGGTATTTCCTTCCGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTa         >  1:1353551/7‑150 (MQ=255)
                                                                                                                                           gTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTAct       >  1:324085/1‑150 (MQ=255)
                                                                                                                                            tAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTActc      <  2:270464/150‑1 (MQ=255)
                                                                                                                                             gcaCGCTGAAGGTATTTCCTTCCGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTActct     >  1:195432/3‑150 (MQ=255)
                                                                                                                                                cgcgGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTgaa  >  2:861853/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGTCAACATCATCCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAA  >  NC_000913/103427‑103720

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: