Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,373,274 A→C 5.3% intergenic (+112/+91) yjeI → / ← yjeJ DUF4156 domain‑containing lipoprotein YjeI/protein YjeJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,373,2740AC 5.3% 236.9 / 4.2 75intergenic (+112/+91)yjeI/yjeJDUF4156 domain‑containing lipoprotein YjeI/protein YjeJ
Reads supporting (aligned to +/- strand):  ref base A (34/37);  new base C (2/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGC  >  NC_000913/4373125‑4373423
                                                                                                                                                     |                                                                                                                                                     
gCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATaa                                                                                                                                                       >  1:140137/1‑150 (MQ=255)
 ccAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAg                                                                                                                                                      >  1:380948/1‑150 (MQ=255)
 ccAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAg                                                                                                                                                      >  1:1629748/1‑150 (MQ=255)
  cAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAgg                                                                                                                                                     >  1:1236495/1‑150 (MQ=255)
  cAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAgg                                                                                                                                                     <  2:542770/150‑1 (MQ=255)
     aTTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGt                                                                                                                                                  <  2:1153477/150‑1 (MQ=255)
     aTTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGt                                                                                                                                                  <  1:1436348/150‑1 (MQ=255)
       tATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGttt                                                                                                                                                <  2:140137/150‑1 (MQ=255)
           ggTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                            >  2:538465/1‑150 (MQ=255)
            gTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGATGCTGTGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGTCCTACATAAGGCATTGCAACTGATTGAATGCTTGATAATTTGTAGGCCGGATAAGGCGTTTActcc                                                                                                                                           <  2:1050521/150‑1 (MQ=39)
            gTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGTGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                           <  2:1392754/150‑1 (MQ=255)
             tCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                          <  2:553090/150‑1 (MQ=255)
                ggTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCAt                                                                                                                                       >  1:476348/1‑150 (MQ=255)
                  tttATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                     <  1:1057749/150‑1 (MQ=255)
                    tttaAGTGCCCGTACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCTCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                   <  2:1509411/148‑1 (MQ=255)
                      taAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                 <  1:1498736/150‑2 (MQ=39)
                       aaGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcct                                                                                                                                <  2:1075713/150‑3 (MQ=39)
                                   tGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCa                                                                                                                    >  2:151970/1‑150 (MQ=255)
                                    gATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCAc                                                                                                                   >  2:1623530/1‑150 (MQ=255)
                                    gATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCAc                                                                                                                   >  1:387226/1‑150 (MQ=255)
                                    gATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCAc                                                                                                                   <  1:538465/150‑1 (MQ=255)
                                    gATGTGAACAAGAGCCTGGGTGACAAACGTATAATTGCATGATGCGATTCGCTTATCAGGCCTACATAATGCATTCCAACTGATTTAATGCTTGAGAAATTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCaa                                                                                                                   >  2:682422/1‑149 (MQ=37)
                                       gtgATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGtt                                                                                                                >  2:466549/1‑150 (MQ=255)
                                         gATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTAt                                                                                                              >  1:1261412/1‑150 (MQ=255)
                                           tCCAGAGGCTGTGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACTCCGCATCCGGCATGACCAACGATCACGTTATCa                                                                                                            <  2:973262/150‑1 (MQ=255)
                                           tCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCa                                                                                                            <  2:1138444/150‑1 (MQ=255)
                                            ccAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAg                                                                                                           <  1:1434818/150‑1 (MQ=255)
                                              agagGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTc                                                                                                         >  2:290100/1‑150 (MQ=255)
                                              agagGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTc                                                                                                         >  1:1230612/1‑150 (MQ=255)
                                                    tGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCa                                                                                                   >  2:1053881/1‑150 (MQ=255)
                                                    tGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCa                                                                                                   <  1:151970/150‑1 (MQ=255)
                                                     gggTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAg                                                                                                  >  2:157561/1‑150 (MQ=255)
                                                       gTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGc                                                                                                <  1:951444/150‑1 (MQ=255)
                                                        tGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCt                                                                                               >  2:358071/1‑150 (MQ=255)
                                                            aaaCGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGAt                                                                                           >  1:1079785/1‑150 (MQ=255)
                                                               cGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGc                                                                                        <  2:387226/150‑1 (MQ=255)
                                                                   taATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCa                                                                                    <  2:201065/150‑1 (MQ=255)
                                                                     aTTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGc                                                                                  >  1:851054/1‑150 (MQ=255)
                                                                        gCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAg                                                                                                  >  1:948011/1‑131 (MQ=255)
                                                                        gCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAg                                                                                                  <  2:948011/131‑1 (MQ=255)
                                                                        gCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGAGCCATCCAGCCtg                                                                               <  2:1386277/150‑2 (MQ=255)
                                                                         ccTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGAGCCATCCAGCCtga                                                                              <  1:13398/150‑3 (MQ=255)
                                                                              tGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTa                                                                         >  1:1021188/1‑150 (MQ=255)
                                                                               gcgcTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAg                                                                        <  1:1033393/150‑1 (MQ=255)
                                                                                cgcttcgcttATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGTCTGGATGGCTCCAGCCTCGCATTAGc                                                                       >  1:1614489/1‑150 (MQ=255)
                                                                                 gcttcgcttATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCa                                                                      <  2:1171277/150‑1 (MQ=255)
                                                                                  cttcgcttATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCaa                                                                     >  1:1406509/1‑150 (MQ=255)
                                                                                     cgcttATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGa                                                                  <  2:865418/150‑1 (MQ=255)
                                                                                         tATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATgag                                                              >  1:1297048/1‑150 (MQ=255)
                                                                                         tATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATgag                                                              >  2:1429192/1‑150 (MQ=255)
                                                                                         tATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATgag                                                              >  1:636557/1‑150 (MQ=255)
                                                                                            cAGGCCTACATAAGGCATTGCAATTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGatt                                                           <  2:96309/150‑1 (MQ=255)
                                                                                            cAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGatt                                                           <  2:1629748/150‑1 (MQ=255)
                                                                                                      tAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTACTTCGCCtg                                                 >  2:1007830/1‑150 (MQ=255)
                                                                                                         ggCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTgtgg                                              <  1:1429192/150‑1 (MQ=255)
                                                                                                              tGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAg                                         <  2:318752/150‑1 (MQ=255)
                                                                                                              tGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAg                                         <  1:321869/150‑1 (MQ=255)
                                                                                                               gCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGt                                        >  2:1241519/1‑150 (MQ=255)
                                                                                                                 aaCTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGttt                                      <  2:1145676/150‑1 (MQ=255)
                                                                                                                     gATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCt                                  <  2:636557/150‑1 (MQ=255)
                                                                                                                       ttGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCtgg                                <  2:1260118/150‑1 (MQ=255)
                                                                                                                         gAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCtggtg                              >  1:319035/1‑150 (MQ=255)
                                                                                                                             gCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATTGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGAtt                          <  2:1297048/150‑1 (MQ=255)
                                                                                                                                tGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGaa                       >  1:556967/1‑150 (MQ=255)
                                                                                                                                tGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGaa                       >  1:802531/1‑150 (MQ=255)
                                                                                                                                    aaTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACg                   <  2:1236495/150‑1 (MQ=255)
                                                                                                                                      tttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGAGGAGATTCTTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTa                 <  2:1569109/150‑1 (MQ=255)
                                                                                                                                        tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGcc                                                                                 >  1:1571223/1‑84 (MQ=255)
                                                                                                                                        tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGcc                                                                                 <  2:1571223/84‑1 (MQ=255)
                                                                                                                                        tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAg               >  2:752057/1‑150 (MQ=255)
                                                                                                                                            tgCCGGATGCGGCGTAAACGCCTTATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGAt           >  1:747481/2‑150 (MQ=255)
                                                                                                                                            tgCCGGATGCGGCGTAAACGCCTTATACGGCATGACCAACGATCACGTTATAAGTCATTTCAGGCTGGATGGCGCCAGCCTCGCCTTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGGGGATTGAACACGTAAGTGAt           >  2:406734/2‑150 (MQ=255)
                                                                                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGAt           >  1:496707/1‑150 (MQ=255)
                                                                                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATg          >  2:1204995/1‑150 (MQ=255)
                                                                                                                                              ccGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGc         <  2:380948/150‑1 (MQ=255)
                                                                                                                                                ggATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAg       <  2:314896/150‑1 (MQ=255)
                                                                                                                                                     aGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAgcgcgc  >  2:104583/1‑150 (MQ=255)
                                                                                                                                                     aGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAgcgcgc  >  2:829087/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGC  >  NC_000913/4373125‑4373423

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