Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,373,281 T→A 5.5% intergenic (+119/+84) 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,2810TA 5.5% 231.7 / 2.4 73intergenic (+119/+84)yjeI/yjeJDUF4156 domain‑containing lipoprotein YjeI/protein YjeJ
Reads supporting (aligned to +/- strand):  ref base T (34/35);  new base A (2/2);  total (36/37)
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.75e-01

TATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAA  >  NC_000913/4373132‑4373427
                                                                                                                                                     |                                                                                                                                                  
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                                                                                                                       <  1:538465/150‑1 (MQ=255)
                             gATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCAc                                                                                                                       >  1:387226/1‑150 (MQ=255)
                             gATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCAc                                                                                                                       >  2:1623530/1‑150 (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                                                                  >  2:1429192/1‑150 (MQ=255)
                                                                                  tATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATgag                                                                  >  1:1297048/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                                                                                     <  2:1571223/84‑1 (MQ=255)
                                                                                                                                 tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGcc                                                                                     >  1:1571223/1‑84 (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:829087/1‑150 (MQ=255)
                                                                                                                                              aGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAgcgcgc      >  2:104583/1‑150 (MQ=255)
                                                                                                                                                 cGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTa   <  1:157561/150‑1 (MQ=255)
                                                                                                                                                 cGTTTACGCCGCATCCGGCATGACCAACGATCACGTTAACAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTa   >  2:743255/1‑150 (MQ=255)
                                                                                                                                                  gTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTaa  <  1:1053881/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAA  >  NC_000913/4373132‑4373427

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: