Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,340,547 A→G 74.4% E102G (GAA→GGA)  yciS → DUF1049 family inner membrane protein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,340,5470AG74.4% 222.0 / 80.9 134E102G (GAA→GGA) yciSDUF1049 family inner membrane protein, function unknown
Reads supporting (aligned to +/- strand):  ref base A (16/18);  new base G (31/67);  total (47/87)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

Only 100 of 135 total aligned reads displayed.

TTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAACAAGCA  >  NC_000913/1340448‑1340643
                                                                                                   |                                                                                                
tttCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGa                                                                                                  >  1:2051136/1‑100 (MQ=255)
  tCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAAt                                                                                                >  1:1517758/1‑100 (MQ=255)
    ccTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATaa                                                                                              >  1:2464617/1‑100 (MQ=255)
     cTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATaa                                                                                              <  2:2603246/99‑1 (MQ=255)
         cgcgcgcTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTc                                                                                         <  1:2187446/100‑1 (MQ=255)
         cgcgcgcTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTc                                                                                         <  2:3967579/100‑1 (MQ=255)
         cgcgcgcTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTc                                                                                         <  1:2841114/100‑1 (MQ=255)
         cgcgcgcTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTc                                                                                         >  2:3562127/1‑100 (MQ=255)
              gcTGAACGTAAAATAAAGCGACTGGAAAACCAGCGTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTa                                                                                       <  2:2898585/97‑1 (MQ=255)
               cTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCt                                                                                   <  1:3187282/100‑1 (MQ=255)
                tGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTg                                                                                  >  1:24866/1‑100 (MQ=255)
                 gAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTgg                                                                                 <  1:1444865/100‑1 (MQ=255)
                 gAACGTAAAACAAAGCGAGTGGAAAACCAGCGTTCACCCGCGACTGACGTGGCTGTATTGCCACACTCGTCAGCGGCGAAGGTATAACTTTCTa                                                                                       <  2:3874372/94‑1 (MQ=255)
                  aaCGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTg                                                                                  <  1:882937/98‑1 (MQ=255)
                   aCGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAg                                                                               >  1:2924708/1‑100 (MQ=255)
                    cGTACATTAAAGCGACTGGAAAACCAGCTTTCACCCGCTACTGACGTGTCTGTAGTGCCGCACTCGTCAGCGGCGAAGTGATAACTTTCTATGCTGGan                                                                               <  2:2152680/99‑2 (MQ=255)
                    cGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAgt                                                                              <  1:3243654/100‑1 (MQ=255)
                    cGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAgt                                                                              <  1:1033760/100‑1 (MQ=255)
                    cGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAgt                                                                              <  1:747617/100‑1 (MQ=255)
                    cGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAg                                                                               >  2:3882981/1‑99 (MQ=255)
                     gTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAgtt                                                                             >  1:1271897/1‑100 (MQ=255)
                      tAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttg                                                                            >  2:840967/1‑100 (MQ=255)
                        aaaTAAAGCGACTGGAAAACCAGCTTTCACACGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAgttgtt                                                                          <  1:1041369/100‑1 (MQ=255)
                            aaaGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGttt                                                                      <  1:1974469/100‑1 (MQ=255)
                             aaGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTc                                                                     <  1:2020348/100‑1 (MQ=255)
                             aaGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTc                                                                     <  2:2838837/100‑1 (MQ=255)
                             aaGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTc                                                                     >  1:3606395/1‑100 (MQ=255)
                              aGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCt                                                                    >  1:1939159/1‑100 (MQ=255)
                              aGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCt                                                                    <  2:733752/100‑1 (MQ=255)
                              aGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGAAGTGCCGCACTAGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCt                                                                    >  2:1010444/1‑100 (MQ=255)
                                cGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGc                                                                  <  2:1014120/100‑1 (MQ=255)
                                cGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGc                                                                  >  1:3411046/1‑100 (MQ=255)
                                 gACTGGAAAACCAGCTTTCGCCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTTAGCGGCTAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCt                                                                    <  2:487827/97‑1 (MQ=255)
                                   tccgatctccAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCttt                                                               <  2:2019915/92‑1 (MQ=255)
                                   cTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCttt                                                               >  1:2438490/1‑100 (MQ=255)
                                   cTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCttt                                                               <  1:1576742/100‑1 (MQ=255)
                                     ggAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTg                                                             >  1:3627224/1‑100 (MQ=255)
                                     ggAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTg                                                             <  1:2771657/100‑1 (MQ=255)
                                     ggAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTg                                                             >  1:1127/1‑100 (MQ=255)
                                     ggAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTg                                                             >  2:1207038/1‑100 (MQ=255)
                                       aaaaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTTGTTAGCGGCGAAGGGATTACTTTCTATGCTGGACTTGTTGTTTCTGCTTTTGcc                                                           <  2:2547557/100‑1 (MQ=255)
                                       aaaaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGcc                                                           <  1:4076749/100‑1 (MQ=255)
                                       aaaaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGcc                                                           <  2:153113/100‑1 (MQ=255)
                                       aaaaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGcc                                                           <  1:2576440/100‑1 (MQ=255)
                                       aaaaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGcc                                                           <  1:1580501/100‑1 (MQ=255)
                                       aaaaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGcc                                                           <  2:377290/100‑1 (MQ=255)
                                       aaaaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGcc                                                           <  2:3385188/100‑1 (MQ=255)
                                        aaaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCt                                                          >  1:1564521/1‑100 (MQ=255)
                                        aaaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCt                                                          >  1:1031798/1‑100 (MQ=255)
                                          aCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGt                                                        <  1:1066251/100‑1 (MQ=255)
                                          aCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGt                                                        >  1:1475744/1‑100 (MQ=255)
                                           ccAGCTTTCACCCGCGACTGACGTGGCTGTAGTGTCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTa                                                       <  1:2557883/100‑1 (MQ=255)
                                           ccAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTa                                                       >  2:1595746/1‑100 (MQ=255)
                                           ccAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTa                                                       >  1:4073669/1‑100 (MQ=255)
                                           ccAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTa                                                       >  2:4043418/1‑100 (MQ=255)
                                           ccAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTa                                                       <  1:3808375/100‑1 (MQ=255)
                                                 ttCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCt                                                 <  2:3205403/100‑1 (MQ=255)
                                                    aCCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGcc                                              <  2:2859649/100‑1 (MQ=255)
                                                    aCCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGcc                                              <  1:142939/100‑1 (MQ=255)
                                                    aCCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGcc                                              <  2:406916/100‑1 (MQ=255)
                                                    aCCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGcc                                              <  2:356214/100‑1 (MQ=255)
                                                     cccGCGACTGACGTGTCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGcc                                              <  2:793266/99‑1 (MQ=255)
                                                     cccGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATCACTTTCTATGCTGGAGCTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCt                                             >  2:3063124/1‑100 (MQ=255)
                                                     cccGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGcc                                              <  1:2352633/99‑1 (MQ=255)
                                                        gcgACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATg                                          <  2:3818167/100‑1 (MQ=255)
                                                          gACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGc                                        <  1:1864309/100‑1 (MQ=255)
                                                          ctcTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGc                                        <  2:3011757/98‑1 (MQ=255)
                                                            cTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGcag                                      >  1:3011775/1‑98 (MQ=255)
                                                            cTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTg                                      <  2:3669960/100‑1 (MQ=255)
                                                             tGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTgg                                     <  2:1745190/100‑1 (MQ=255)
                                                             tGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTgg                                     <  1:1815955/100‑1 (MQ=255)
                                                                cGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGtat                                  <  2:3359837/100‑1 (MQ=255)
                                                                 gTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGtata                                 <  1:956969/100‑1 (MQ=255)
                                                                  tGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTatat                                <  1:3068663/100‑1 (MQ=255)
                                                                    gCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATgg                              <  2:3627197/100‑1 (MQ=255)
                                                                    gCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATgg                              <  1:532861/100‑1 (MQ=255)
                                                                      tGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGc                            >  1:2580878/1‑100 (MQ=255)
                                                                       gTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGcc                           <  1:1005643/100‑1 (MQ=255)
                                                                          gTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCa                        >  2:471716/1‑100 (MQ=255)
                                                                           tGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAg                       <  1:667409/100‑1 (MQ=255)
                                                                           tGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAg                       >  1:3900184/1‑100 (MQ=255)
                                                                            gCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGa                      >  1:3780606/1‑100 (MQ=255)
                                                                             ccGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGaa                     <  1:407559/100‑1 (MQ=255)
                                                                              cGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAg                    <  1:2332014/100‑1 (MQ=255)
                                                                              cGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAg                    <  2:748764/100‑1 (MQ=255)
                                                                               gCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTCTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAg                    >  2:255196/1‑99 (MQ=255)
                                                                                  cTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTgcg                >  1:859796/1‑100 (MQ=255)
                                                                                   tCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTgcgc               <  1:490442/100‑1 (MQ=255)
                                                                                    cGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGca              >  2:3907290/1‑100 (MQ=255)
                                                                                     gTCAGCGGCGAAGGGCTAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGca              <  2:2438531/99‑1 (MQ=255)
                                                                                     gTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGcaa             <  2:563162/100‑1 (MQ=255)
                                                                                     gTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGcaa             >  2:2321013/1‑100 (MQ=255)
                                                                                       cAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCAATGGCTGGTATATGGGCCGCAGAAATGCGcaaca           >  2:2107450/1‑100 (MQ=255)
                                                                                          cggcgAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCaacaaaa        >  2:448587/1‑100 (MQ=255)
                                                                                           ggcgAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCaacaaaac       <  1:1632362/100‑1 (MQ=255)
                                                                                             cgAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCaacaaaaca      <  1:233106/99‑1 (MQ=255)
                                                                                             cgAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTACCTGTAGCCGCTGCCTAGGGCTGGTATATGGGCCGCAGAAGAGCCCaacaaaccca     >  2:2261523/1‑96 (MQ=255)
                                                                                               aaGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAAcaagc   <  1:3607459/100‑1 (MQ=255)
                                                                                               aaGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAAcaag    <  1:1840430/99‑1 (MQ=255)
                                                                                                aGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAAcaagca  <  2:1451974/100‑1 (MQ=255)
                                                                                                   |                                                                                                
TTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAACAAGCA  >  NC_000913/1340448‑1340643

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