Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,423,776 G→T 23.8% noncoding (12/120 nt) rrfD ← 5S ribosomal RNA of rrnD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,7760GT23.8% 217.8 / 55.6 99noncoding (12/120 nt)rrfD5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base G (19/56);  new base T (2/22);  total (21/78)
Fisher's exact test for biased strand distribution p-value = 9.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATC  >  NC_000913/3423700‑3423846
                                                                            |                                                                      
gACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTA                                                 >  1:2690779‑M1/1‑100 (MQ=255)
gACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTA                                                 <  2:2161906‑M1/100‑1 (MQ=255)
gACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTA                                                 <  1:83714‑M1/100‑1 (MQ=255)
 aCCCCACACTACCATCGGCGCTACTGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTAT                                                <  1:696278‑M1/100‑1 (MQ=255)
 aCCCCACACTACCATCGGCGCTACTGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTAT                                                >  1:126827‑M1/1‑100 (MQ=255)
 aCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTAT                                                <  2:2112433‑M1/100‑1 (MQ=255)
  ccccACACTACCATCGGCGCTACTGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATC                                               >  1:2315659‑M1/1‑100 (MQ=255)
  ccccACACTACCATCGGCGCTACTGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATC                                               >  2:2078672‑M1/1‑100 (MQ=255)
  ccccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATC                                               <  1:527800‑M1/100‑1 (MQ=255)
  ccccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATC                                               >  1:1597001‑M1/1‑100 (MQ=255)
  ccccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATC                                               >  1:3481187‑M1/1‑100 (MQ=255)
  ccccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATC                                               <  2:473219‑M1/100‑1 (MQ=255)
  ccccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATC                                               <  1:1031692‑M1/100‑1 (MQ=255)
  ccccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATC                                               >  1:1489306‑M1/1‑100 (MQ=255)
  ccccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATC                                               >  1:2100341‑M1/1‑100 (MQ=255)
  ccccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATC                                               <  2:1114989‑M1/100‑1 (MQ=255)
  ccccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATC                                               >  1:2807761‑M1/1‑100 (MQ=255)
  ccccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATC                                               >  2:139535‑M1/1‑100 (MQ=255)
   cccACACTACCATCGGCGCTACTGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCA                                              <  2:3220584‑M1/100‑1 (MQ=255)
   cccACACTACCATCGGCGCTACTGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCA                                              >  1:353210‑M1/1‑100 (MQ=255)
   cccACACTACCATCGGCGCTACTGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCA                                              <  1:2003973‑M1/100‑1 (MQ=255)
   cccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCA                                              >  1:939818‑M1/1‑100 (MQ=255)
   cccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCA                                              >  2:1961720‑M1/1‑100 (MQ=255)
   cccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCA                                              <  1:1309116‑M1/100‑1 (MQ=255)
   cccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCA                                              <  1:1653681‑M1/100‑1 (MQ=255)
   ccaacactacCATCGGCGCTACTGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTAT                                                >  1:2429298‑M1/4‑98 (MQ=255)
    ccACACTACCATCGGCGCTACTGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAG                                             <  2:3401192‑M1/100‑1 (MQ=255)
    ccACACTACCATCGGCGCTACTGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAG                                             <  1:2558584‑M1/100‑1 (MQ=255)
    ccACACTACCATCGGCGCTACTGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAG                                             <  2:343605‑M1/100‑1 (MQ=255)
    ccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCCTGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAG                                             >  2:3675087‑M1/1‑100 (MQ=255)
    ccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAG                                             <  2:3644306‑M1/100‑1 (MQ=255)
    ccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAG                                             >  1:1859260‑M1/1‑100 (MQ=255)
    ccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTAAGTTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAG                                             <  2:560033‑M1/100‑1 (MQ=255)
     cacaCTACCATCGGCGCTACTGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGA                                            >  1:3740415‑M1/1‑100 (MQ=255)
     cacaCTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGA                                            >  1:1687081‑M1/1‑100 (MQ=255)
      acacTACCATCGGCGCTACTGCGTTTCACTTCGGAGTTCGGCATGGGGTCAGGTTGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGA                                            >  2:3749574‑M1/1‑99 (MQ=255)
              aTCGGCGCTACTGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaag                                   <  1:62547‑M1/100‑10 (MQ=255)
                   cgcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtc                              <  1:1032321‑M1/100‑15 (MQ=255)
                      tACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctga                           <  1:563784‑M1/100‑18 (MQ=255)
                       actgcgtTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgat                          <  1:2181094‑M1/97‑19 (MQ=255)
                       aCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgat                          <  1:1400853‑M1/100‑19 (MQ=255)
                       aCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgat                          <  2:647242‑M1/100‑19 (MQ=255)
                       aCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgat                          <  1:973792‑M1/100‑19 (MQ=255)
                       aCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgat                          <  1:530083‑M1/100‑19 (MQ=255)
                        cGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgata                         <  1:2225596‑M1/100‑20 (MQ=255)
                        cGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgata                         <  1:3386892‑M1/100‑20 (MQ=255)
                        cGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgata                         <  2:3851506‑M1/100‑20 (MQ=255)
                        cGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgata                         <  1:790512‑M1/100‑20 (MQ=255)
                        cGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgata                         <  2:1264569‑M1/100‑20 (MQ=255)
                        cGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTCATCAGAacgcagaagcggtctgatg                         <  2:4083894‑M1/100‑20 (MQ=255)
                        cGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTCATCAGAacgcagaagcggtctgatg                         <  2:3592142‑M1/100‑20 (MQ=255)
                          gtgttTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaa                       <  2:2493270‑M1/98‑22 (MQ=255)
                          gCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaa                       <  1:2391658‑M1/100‑22 (MQ=255)
                          gCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaa                       <  1:3498777‑M1/100‑22 (MQ=255)
                          gCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaa                       <  1:1856487‑M1/100‑22 (MQ=255)
                          gCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTCATCAGAacgcagaagcggtctgatgaa                       <  1:3991682‑M1/100‑22 (MQ=255)
                          gCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTCATCAGAacgcagaagcggtctgatgaa                       <  1:1929224‑M1/100‑22 (MQ=255)
                           cGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaa                      <  2:170921‑M1/100‑23 (MQ=255)
                           cGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaa                      <  1:3141497‑M1/100‑23 (MQ=255)
                           cGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTCATCAGAacgcagaagcggtctgatgaaa                      <  1:450389‑M1/100‑23 (MQ=255)
                            gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaac                     <  2:2829741‑M1/100‑24 (MQ=255)
                            gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaac                     <  1:4126712‑M1/100‑24 (MQ=255)
                            gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaac                     <  1:3327334‑M1/100‑24 (MQ=255)
                            gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaac                     <  1:3134464‑M1/100‑24 (MQ=255)
                            gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaac                     <  1:2346093‑M1/100‑24 (MQ=255)
                            gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaac                     <  2:1158147‑M1/100‑24 (MQ=255)
                            gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaac                     <  1:1205608‑M1/100‑24 (MQ=255)
                            gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTCTATCAGAacgcagaagcggtctgatagagg                     >  2:3498750‑M1/1‑77 (MQ=255)
                            gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaac                     <  1:882174‑M1/100‑24 (MQ=255)
                            gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaac                     <  1:2900164‑M1/100‑24 (MQ=255)
                            gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaac                     <  1:1640021‑M1/100‑24 (MQ=255)
                            gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaac                     <  1:1059703‑M1/100‑24 (MQ=255)
                            gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaac                     <  2:1667922‑M1/100‑24 (MQ=255)
                            gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTCATCAGAacgcagaagcggtctgatgaaac                     <  2:907961‑M1/100‑24 (MQ=255)
                            gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTCATCAGAacgcagaagcggtctgatgaaac                     <  1:257576‑M1/100‑24 (MQ=255)
                            gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTCATCAGAacgcagaagcggtctgatgaaac                     <  1:4129049‑M1/100‑24 (MQ=255)
                            gTTTCACTTCTGAGTTCGGAATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTTTGTTTTATCAGAacgcagaagcggtctgataaaac                     <  2:2328108‑M1/100‑24 (MQ=255)
                             tttCTCTTCTGAGTTCGGCATGGGGTCAGGTGGGACCAACGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataagaca                    <  2:1142947‑M1/100‑25 (MQ=255)
                             tttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaaca                    <  1:1939759‑M1/100‑25 (MQ=255)
                             tttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaaca                    <  1:2047537‑M1/100‑25 (MQ=255)
                              ttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCCCCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaacag                   <  2:188068‑M1/100‑26 (MQ=255)
                              ttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaacag                   <  1:1313336‑M1/100‑26 (MQ=255)
                              ttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaacag                   <  1:939395‑M1/100‑26 (MQ=255)
                              ttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaacag                   <  1:2589130‑M1/100‑26 (MQ=255)
                              ttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaacag                   <  1:499729‑M1/100‑26 (MQ=255)
                              ttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaacag                   <  1:4149693‑M1/100‑26 (MQ=255)
                              ttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaacag                   <  1:3038212‑M1/100‑26 (MQ=255)
                              ttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTCATCAGAacgcagaagcggtctgatgaaacag                   <  1:3139213‑M1/100‑26 (MQ=255)
                              ttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTCATCAGAacgcagaagcggtctgatgaaacag                   <  2:3044029‑M1/100‑26 (MQ=255)
                              ttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTCATCAGAacgcagaagcggtctgatgaaacag                   <  2:153094‑M1/100‑26 (MQ=255)
                               tCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaacaga                  <  1:1926574‑M1/100‑27 (MQ=255)
                               tCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaacaga                  <  1:1569775‑M1/100‑27 (MQ=255)
                                cACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaacagaa                 <  1:900746‑M1/100‑28 (MQ=255)
                                cACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaacaga                  <  1:1815438‑M1/99‑27 (MQ=255)
                                cACTTCTGAGTTCGGCATGGCGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTATGTTTTATCAGAacgcagaagcggtctgataaaacagaa                 <  2:3613566‑M1/100‑28 (MQ=255)
                                    tCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaacag                   <  2:3974031‑M1/94‑26 (MQ=255)
                                      tGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTGTTTTACCAGAacgcagaagcggtctgataaaacaga                  >  2:1059699‑M1/1‑67 (MQ=255)
                                      tGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaacagaatttgcc           <  1:4102857‑M1/100‑34 (MQ=255)
                                               cATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGAacgcagaagcggtctgataaaacagaatttgcctggcggcag  <  1:1549492‑M1/100‑43 (MQ=255)
                                                                            |                                                                      
GACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATC  >  NC_000913/3423700‑3423846

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