Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,390,512 A→G 5.8% intergenic (+63/+71) aaeR → / ← tldD transcriptional regulator for aaeXAB operon/putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,390,5120AG 5.8% 333.0 / 8.2 103intergenic (+63/+71)aaeR/tldDtranscriptional regulator for aaeXAB operon/putative peptidase
Reads supporting (aligned to +/- strand):  ref base A (49/48);  new base G (4/2);  total (53/50)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.92e-01

Only 100 of 108 total aligned reads displayed.

AGGTCGTGATCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGC  >  NC_000913/3390364‑3390655
                                                                                                                                                    |                                                                                                                                               
aGGTCGTGATCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCa                                                                                                                                                 >  1:1105248/1‑149 (MQ=255)
   tCGTGATCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCAccc                                                                                                                                              <  1:759149/149‑1 (MQ=255)
        aTCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCTACGGGTTCCCCTCACCCTAACCCTCACCCCAAAg                                                                                                                                         >  1:1319749/1‑149 (MQ=255)
        aTCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAg                                                                                                                                         >  2:113377/1‑149 (MQ=255)
         tCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAgg                                                                                                                                        <  2:1105248/149‑1 (MQ=255)
            aCTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAgg                                                                                                                                        >  3:1122594/1‑146 (MQ=255)
            aCTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAgg                                                                                                                                        <  4:1122594/146‑1 (MQ=255)
             cTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCAc                                                                                                                                                <  1:1436352/137‑1 (MQ=255)
             cTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCAc                                                                                                                                                >  2:1436352/1‑137 (MQ=255)
              tCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCTACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGa                                                                                                                                   <  2:1319749/149‑1 (MQ=255)
               cGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAg                                                                                                                                  <  4:108758/149‑1 (MQ=255)
                gCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAgg                                                                                                                                 >  3:798033/1‑149 (MQ=255)
                 cTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAggg                                                                                                                                <  1:113377/149‑1 (MQ=255)
                    aCGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAg                                                                                                                                  >  3:1146906/1‑144 (MQ=255)
                    aCGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAg                                                                                                                                  <  4:1146906/144‑1 (MQ=255)
                      ggATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAgg                                                                                                                                        <  4:462663/136‑1 (MQ=255)
                      ggATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAgg                                                                                                                                        >  3:462663/1‑136 (MQ=255)
                       gATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGTGAgg                                                                                                                                 >  1:1391719/1‑142 (MQ=255)
                       gATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCaaa                                                                                                                                          <  3:677231/133‑1 (MQ=255)
                       gATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCaaa                                                                                                                                          >  4:677231/1‑133 (MQ=255)
                       gATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGc                                                                                                                                     >  2:324427/1‑138 (MQ=255)
                       gATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGt                                                                                                                          >  2:208585/1‑149 (MQ=255)
                       gATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGt                                                                                                                          <  1:208585/149‑1 (MQ=255)
                       gATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTCTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGt                                                                                                                          >  1:127495/1‑149 (MQ=255)
                         ttattTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCaaa                                                                                                                                          <  4:1218307/131‑1 (MQ=255)
                         ttattTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCaaa                                                                                                                                          >  3:1218307/1‑131 (MQ=255)
                                 ttgAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGt                                                                                                                >  4:160628/1‑149 (MQ=255)
                                   gAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAggg                                                                                                                                       <  2:781326/124‑1 (MQ=255)
                                   gAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAggg                                                                                                                                       >  1:781326/1‑124 (MQ=255)
                                       tCGGTAAATAGTTTCAGGAGATGCACGGGCGCGGGAAAGAGGAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAg                                                                                                          >  2:217929/1‑149 (MQ=255)
                                            aaaTTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGattaa                                                                                                     <  4:798033/149‑1 (MQ=255)
                                              aTTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCGccccttt                                                                                                                                          <  3:879057/110‑4 (MQ=255)
                                              aTTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCGccccttt                                                                                                                                          >  4:879057/1‑107 (MQ=255)
                                              aTTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAg                                                                                                                                         <  4:392922/111‑1 (MQ=255)
                                               ttGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGTGAgggtt                                                                                                                              >  4:1362845/1‑119 (MQ=255)
                                                  tttCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAg                                                                                                                                         >  4:1456118/1‑107 (MQ=255)
                                                  tttCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAg                                                                                                                                         <  3:1456118/107‑1 (MQ=255)
                                                   ttCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTCTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCt                                                                                              <  2:127495/149‑1 (MQ=255)
                                                     cAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGAACGTAGAAAGATTAATTATCCTTc                                                                                            <  1:217929/149‑1 (MQ=255)
                                                      aGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCAc                                                                                                                  >  3:242584/1‑126 (MQ=255)
                                                               aCGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGt                                                                                  <  3:160628/149‑1 (MQ=255)
                                                                    cgcgGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAat                                                                             >  1:28401/1‑149 (MQ=255)
                                                                        ggAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGaaaa                                                                                      >  4:1004633/1‑136 (MQ=255)
                                                                        ggAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGaaaa                                                                                      <  3:1004633/136‑1 (MQ=255)
                                                                                aaGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGaaaa                                                                                      >  1:800772/1‑128 (MQ=255)
                                                                                aaGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGaaaa                                                                                      <  2:800772/128‑1 (MQ=255)
                                                                                   tAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGattaat                                                                                                    <  3:1373618/111‑1 (MQ=255)
                                                                                   tAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGattaat                                                                                                    >  4:1373618/1‑111 (MQ=255)
                                                                                       tCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAAttatta                                                                        <  2:972993/135‑1 (MQ=255)
                                                                                       tCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAAttatta                                                                        >  1:972993/1‑135 (MQ=255)
                                                                                         atatTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGa                                                                                                                                   <  4:510665/74‑1 (MQ=255)
                                                                                         atatTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGa                                                                                                                                   >  3:510665/1‑74 (MQ=255)
                                                                                         atatTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGaaaa                                                                                      <  1:552121/119‑1 (MQ=255)
                                                                                         atatTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGaaaa                                                                                      >  2:552121/1‑119 (MQ=255)
                                                                                          tatTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAg                                                                                                                                  <  2:393971/74‑1 (MQ=255)
                                                                                          tatTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAg                                                                                                                                  >  1:393971/1‑74 (MQ=255)
                                                                                             tGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAg                                                                                                                                         >  4:523869/1‑64 (MQ=255)
                                                                                             tGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAg                                                                                                                                         <  3:523869/64‑1 (MQ=255)
                                                                                                    ttACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGattaa                                                                                                     <  1:466401/93‑1 (MQ=255)
                                                                                                    ttACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGattaa                                                                                                     >  2:466401/1‑93 (MQ=255)
                                                                                                    ttACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACa                                                        <  4:1313403/138‑1 (MQ=255)
                                                                                                    ttACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACa                                                        >  3:1313403/1‑138 (MQ=255)
                                                                                                        gTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCt                                                                                              <  1:1144234/96‑1 (MQ=255)
                                                                                                        gTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCt                                                                                              >  2:1144234/1‑96 (MQ=255)
                                                                                                            tACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCa                                     >  1:376149/1‑149 (MQ=255)
                                                                                                               aaaCCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGatt                                                                                                       >  2:1429292/1‑80 (MQ=255)
                                                                                                               aaaCCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGatt                                                                                                       <  1:1429292/80‑1 (MQ=255)
                                                                                                               aaaCCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATGGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACg                                  >  3:513424/1‑149 (MQ=255)
                                                                                                                  ccTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTtg                               >  1:46354/1‑149 (MQ=255)
                                                                                                                         cAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAgg                                                  <  2:629259/123‑1 (MQ=255)
                                                                                                                         cAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAgg                                                  >  1:629259/1‑123 (MQ=255)
                                                                                                                          aaCGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGccc                       <  2:46354/149‑1 (MQ=255)
                                                                                                                           aCGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGa                                           <  2:613684/128‑1 (MQ=255)
                                                                                                                           aCGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGa                                           >  1:613684/1‑128 (MQ=255)
                                                                                                                               gTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTGTCCTTCTGAAAATa                                                                                    <  3:106602/83‑1 (MQ=255)
                                                                                                                               gTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATa                                                                                    >  4:106602/1‑83 (MQ=255)
                                                                                                                               gTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGAGCACGTAGAAAGTTTAATTATCCTTCTGAAAATAGTGAAATTATTTCGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGcc                  >  4:1076855/1‑149 (MQ=255)
                                                                                                                                   cccTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAAcc              >  2:687846/1‑149 (MQ=255)
                                                                                                                                    ccTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCGCGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCa                                                                    >  3:926307/1‑94 (MQ=255)
                                                                                                                                    ccTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCa                                                                    <  4:926307/94‑1 (MQ=255)
                                                                                                                                    ccTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTtggctggc                         <  1:990187/137‑1 (MQ=255)
                                                                                                                                    ccTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTtggctggc                         >  2:990187/1‑137 (MQ=255)
                                                                                                                                          ccTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACa                                                        <  2:1090101/100‑1 (MQ=255)
                                                                                                                                          ccTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACa                                                        >  1:1090101/1‑100 (MQ=255)
                                                                                                                                          ccTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAg                                                   >  3:702334/1‑105 (MQ=255)
                                                                                                                                          ccTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAg                                                   <  2:1306272/105‑1 (MQ=255)
                                                                                                                                          ccTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAg                                                   >  1:1306272/1‑105 (MQ=255)
                                                                                                                                          ccTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAg                                                   <  4:702334/105‑1 (MQ=255)
                                                                                                                                           cTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTtg                               <  2:728899/124‑1 (MQ=255)
                                                                                                                                           cTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTtg                               >  1:728899/1‑124 (MQ=255)
                                                                                                                                            tAACCCTCGCCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATc                                             >  4:1151281/1‑109 (MQ=255)
                                                                                                                                            tAACCCTCGCCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATc                                             <  3:1151281/109‑1 (MQ=255)
                                                                                                                                                ccTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAg                                                                   >  3:415793/1‑83 (MQ=255)
                                                                                                                                                ccTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAg                                                                   <  4:415793/83‑1 (MQ=255)
                                                                                                                                                 cTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATc                                             >  2:703374/1‑104 (MQ=255)
                                                                                                                                                 cTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATc                                             <  1:703374/104‑1 (MQ=255)
                                                                                                                                                    aCCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGc                                                             >  4:455889/1‑85 (MQ=255)
                                                                                                                                                    aCCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGc                                                             <  3:455889/85‑1 (MQ=255)
                                                                                                                                                    aCCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGc  <  4:516289/144‑1 (MQ=255)
                                                                                                                                                    aCCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGc  >  3:516289/1‑144 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
AGGTCGTGATCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGC  >  NC_000913/3390364‑3390655

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