Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,597,272 T→A 5.6% pseudogene (815/5422 nt) yneO ← AIDA‑I family autotransporter YneO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,597,2720TA 5.6% 284.0 / 3.3 88pseudogene (815/5422 nt)yneOAIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base T (43/40);  new base A (2/3);  total (45/43)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCCCACCAATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCC  >  NC_000913/1597125‑1597420
                                                                                                                                                   |                                                                                                                                                    
aCGCCCACCAATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGGCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTCAGGTCAGAGTCTTGATCCTGCACATTCAGGGTTCCGTTTCCGGTGGCAttt                                                                                                                                                    <  2:986687/150‑2 (MQ=255)
aCGCCCACCAATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTg                                                                                                                                                    >  1:285569/1‑150 (MQ=255)
 agCCCACCAATGAATATTCTTTCTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGc                                                                                                                                                   >  2:1113562/2‑150 (MQ=255)
  gCCCACCAATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGcc                                                                                                                                                  >  2:154797/1‑150 (MQ=255)
  gCCCACCAATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGcc                                                                                                                                                  <  1:602233/150‑1 (MQ=255)
   cccACCAATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCa                                                                                                                                                 >  2:1147218/1‑150 (MQ=255)
     caccaATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCTAAATTTCCATTATATAGTCTGCGTAACTTTATGACTGAGTCTTTATCCTGCACATTCTGGGTTCCGATACCGGTGGCATTGCCAcc                                                                                                                                               <  2:552008/150‑1 (MQ=255)
     caccaATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCAcc                                                                                                                                               >  1:800954/1‑150 (MQ=255)
          aTGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTat                                                                                                                                          <  2:885731/150‑2 (MQ=255)
           tGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAt                                                                                                                                         >  2:1050913/1‑150 (MQ=255)
            gAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATc                                                                                                                                        >  1:682873/1‑150 (MQ=255)
            gAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATc                                                                                                                                        >  2:973270/1‑150 (MQ=255)
              atatTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTatcgga                                                                                                                                      <  1:308321/150‑6 (MQ=255)
                atTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCATGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAt                                                                                                                                    >  1:869083/1‑150 (MQ=255)
                      tGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGAGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGTTGGCATTGCCACCGATAATCGTATTTTCAg                                                                                                                              >  1:436416/1‑150 (MQ=255)
                        ttattaATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATTACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGcc                                                                                                                            >  2:54405/1‑150 (MQ=255)
                        ttattaATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGcc                                                                                                                            <  1:1145026/150‑1 (MQ=255)
                         tattaATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCAGATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCg                                                                                                                           >  1:112914/1‑150 (MQ=255)
                         tattaATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCg                                                                                                                           <  1:787910/150‑1 (MQ=255)
                         tattaATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCg                                                                                                                           >  2:134429/1‑150 (MQ=255)
                          attaATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGt                                                                                                                          >  2:1396904/1‑150 (MQ=255)
                          attaATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGt                                                                                                                          >  2:213120/1‑150 (MQ=255)
                          attaATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGt                                                                                                                          <  2:58413/150‑1 (MQ=255)
                            taATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATTACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaa                                                                                                                        <  1:54405/150‑1 (MQ=255)
                             aaTCAGTCCATTATTTGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGt                                                                                                                          >  2:133089/1‑147 (MQ=255)
                             aaTCAGTCCATTATTTGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGt                                                                                                                          <  1:133089/147‑1 (MQ=255)
                             aaTCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaa                                                                                                                       >  1:1092621/1‑149 (MQ=255)
                                  gTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaacc                                                                                                                  >  2:1077523/1‑150 (MQ=255)
                                    ccATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGAGAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTACCCAAcccc                                                                                                                >  2:818963/1‑150 (MQ=255)
                                     cATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAccccc                                                                                                               <  2:1574600/150‑1 (MQ=255)
                                     cATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAccccc                                                                                                               >  1:1318764/1‑150 (MQ=255)
                                           tGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCCCCCCCCGCTcgc                                                                                                         >  1:622265/1‑150 (MQ=255)
                                           tGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTcgc                                                                                                         >  1:1507289/1‑150 (MQ=255)
                                            ggAAATATTGACTGTGCCATTACCGAAATATCCATTATAGCGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACTTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTcgcg                                                                                                        <  2:1116123/150‑1 (MQ=255)
                                            ggAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTcgcg                                                                                                        >  2:706420/1‑150 (MQ=255)
                                             gAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGa                                                                                                       >  2:836778/1‑150 (MQ=255)
                                                atatTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATa                                                                                                    <  1:927257/150‑1 (MQ=255)
                                                  atTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGt                                                                                                  >  2:1208109/1‑150 (MQ=255)
                                                    tGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCg                                                                                                <  2:847851/150‑1 (MQ=255)
                                                     gACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGc                                                                                               >  2:1208000/1‑150 (MQ=255)
                                                     gACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGc                                                                                               >  1:1104519/1‑150 (MQ=255)
                                                      aCTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGcc                                                                                              <  2:1501001/150‑1 (MQ=255)
                                                              aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGt                                                                                      >  1:1423101/1‑150 (MQ=255)
                                                              aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGt                                                                                      <  1:424729/150‑1 (MQ=255)
                                                                  ccGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCt                                                                                  <  2:1020731/150‑1 (MQ=255)
                                                                    gAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                                                <  1:836778/150‑1 (MQ=255)
                                                                    gAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                                                >  2:836636/1‑150 (MQ=255)
                                                                       atatCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                                             <  1:1147218/150‑1 (MQ=255)
                                                                        tatCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                                                                            <  1:102469/150‑1 (MQ=255)
                                                                        tatCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                                                                            <  2:1281446/150‑1 (MQ=255)
                                                                        tatCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                                                                            <  1:875873/150‑1 (MQ=255)
                                                                        tatCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                                                                            <  1:60939/150‑1 (MQ=255)
                                                                          tCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAAACAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAcc                                                                          <  2:1092621/150‑1 (MQ=255)
                                                                           ccATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCACCCCCCCCTCGCGAATAGTCGCGTCCCCAGTTCCTTGATTACCa                                                                         >  2:932783/1‑150 (MQ=255)
                                                                           ccATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCa                                                                         >  2:1224470/1‑150 (MQ=255)
                                                                              ttATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                                                                      <  2:1081588/150‑1 (MQ=255)
                                                                               tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt                                                                     <  1:1519002/150‑1 (MQ=255)
                                                                                  agagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                                                                  <  1:134429/150‑1 (MQ=255)
                                                                                   gagGCTGCGTACCGTGATGACTGACTCTTGATCCTGCACATTCAGGGTTGCGCTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAGCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                                                 <  2:635106/150‑1 (MQ=255)
                                                                                    agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                                                <  1:1371783/150‑1 (MQ=255)
                                                                                            tACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGa                                                        >  2:15059/1‑150 (MQ=255)
                                                                                             aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAg                                                       <  1:1050913/150‑1 (MQ=255)
                                                                                             aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAg                                                       <  1:1113562/150‑1 (MQ=255)
                                                                                              ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                                      <  2:285569/150‑1 (MQ=255)
                                                                                                gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCACTTCCTTGATTACCAATTTGAAATTCAATTGAGGaa                                                    >  2:310229/1‑150 (MQ=255)
                                                                                                  gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGACACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATc                                                  >  2:743586/1‑150 (MQ=255)
                                                                                                       cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTGCCTTGATTACCAATTTGAAATTCAATTGAGGAATCattat                                             >  1:203842/1‑150 (MQ=255)
                                                                                                         gagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAttt                                           <  2:1358585/150‑1 (MQ=255)
                                                                                                             cTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGGCGCCGCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAAGCATTATTTTtta                                       <  1:675761/150‑1 (MQ=255)
                                                                                                             cTTGATCCTGCACAATCAGGGGTCCGATACCGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                                       <  1:356174/150‑1 (MQ=255)
                                                                                                              ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                                      <  1:1077523/150‑1 (MQ=255)
                                                                                                              ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                                      <  1:973270/150‑1 (MQ=255)
                                                                                                              ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCAACGCGATTAGTCGCGGCCCCCGTTCGTGGATTACCAATTTAAAATTAAATTGAGGAAGAATTATTTTttat                                      >  2:234119/1‑150 (MQ=255)
                                                                                                                      gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACt                              >  2:954293/1‑150 (MQ=255)
                                                                                                                        acaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCg                            <  1:836636/150‑1 (MQ=255)
                                                                                                                         caTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGc                           >  2:233648/1‑150 (MQ=255)
                                                                                                                         caTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGc                           <  1:1208000/150‑1 (MQ=255)
                                                                                                                            tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCCACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTGCTTGATTACCAATTTGAAATTCAATTGAGGAGTCATTATTTTTTATTAGCCACTCGcccc                        >  1:1010106/1‑148 (MQ=255)
                                                                                                                            tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCACCCCCCCCTCGCGCATAGTCGCCGCCCCCGGTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGcccc                        >  1:199125/1‑148 (MQ=255)
                                                                                                                            tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATAATTATTTTTTATTAGCCCCTCGacac                        >  2:6385/1‑150 (MQ=255)
                                                                                                                                    ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCa                <  2:1423101/150‑1 (MQ=255)
                                                                                                                                     gattatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTGCTTGAGTACCAATTTGAAATTCAATTGAGGAATCAGTATTTTTTATTGGCCACTCGCCGCCCTTTTCaa               >  2:1577344/7‑150 (MQ=255)
                                                                                                                                      attatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTGCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAAc              >  2:738001/6‑150 (MQ=255)
                                                                                                                                       aTACCGTTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACg             <  2:436416/150‑1 (MQ=255)
                                                                                                                                            ggTGGCATTGCCACCGATAATCGTATTTTAAGCCGTAACCAAACCCCCCTCGCGAATAGTCGCATCCCCAATTACTTGATTACAAAATTTAAATTAAATTGAGCAATCATAATTTTTTATTAGCCACTCCCCACCCTTTTCAAccaaaaa        >  2:823340/1‑144 (MQ=255)
                                                                                                                                                gCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTgg    <  1:154797/150‑1 (MQ=255)
                                                                                                                                                gCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTgg    <  1:823340/150‑1 (MQ=255)
                                                                                                                                                 cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGGTCCTTGATTAGCAATTTCAAATTCAATTGAGGAAACATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGGCAACCTGGc   >  1:594239/1‑150 (MQ=255)
                                                                                                                                                  aTTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGcc  <  1:706420/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACGCCCACCAATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCC  >  NC_000913/1597125‑1597420

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