Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,101,039 (C)7→8 coding (585/720 nt) yggN ← DUF2884 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,101,0321.C93.5% 201.0 / 5.8 62coding (592/720 nt)yggNDUF2884 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base C (27/31);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

GAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGGCG  >  NC_000913/3100899‑3101165
                                                                                                                                      |                                                                                                                                       
gagGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑ccccc                                                                                                                                  >  2:1955877/1‑139 (MQ=255)
gagGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑ccccc                                                                                                                                  >  1:1394742/1‑139 (MQ=255)
      tattTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑ccccccccagg                                                                                                                            >  1:421110/1‑135 (MQ=255)
          ttAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTc                                                                                                                         <  1:765071/139‑1 (MQ=255)
                 cccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGcac                                                                                                                  <  1:1581503/139‑1 (MQ=255)
                  ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAGTCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑c                                                                                                                                       >  1:1181972/1‑117 (MQ=255)
                  ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGcaca                                                                                                                 <  2:1366582/139‑1 (MQ=255)
                  ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGcaca                                                                                                                 >  2:635845/1‑139 (MQ=255)
                  ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑c                                                                                                                                       <  2:1181972/117‑1 (MQ=255)
                   cGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCCAGGCTTCCCAGcaca                                                                                                                 >  2:1347569/1‑139 (MQ=255)
                   cGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACAt                                                                                                                <  1:635845/139‑1 (MQ=255)
                    gACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACAtt                                                                                                               <  2:421110/139‑1 (MQ=255)
                          ggCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGc                                                                                                                    <  2:404665/128‑1 (MQ=255)
                          ggCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTaa                                                                                                         >  1:2042398/1‑139 (MQ=255)
                          ggCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTaa                                                                                                         >  1:1287876/1‑139 (MQ=255)
                          ggCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cccccc                                                                                                                                 <  2:419945/114‑1 (MQ=255)
                          ggCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cccccc                                                                                                                                 >  1:419945/1‑114 (MQ=255)
                          ggCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑CCCCACCCAGGCTTCCCAGc                                                                                                                    >  1:404665/1‑128 (MQ=255)
                               tGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGt                                                                                                    >  2:1580107/1‑139 (MQ=255)
                               tGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGt                                                                                                    >  1:937708/1‑139 (MQ=255)
                                gCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGtt                                                                                                   >  1:1025551/1‑139 (MQ=255)
                                    cTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccg                                                                                               >  2:995154/1‑139 (MQ=255)
                                      aTCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑ccccccccagg                                                                                                                            >  1:916177/1‑103 (MQ=255)
                                      aTCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑ccccccccagg                                                                                                                            <  2:916177/107‑5 (MQ=255)
                                      aTCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cccccc                                                                                                                                 <  1:224907/102‑1 (MQ=255)
                                      aTCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cccccc                                                                                                                                 >  2:224907/1‑102 (MQ=255)
                                         ttCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTccc                                                                                                                       <  1:244577/110‑1 (MQ=255)
                                         ttCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTccc                                                                                                                       >  2:244577/1‑110 (MQ=255)
                                            cAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTc                                                                                      >  1:1065215/1‑139 (MQ=255)
                                                gTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCa                                                                                  <  2:1065215/139‑1 (MQ=255)
                                                 tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCa                                                                                  >  2:208025/1‑139 (MQ=255)
                                                   acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAcc                                                                                >  1:1471689/1‑139 (MQ=255)
                                                   acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAcc                                                                                >  2:1594346/1‑139 (MQ=255)
                                                   acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAcc                                                                                >  1:381743/1‑139 (MQ=255)
                                                         cgGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTc                                                                          >  2:1601031/1‑139 (MQ=255)
                                                            cTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcgcg                                                                       <  2:1025551/139‑1 (MQ=255)
                                                               cAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCa                                                                   <  1:1955877/139‑1 (MQ=255)
                                                                 aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  2:937708/139‑1 (MQ=255)
                                                                 aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  1:1580107/139‑1 (MQ=255)
                                                                 aaCATCTTTGCCAAACAGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTTAGCACCGCTTTCGCGCCCAt                                                                  <  2:1394742/139‑1 (MQ=255)
                                                                  aCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCTCGCCCAt                                                                  <  1:1347569/139‑1 (MQ=255)
                                                                      cTTTGCCAAACTGCTTGAAATTTTTTTCCTCCTTTTTCCACGCGGTTTGGATTCAGGCTTGCAGCCCCCCCCCGGCTTCCCAGCACATTCTGTAATGGGTTCCCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAt                                                             <  1:953997/139‑1 (MQ=255)
                                                                      cTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAt                                                             <  1:1594346/139‑1 (MQ=255)
                                                                         tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGc                                                                                                                    >  1:834809/1‑81 (MQ=255)
                                                                         tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGc                                                                                                                    <  2:834809/81‑1 (MQ=255)
                                                                         tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTaa                                                          >  1:1441410/1‑139 (MQ=255)
                                                                          gCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCa                                                                                                                      >  2:911081/1‑78 (MQ=255)
                                                                          gCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCa                                                                                                                      <  1:911081/78‑1 (MQ=255)
                                                                          gCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCCCTTTCGCGCCCATTTCATTatt                                                         >  2:953997/1‑137 (MQ=255)
                                                                           ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATg                                                        <  2:1287876/139‑1 (MQ=255)
                                                                           ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATg                                                        <  2:381743/139‑1 (MQ=255)
                                                                             aaaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCt                                                      <  2:1471689/139‑1 (MQ=255)
                                                                              aaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  >  2:137388/1‑126 (MQ=255)
                                                                              aaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  >  1:332282/1‑126 (MQ=255)
                                                                              aaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  1:137388/126‑1 (MQ=255)
                                                                              aaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  2:332282/126‑1 (MQ=255)
                                                                                   ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑ccccccccagg                                                                                                                            >  2:174597/1‑58 (MQ=255)
                                                                                   ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑ccccccccagg                                                                                                                            <  1:174597/62‑5 (MQ=255)
                                                                                     tgAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTa                                              <  1:747309/138‑1 (MQ=255)
                                                                                             ttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cc                                                                                                                                      <  1:835643/43‑1 (MQ=255)
                                                                                             ttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cc                                                                                                                                      >  2:835643/1‑43 (MQ=255)
                                                                                                      tttttCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGcc                             <  1:995154/139‑1 (MQ=255)
                                                                                                               tCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt                    <  2:1441410/139‑1 (MQ=255)
                                                                                                                cGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                                            >  1:536592/1‑66 (MQ=255)
                                                                                                                cGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                                            <  2:536592/66‑1 (MQ=255)
                                                                                                                cGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATg                                                        <  2:854449/102‑1 (MQ=255)
                                                                                                                cGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATg                                                        >  1:854449/1‑102 (MQ=255)
                                                                                                                           gAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAAGGGTTACCGCCGCTTTTCAGCACCGCTTTc                                                                          <  1:1403311/72‑1 (MQ=255)
                                                                                                                           gAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAAGGGTTACCGCCGCTTTTCAGCACCGCTTTc                                                                          >  2:1403311/1‑72 (MQ=255)
                                                                                                                            aGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGcc                             >  2:445285/1‑117 (MQ=255)
                                                                                                                            aGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGcc                             <  1:445285/117‑1 (MQ=255)
                                                                                                                            aGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCtt       >  2:232963/1‑139 (MQ=255)
                                                                                                                                 tGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTaa                                                                                                         <  2:182236/36‑1 (MQ=255)
                                                                                                                                 tGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTaa                                                                                                         >  1:182236/1‑36 (MQ=255)
                                                                                                                                 tGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  1:446837/75‑1 (MQ=255)
                                                                                                                                 tGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  >  2:446837/1‑75 (MQ=255)
                                                                                                                                 tGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGgcg  <  2:2042398/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
GAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGGCG  >  NC_000913/3100899‑3101165

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