Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,466,123 4 bp→TGCT pseudogene (732‑735/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,1230CT98.4% 143.6 / ‑4.3 61pseudogene (732/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (28/32);  total (28/33)
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.98e-01
*NC_0009131,466,1240AG98.4% 142.3 / ‑4.4 61pseudogene (733/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (28/32);  total (28/33)
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.98e-01
*NC_0009131,466,1250AC98.4% 144.7 / ‑4.3 61pseudogene (734/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (28/32);  total (28/33)
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.91e-01
*NC_0009131,466,1260CT98.4% 140.6 / ‑4.3 61pseudogene (735/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (28/32);  total (28/33)
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.86e-01

TGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGT  >  NC_000913/1465979‑1466201
                                                                                                                                                ||||                                                                           
tGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                                                           <  1:1457914/150‑1 (MQ=255)
        aaGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACgg                                                                  >  2:541784/1‑151 (MQ=255)
        aaGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACgg                                                                  >  2:950636/1‑151 (MQ=255)
         aGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACggg                                                                 >  1:13853/1‑151 (MQ=255)
           gggTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGaa                                                               <  2:1198200/151‑1 (MQ=255)
           gggTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGa                                                                >  1:831384/1‑150 (MQ=255)
            ggTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGa                                                                >  1:1019128/1‑149 (MQ=255)
             gTGAGCGCAGATGGCGCAACGATGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGaagg                                                             <  2:207825/151‑3 (MQ=25)
             gTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGaagg                                                             <  1:274863/151‑3 (MQ=25)
               gAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTg                                                           <  2:1195854/151‑1 (MQ=25)
                 gcgcAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAt                                                         <  2:1375379/151‑1 (MQ=25)
                   gcAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTc                                                       >  1:1081654/1‑151 (MQ=21)
                   gcAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTc                                                       <  1:896362/151‑1 (MQ=21)
                    cAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCa                                                      <  2:890565/151‑1 (MQ=21)
                    cAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCa                                                      <  2:457981/151‑1 (MQ=21)
                     agttGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTc                                                       >  1:1227052/4‑149 (MQ=21)
                         ggCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGAc                                                  <  1:1102287/150‑1 (MQ=21)
                           cgcAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGc                                               <  2:1049795/151‑1 (MQ=21)
                            gcAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGAGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAAAAACGGGAAGGTGATTCAGGTCGGCg                                              >  1:488998/1‑151 (MQ=21)
                             cAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCg                                              >  1:695944/1‑150 (MQ=21)
                             cAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGa                                             <  2:458779/151‑1 (MQ=21)
                             cAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGa                                             <  2:5267/151‑1 (MQ=21)
                              aaCGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGAt                                            >  2:933414/1‑151 (MQ=18)
                              aaCGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGAt                                            >  2:1363351/1‑151 (MQ=18)
                              aaCGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGAt                                            <  2:1329948/151‑1 (MQ=18)
                                 gAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTg                                         <  1:839475/151‑1 (MQ=18)
                                  aaGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTgg                                        <  2:1139218/151‑1 (MQ=18)
                                    gTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                                      <  1:446382/151‑1 (MQ=18)
                                        cAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTc                                   >  2:750428/1‑150 (MQ=18)
                                        cAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCa                                  <  2:916323/151‑1 (MQ=18)
                                         aaTATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTc                                   >  2:810935/1‑149 (MQ=18)
                                          atatAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgc                                <  1:451810/151‑1 (MQ=17)
                                           tataATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcg                               >  2:730547/1‑151 (MQ=17)
                                           tataATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcg                               >  1:975288/1‑151 (MQ=17)
                                            ataATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcg                               <  2:1461530/150‑1 (MQ=17)
                                             taATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggc                             >  1:417005/1‑151 (MQ=17)
                                             taATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggc                             >  2:978190/1‑151 (MQ=17)
                                              aaTCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcg                            <  1:908282/151‑1 (MQ=17)
                                               aTCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcgg                           >  1:939871/1‑151 (MQ=17)
                                               aTCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcgg                           >  2:48723/1‑151 (MQ=17)
                                               aTCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcgg                           <  2:465544/151‑1 (MQ=17)
                                                 cacaACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcggcg                         <  1:1271972/151‑1 (MQ=17)
                                                 cacaACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcggcg                         >  1:262001/1‑151 (MQ=17)
                                                    aaCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTc                      <  2:1266652/151‑1 (MQ=17)
                                                     aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGCTTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTc                      >  1:314040/1‑150 (MQ=17)
                                                     aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCa                     <  1:321327/151‑1 (MQ=14)
                                                       ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg                   <  2:1237300/151‑1 (MQ=14)
                                                        gTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACgg                  <  1:80481/151‑1 (MQ=14)
                                                           aGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAt               <  2:589503/151‑1 (MQ=14)
                                                               tGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGAGACggga                >  2:857538/1‑144 (MQ=14)
                                                                gTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGat           <  1:35044/150‑1 (MQ=14)
                                                                  tattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAGCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTAGATGTCAGCGGCGGCGGTCACGGGATTGat           >  2:413771/1‑148 (MQ=11)
                                                                  tattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGatat         >  2:465645/1‑150 (MQ=11)
                                                                   attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATc        <  1:759576/150‑1 (MQ=11)
                                                                    ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAc      >  1:1000217/1‑151 (MQ=11)
                                                                      aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAcc     >  2:114709/1‑150 (MQ=11)
                                                                      aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCg    <  1:1355095/151‑1 (MQ=11)
                                                                       cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg   <  2:416949/151‑1 (MQ=11)
                                                                       cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg   <  1:1247174/151‑1 (MQ=11)
                                                                       cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg   >  1:584863/1‑151 (MQ=11)
                                                                             aTAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGGGCTTGAGGCCGGCGTTCTGGATGTCAACGGGAGCGGTCACGGTATTGatagcctgggc  >  1:629901/1‑138 (MQ=11)
                                                                                                                                                ||||                                                                           
TGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGT  >  NC_000913/1465979‑1466201

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