Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,221,238 T→G pseudogene (1638/2574 nt) ycgH → pseudogene; putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,221,2380TG96.5% 317.9 / 1.7 85pseudogene (1638/2574 nt)ycgHpseudogene; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base G (39/43);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GATGAATCTCCAGCAGATGCCTTCAGATAAGGTGATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAATTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAAC  >  NC_000913/1221092‑1221383
                                                                                                                                                  |                                                                                                                                                 
gATGAATCTCCAGCAGATGCCTTCAGATAAGGTGATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTc                                                                                                                                                 >  2:3378/1‑149 (MQ=255)
gATGAATCTCCAGCAGATGCCTTCAGATAAGGTGATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAg                                                                                                                                               <  2:1076663/151‑1 (MQ=255)
 aTGAATCTCCAGCAGATGCCTTCAGATAAGGTGATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGa                                                                                                                                              <  1:33290/151‑1 (MQ=255)
         ccAGCAGATGCCTTCAGATAAGGTGATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACa                                                                                                                                       <  1:1066791/150‑1 (MQ=255)
         ccAGCAGATGCCTTCAGATAAGGTGATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACa                                                                                                                                       <  2:752639/150‑1 (MQ=255)
         ccAGCAGATGCCTTCAGATAAGGTGATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAg                                                                                                                                      <  1:947333/151‑1 (MQ=255)
           agcagATGCCTTCAGATAAGGTGATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGa                                                                                                                                     >  1:135260/1‑150 (MQ=255)
                aTGCCTTCAGATAAGGTGATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGcc                                                                                                                               >  1:387201/1‑151 (MQ=255)
                aTGCCTTCAGATAAGGTGATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGcc                                                                                                                               >  1:225087/1‑151 (MQ=255)
                  gCCTTCAGATAAGGTGATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTa                                                                                                                             <  1:21015/151‑1 (MQ=255)
                     ttCAGATAAGGTGATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTAccc                                                                                                                          >  2:1123009/1‑151 (MQ=255)
                      tCAGATAAGGTGATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTAcc                                                                                                                           <  1:1133452/149‑1 (MQ=255)
                       cAGATAAGGTGATGCAAGGATTGCAGCTGTAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTAcccc                                                                                                                         <  2:728333/150‑1 (MQ=255)
                          aTAAGGTGATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTc                                                                                                                     >  2:742386/1‑151 (MQ=255)
                           tAAGGTGATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTcc                                                                                                                    >  2:1223252/1‑151 (MQ=255)
                            aaGGTGATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCt                                                                                                                   <  2:1013550/151‑1 (MQ=255)
                                 gATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTgcgc                                                                                                               <  2:320060/150‑1 (MQ=255)
                                    gCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGcccc                                                                                                            <  2:1065849/150‑1 (MQ=255)
                                    gCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCg                                                                                                           <  2:448518/151‑1 (MQ=255)
                                     cAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCg                                                                                                           >  1:522010/1‑150 (MQ=255)
                                       aGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAg                                                                                                        <  2:60769/151‑1 (MQ=255)
                                       aGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAg                                                                                                        <  2:296199/151‑1 (MQ=255)
                                        ggATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGc                                                                                                       <  1:325714/151‑1 (MQ=255)
                                          aTTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTg                                                                                                     >  2:138666/1‑151 (MQ=255)
                                           ttGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTgg                                                                                                    >  2:697634/1‑151 (MQ=255)
                                            tGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTgg                                                                                                    >  1:1044084/1‑150 (MQ=255)
                                             gCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGaa                                                                                                  <  2:201787/151‑1 (MQ=255)
                                              cAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAATTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGa                                                                                                   >  2:1269195/1‑149 (MQ=255)
                                               aGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGt                                                                                                >  2:807797/1‑151 (MQ=255)
                                                gCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGt                                                                                                <  2:988924/150‑1 (MQ=255)
                                                 cTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGtt                                                                                               >  1:580824/1‑150 (MQ=255)
                                                     aGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTc                                                                                          >  2:90081/1‑151 (MQ=255)
                                                     aGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTc                                                                                          <  2:1001932/151‑1 (MQ=255)
                                                        tttCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCt                                                                                       >  1:1226653/1‑151 (MQ=255)
                                                          tCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTa                                                                                      >  1:1197118/1‑150 (MQ=255)
                                                           cGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTaa                                                                                     <  1:1244898/150‑1 (MQ=255)
                                                               tatCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAATTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTgg                                                                                >  1:1250332/1‑151 (MQ=255)
                                                                  cAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCa                                                                              <  1:453513/150‑1 (MQ=255)
                                                                  cAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAg                                                                             <  2:1196627/151‑1 (MQ=255)
                                                                      tAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAgctgc                                                                         >  1:860182/1‑151 (MQ=255)
                                                                      tAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAgctgc                                                                         >  1:44128/1‑151 (MQ=255)
                                                                       aaaaCAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAgctgc                                                                         <  1:594698/150‑1 (MQ=255)
                                                                           caacaaGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAaca                                                                    <  2:1103922/151‑1 (MQ=255)
                                                                             acaaGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAacac                                                                   >  1:139008/1‑150 (MQ=255)
                                                                                   cAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCttttt                                                            >  2:1062745/1‑151 (MQ=255)
                                                                                    aaaaaCTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTg                                                           <  1:257146/151‑1 (MQ=255)
                                                                                     aaaaCTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCATCAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGt                                                          <  1:1190513/151‑1 (MQ=255)
                                                                                      aaaCTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGtt                                                         <  1:1158022/151‑1 (MQ=255)
                                                                                       aaCTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTa                                                        <  2:595924/151‑1 (MQ=255)
                                                                                         cTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAATTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATg                                                      >  1:1177666/1‑151 (MQ=255)
                                                                                           ggTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCg                                                    <  2:387211/151‑1 (MQ=255)
                                                                                                  gACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAAc                                             <  1:609547/151‑1 (MQ=255)
                                                                                                   aCAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACg                                            >  2:944027/1‑151 (MQ=255)
                                                                                                   aCAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAAACAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACg                                            <  2:1052265/151‑1 (MQ=255)
                                                                                                    cAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGa                                           <  2:522025/151‑1 (MQ=255)
                                                                                                     aaGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGAc                                          >  1:128686/1‑151 (MQ=255)
                                                                                                      aGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGAcc                                         <  1:999201/151‑1 (MQ=255)
                                                                                                      aGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGAcc                                         <  2:201895/151‑1 (MQ=255)
                                                                                                       gTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCg                                        >  2:998661/1‑151 (MQ=255)
                                                                                                       gTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCg                                        <  2:1260970/151‑1 (MQ=255)
                                                                                                         tATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTg                                                >  1:609954/1‑141 (MQ=255)
                                                                                                         tATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTg                                                <  2:609972/141‑1 (MQ=255)
                                                                                                         tATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGt                                       <  1:104885/150‑1 (MQ=255)
                                                                                                          aTCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGc                                     <  2:1122472/151‑1 (MQ=255)
                                                                                                            cagtcagtGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATTCGTCTGAACGACCGTGCg                                    >  2:348089/1‑150 (MQ=255)
                                                                                                            cagtcagtGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCgg                                   >  2:234943/1‑151 (MQ=255)
                                                                                                              gtcagtGAATGAGGCGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCggg                                  >  1:77479/1‑150 (MQ=255)
                                                                                                                   tGAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAAc                            >  2:693071/1‑151 (MQ=255)
                                                                                                                    gAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAAcg                           >  2:1003771/1‑151 (MQ=255)
                                                                                                                    gAATGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTACGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAAc                            <  1:53379/150‑1 (MQ=255)
                                                                                                                     aaTGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAAcgc                          <  1:1136244/151‑1 (MQ=255)
                                                                                                                     aaTGAGGAGGATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGCTATGCGTCTGAACGACCGTGCGGGTGAAAcgc                          >  2:596809/1‑151 (MQ=255)
                                                                                                                         aggaggATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCt                       >  2:970612/1‑150 (MQ=255)
                                                                                                                         aggaggATAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTa                      <  2:477430/151‑1 (MQ=255)
                                                                                                                          ggaggaTAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTAc                     <  1:1122939/151‑1 (MQ=255)
                                                                                                                          ggaggaTAACAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTAc                     <  1:537429/151‑1 (MQ=255)
                                                                                                                                  aCAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATcc             >  1:754455/1‑151 (MQ=255)
                                                                                                                                   cAGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCt            <  2:391310/151‑1 (MQ=255)
                                                                                                                                    aGCTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTg           >  2:1294918/1‑151 (MQ=255)
                                                                                                                                      cTCAGAGGGAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGt          >  1:547243/1‑150 (MQ=255)
                                                                                                                                            gggAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGa    >  2:934977/1‑150 (MQ=255)
                                                                                                                                             ggAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAAc  >  2:1107101/1‑151 (MQ=255)
                                                                                                                                              gAAAGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAAc  >  2:663507/1‑150 (MQ=255)
                                                                                                                                               aaaGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGa    >  2:1149902/1‑147 (MQ=255)
                                                                                                                                               aaaGTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAAc  >  2:1037308/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GATGAATCTCCAGCAGATGCCTTCAGATAAGGTGATGCAAGGATTGCAGCTGGAGCTTTCGAATATCAACTAAAACAACAAGGCAAAAACTGGTATTTGACAAGTTATCAGTCAGTGAATGAGGAGGATAACAGCTCAGAGGGAAATTCAGAATCAACAGAAACGCCTACCCCCGTCCTGCGCCCCGAAGCTGGAAGTTACGTCGCTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAAC  >  NC_000913/1221092‑1221383

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