Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,681,969 C→G 5.2% Q259H (CAG→CAC yphG ← DUF5107 domain‑containing protein YphG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,681,9690CG 5.2% 326.2 / 5.8 96Q259H (CAG→CACyphGDUF5107 domain‑containing protein YphG
Reads supporting (aligned to +/- strand):  ref base C (45/46);  new base G (3/2);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

CCGTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCG  >  NC_000913/2681820‑2682117
                                                                                                                                                     |                                                                                                                                                     
ccGTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTAc                                                                                                                                                       <  1:1070/150‑1 (MQ=255)
    tATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCg                                                                                                                                                   <  1:244150/150‑1 (MQ=255)
     aTTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGa                                                                                                                                                  >  2:1603309/1‑150 (MQ=255)
       tGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGAtt                                                                                                                                                >  1:770587/1‑150 (MQ=255)
        gTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGAttt                                                                                                                                               >  1:1263085/1‑150 (MQ=255)
          cGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGAttttt                                                                                                                                            >  2:1477997/1‑150 (MQ=255)
          cGGTCAGACTCTTATCCCACGCCTGGACAAAATCACTGTGTCCCCAACTCCACTGTTTTTTACCAGGCGCAATATTGTGGTTGGCAACGTGCACCAAACCGCCATATACATCCTCACAACACGCGCCAACAAAATCGAACTGCGAttttt                                                                                                                                            >  2:1501296/1‑150 (MQ=255)
            gTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCg                                                                                                                                          <  1:1358122/150‑1 (MQ=255)
             tCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGTTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCgg                                                                                                                                         <  1:1099340/150‑1 (MQ=255)
             tCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCgg                                                                                                                                         <  2:1498159/150‑1 (MQ=255)
             tCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCgg                                                                                                                                         <  1:1139174/150‑1 (MQ=255)
             tCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGAAAAATCGCAGTACGAttttgtcgg                                                                                                                                         <  2:274015/150‑6 (MQ=255)
             tCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGAAAAATCGCAGTACGAttttgtcgg                                                                                                                                         <  2:36430/150‑6 (MQ=255)
              cAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGc                                                                                                                                        <  1:1096113/150‑1 (MQ=255)
              cAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGc                                                                                                                                        <  1:87614/150‑1 (MQ=255)
                    ctTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCatata                                                                                                                                  <  2:502156/150‑1 (MQ=255)
                    ctTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCatata                                                                                                                                  <  1:875440/150‑1 (MQ=255)
                    ctTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCatata                                                                                                                                  >  1:1175141/1‑150 (MQ=255)
                        tCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAg                                                                                                                              >  2:975431/1‑150 (MQ=255)
                          ccACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGa                                                                                                                               >  2:607508/1‑147 (MQ=255)
                          ccACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGa                                                                                                                               <  1:607508/147‑1 (MQ=255)
                          ccACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGt                                                                                                                            >  2:200809/1‑150 (MQ=255)
                           cACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGAATTTTCGGCCATATATGAGGtt                                                                                                                           >  1:483190/1‑150 (MQ=255)
                            aCGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTg                                                                                                                          <  2:184850/150‑1 (MQ=255)
                             cGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTgg                                                                                                                         <  2:1253756/150‑1 (MQ=255)
                             cGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTgg                                                                                                                         <  1:307978/150‑1 (MQ=255)
                              gCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGa                                                                                                                        >  1:157915/1‑150 (MQ=255)
                                cTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAAc                                                                                                                      <  2:927892/150‑1 (MQ=255)
                                cTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAAc                                                                                                                      <  1:459908/150‑1 (MQ=255)
                                   gCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAgg                                                                                                                   <  1:1017125/150‑1 (MQ=255)
                                    ccAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGc                                                                                                                  <  2:1052155/150‑1 (MQ=255)
                                    ccAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGc                                                                                                                  <  1:94214/150‑1 (MQ=255)
                                     cAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGca                                                                                                                 >  2:1085401/1‑150 (MQ=255)
                                           aaCTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGTCAACGTTCAGCAAACCTCCATCTTCATCGTTACACCAATCGCCGACAAAATCGTACTGCTATTTTTCGGCCATATATGAGGTTGGAACAGGCACAtttt                                                                                                           <  2:755562/149‑1 (MQ=255)
                                            aCTGTGTCCCCAACTGCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCATCGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACAttttt                                                                                                          <  2:1126553/150‑1 (MQ=255)
                                             cTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTta                                                                                                         <  1:25902/150‑1 (MQ=255)
                                             cTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTta                                                                                                         >  1:738614/1‑150 (MQ=255)
                                                tgtCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAg                                                                                                      <  2:104619/150‑1 (MQ=255)
                                                tgtCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAg                                                                                                      <  2:98710/150‑1 (MQ=255)
                                                 gtacccAACTCCACTGTTTTTTACCTGGCGCAATATTGTGGTTGGCAACGTGCAGCAAACCGCAATCTTCATCGTGACACCACGCTCCGACAAAATCTTACTGAGATTTTTCGGCCATATATGAGGTTGGAACAGTCACATTTTTATAga                                                                                                     >  2:868594/4‑149 (MQ=255)
                                                  tCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTATTATAGCg                                                                                                    >  1:1616252/1‑150 (MQ=255)
                                                    cccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCgag                                                                                                  <  2:58219/150‑1 (MQ=255)
                                                        aCTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAAt                                                                                              >  2:911881/1‑150 (MQ=255)
                                                         cTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATg                                                                                             >  1:1268298/1‑150 (MQ=255)
                                                                   ttttACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAg                                                                                   >  1:330397/1‑150 (MQ=255)
                                                                     ttACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCg                                                                                 >  2:607254/1‑150 (MQ=255)
                                                                      tACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCgg                                                                                >  2:662586/1‑150 (MQ=255)
                                                                       aCCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGa                                                                               >  1:591508/1‑150 (MQ=255)
                                                                       aCCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGAGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCgca                                                                               >  2:1086925/1‑148 (MQ=255)
                                                                        ccTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGag                                                                              <  1:1477997/150‑1 (MQ=255)
                                                                         cTGGGGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGagt                                                                             >  2:1390371/1‑150 (MQ=255)
                                                                          tGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGagta                                                                            >  1:233349/1‑150 (MQ=255)
                                                                           ggCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGagtag                                                                           >  2:347723/1‑150 (MQ=255)
                                                                           ggCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGagtag                                                                           <  1:1297864/150‑1 (MQ=255)
                                                                           ggCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGagtag                                                                           >  2:996364/1‑150 (MQ=255)
                                                                            gcgcAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGagtagt                                                                          >  2:611066/1‑150 (MQ=255)
                                                                             cgcAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTc                                                                         >  1:1519318/1‑150 (MQ=255)
                                                                              gcAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTcc                                                                        <  1:313947/150‑1 (MQ=255)
                                                                                aaTATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCAc                                                                      <  2:770587/150‑1 (MQ=255)
                                                                                    tggtggTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTg                                                                  <  1:1235316/150‑1 (MQ=255)
                                                                                     ggtggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGAGGAGTAGTCCACTTTgt                                                                 >  2:395614/1‑150 (MQ=255)
                                                                                      gtggtTGGCAACTTTCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgta                                                                <  2:467908/150‑1 (MQ=255)
                                                                                           tGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAg                                                           <  2:824380/150‑1 (MQ=255)
                                                                                            ggCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGt                                                          >  2:24783/1‑150 (MQ=255)
                                                                                             gCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTg                                                         >  1:1471901/1‑150 (MQ=255)
                                                                                                aCGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgt                                                                 >  2:1119495/1‑139 (MQ=255)
                                                                                                aCGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgt                                                                 <  1:1119495/139‑1 (MQ=255)
                                                                                                  gTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGt                                                    <  2:1263085/150‑1 (MQ=255)
                                                                                                  gTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGAGGAGTAGTCCACTTTGTAGTAAGTGCCGGt                                                    >  2:994743/1‑150 (MQ=255)
                                                                                                     cagcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGc                                                 <  1:1224230/150‑1 (MQ=255)
                                                                                                     cagcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGc                                                 <  1:1603309/150‑1 (MQ=255)
                                                                                                       gcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGa                                               >  1:957524/1‑150 (MQ=255)
                                                                                                          aaCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgg                                            >  2:997512/1‑150 (MQ=255)
                                                                                                          aaCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgg                                            >  2:1362929/1‑150 (MQ=255)
                                                                                                            ccGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgggg                                          >  1:1293779/1‑150 (MQ=255)
                                                                                                            ccGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgggg                                          >  1:436074/1‑150 (MQ=255)
                                                                                                            ccGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATGTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgggg                                          >  2:1027029/1‑150 (MQ=255)
                                                                                                                cATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAg                                      <  2:1470957/150‑1 (MQ=255)
                                                                                                                  tCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCg                                    <  1:611066/150‑1 (MQ=255)
                                                                                                                   cTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCgg                                   >  1:1417182/1‑150 (MQ=255)
                                                                                                                      cATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGgaga                                >  1:1618441/1‑150 (MQ=255)
                                                                                                                            gACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGccc                          <  2:695576/150‑1 (MQ=255)
                                                                                                                              caccacGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGt                        >  1:952845/1‑150 (MQ=255)
                                                                                                                              caccacGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGt                        <  1:347723/150‑1 (MQ=255)
                                                                                                                                ccacGCGCCGACAAAATCGTACTGCGAATTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGttt                      >  2:1412706/1‑150 (MQ=255)
                                                                                                                                  acGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGCACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGc                    >  1:1114931/1‑150 (MQ=255)
                                                                                                                                   cgcgCCGACAAAATCGTACTGCGATTTTTCGGCCATATATTAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAGTGCAGGTAGAGATGGGGAAAGCGGAGACGGCCCGTTTGCCg                  <  2:970368/150‑1 (MQ=255)
                                                                                                                                   cgcgCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGcc                   <  1:156590/150‑1 (MQ=255)
                                                                                                                                       gccgaaaaaTCGCAGTACGATTTTGTCGGCCATATATGAGGTTTGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGGAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTg                >  2:1426417/6‑150 (MQ=255)
                                                                                                                                       ccGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGa               <  2:1519318/150‑1 (MQ=255)
                                                                                                                                        ccgaaaaaTCGCAGTACGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGAt              >  2:1177435/5‑150 (MQ=255)
                                                                                                                                         gACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATc             <  2:738614/150‑1 (MQ=255)
                                                                                                                                         gACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATc             <  1:975431/150‑1 (MQ=255)
                                                                                                                                           cAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCAGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCaa           >  1:189810/1‑150 (MQ=255)
                                                                                                                                            aaaaTCGCAGTACGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCaaa          >  1:371129/1‑150 (MQ=255)
                                                                                                                                              aaTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAaca        <  1:991231/150‑1 (MQ=255)
                                                                                                                                                    aCTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACAccgccg  <  1:1345420/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCG  >  NC_000913/2681820‑2682117

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: