Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,844,146 T→C G61G (GGT→GGC adeD → cryptic adenine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,844,1460TC98.8% 289.0 / ‑3.7 81G61G (GGT→GGCadeDcryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (36/44);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

TAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGC  >  NC_000913/3843999‑3844293
                                                                                                                                                   |                                                                                                                                                   
tAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGgcgt                                                                                                                                                   <  1:823988/150‑1 (MQ=255)
tAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGtt                                                                                                                                                  >  1:1017508/1‑151 (MQ=255)
 aGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTg                                                                                                                                                 <  1:768949/151‑1 (MQ=255)
     aggCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGgcgc                                                                                                                                             >  2:1011865/2‑151 (MQ=255)
       gCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAg                                                                                                                                           <  1:500460/151‑1 (MQ=255)
         tGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGa                                                                                                                                          >  2:352020/1‑150 (MQ=255)
           aaTACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAAt                                                                                                                                        >  1:206803/1‑150 (MQ=255)
           aaTACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATa                                                                                                                                       <  1:135361/151‑1 (MQ=255)
                cAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTg                                                                                                                                  <  2:547389/151‑1 (MQ=255)
                 aGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTg                                                                                                                                  >  1:1208012/1‑150 (MQ=255)
                 aGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGa                                                                                                                                 <  2:1189215/151‑1 (MQ=255)
                    aaTTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGc                                                                                                                              <  1:439336/151‑1 (MQ=255)
                    aaTTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGc                                                                                                                              >  1:327376/1‑151 (MQ=255)
                     aTTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCt                                                                                                                             <  1:559707/151‑1 (MQ=255)
                       tGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTcc                                                                                                                           >  2:916429/1‑151 (MQ=255)
                       tGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTcc                                                                                                                           >  2:596146/1‑151 (MQ=255)
                         ttAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCgg                                                                                                                         >  1:905752/1‑151 (MQ=255)
                         ttAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCgg                                                                                                                         >  1:726393/1‑151 (MQ=255)
                           aGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGc                                                                                                                        <  1:1104016/150‑1 (MQ=255)
                           aGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCt                                                                                                                       <  2:158511/151‑1 (MQ=255)
                             ccGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCttt                                                                                                                     <  2:417091/151‑1 (MQ=255)
                             ccGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCttt                                                                                                                     <  2:674109/151‑1 (MQ=255)
                                 ttCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGACCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGc                                                                                                                   >  2:1099925/1‑149 (MQ=255)
                                      gTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGa                                                                                                             <  2:206816/150‑1 (MQ=255)
                                      gTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGAt                                                                                                            <  1:298328/151‑1 (MQ=255)
                                        ggCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGAtt                                                                                                           <  1:574379/150‑1 (MQ=255)
                                           gACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGa                                                                                                                                 <  2:198826/125‑1 (MQ=255)
                                           gACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGa                                                                                                                                 >  1:198813/1‑125 (MQ=255)
                                              gCTGTTGCTGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTc                                                                                                    <  1:1148026/151‑1 (MQ=255)
                                                 gTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTcgcg                                                                                                 <  1:1138544/151‑1 (MQ=255)
                                                 gTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTcg                                                                                                   <  2:206139/149‑1 (MQ=255)
                                                   tGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCgcgg                                                                                                >  2:456615/1‑150 (MQ=255)
                                                     ccGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGgcgc                                                                                             <  1:501665/151‑1 (MQ=255)
                                                         ttatattatTGCTAATGTATCTATTCTCGAACTGATCAATGGCGGAGAAATTTACGGCCCACTTGCGATTAAAGGACGTTACATTGCAGGCGGTGGCGCAAAATACACTGATGCTCCGGCTTTGCCGAGGATTGCTGTTCGCGGCGCAAc                                                                                          >  2:945938/1‑150 (MQ=255)
                                                          tatattatTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAcg                                                                                         <  2:635625/150‑1 (MQ=255)
                                                           atattatTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAcggc                                                                                       <  2:806594/151‑1 (MQ=255)
                                                             attattGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAcggcg                                                                                      <  1:848074/150‑1 (MQ=255)
                                                              ttattGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAcggcgg                                                                                     >  1:994755/1‑150 (MQ=255)
                                                                  tGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGcc                                                                                 >  2:943145/1‑150 (MQ=255)
                                                                  tGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGcc                                                                                 >  1:300854/1‑150 (MQ=255)
                                                                    aTAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAgg                                                                              >  1:287465/1‑151 (MQ=255)
                                                                      aaTGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAggg                                                                             >  2:363346/1‑150 (MQ=255)
                                                                          tctcTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGgtt                                                                          >  2:57336/1‑149 (MQ=255)
                                                                           ctctATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGtttttt                                                                       >  2:536113/1‑151 (MQ=255)
                                                                            tctATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTAtt                                                                       >  1:1177728/1‑150 (MQ=255)
                                                                                    cGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCt                                                                 <  2:326645/148‑1 (MQ=255)
                                                                                     gACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCAcc                                                             >  2:125815/1‑151 (MQ=255)
                                                                                       ccTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTg                                                           <  2:185865/151‑1 (MQ=255)
                                                                                        cTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGc                                                          >  1:269156/1‑151 (MQ=255)
                                                                                          gATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCat                                                        <  2:221135/151‑1 (MQ=255)
                                                                                             cAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCatat                                                      <  1:470880/150‑1 (MQ=255)
                                                                                             cAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATAtt                                                     >  1:1083628/1‑151 (MQ=255)
                                                                                                  gCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATc                                                <  2:571802/151‑1 (MQ=255)
                                                                                                  gCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGGAt                                                 >  2:1130310/1‑150 (MQ=255)
                                                                                                   cGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATcc                                               <  2:588190/151‑1 (MQ=255)
                                                                                                    ggAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATcc                                               >  2:1240676/1‑150 (MQ=255)
                                                                                                      agaAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGc                                            <  2:657559/151‑1 (MQ=255)
                                                                                                      agaAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGc                                            <  2:440456/151‑1 (MQ=255)
                                                                                                       gaAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGc                                            <  1:355415/150‑1 (MQ=255)
                                                                                                          aTTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCCTATTGAATCCAGGatg                                         >  2:1086544/1‑150 (MQ=255)
                                                                                                            ttCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCatgat                                       <  1:385360/150‑1 (MQ=255)
                                                                                                              ccGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCatgatg                                      <  2:259158/149‑1 (MQ=255)
                                                                                                              ccGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTAACCTGCATATTGAATCCAGCAtgatga                                     <  1:1099993/150‑1 (MQ=255)
                                                                                                                    cAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGt                              <  1:1099909/151‑1 (MQ=255)
                                                                                                                     aaTTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTTATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGt                              >  1:1119235/1‑150 (MQ=255)
                                                                                                                       ttGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTATTGATGCTCACCCTCATATTGAATCCAGCATGATGGCGCCgg                               >  2:433176/1‑147 (MQ=255)
                                                                                                                        tgtgATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCa                            <  2:205365/149‑1 (MQ=255)
                                                                                                                         gtgATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGACTACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTTTTTATGCTCACCTGCATTTTTGATCCAGCATGATGACGCCGGGCACtt                         >  2:1019225/1‑151 (MQ=255)
                                                                                                                              tAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTTTTGATGCTCACCTGCATATTGAATACAGCATGATGACGCCGGTCACtttt                       >  1:898399/1‑148 (MQ=255)
                                                                                                                                 aGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAAcc                 <  2:852445/151‑1 (MQ=255)
                                                                                                                                 aGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTTATGCCCACCTGCCTATTGAATCCCGCATGATGACGCCGGTCCCTTTTGaaaac                 >  2:1235459/1‑149 (MQ=255)
                                                                                                                                  ggACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCg                <  1:352034/151‑1 (MQ=255)
                                                                                                                                  ggACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGAGTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATCCCGCCGGTCCCTTTTGacgccg                <  2:773025/151‑6 (MQ=255)
                                                                                                                                    aCGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCt              <  2:302910/151‑1 (MQ=255)
                                                                                                                                         aCATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCt         <  2:833265/151‑1 (MQ=255)
                                                                                                                                         aCATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCt         <  1:71319/151‑1 (MQ=255)
                                                                                                                                          cATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCt         <  1:916425/150‑1 (MQ=255)
                                                                                                                                          cATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGGCACTTTTGAAACCGCTAccc          >  2:323975/1‑149 (MQ=255)
                                                                                                                                           aTTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGGCAGGGTTTATTGATTCCCCCCTGGCTTTTTCAACCAGCATAATGACGCCGGTCCATTTTTGAACCGCCACCCt         >  2:171841/1‑149 (MQ=255)
                                                                                                                                           aTTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTTTTGATGCTCACCTGCATATTTAATGCAGCATGATGACGCCGGTCACTTTTGAAACCGc               >  2:127360/1‑143 (MQ=255)
                                                                                                                                             tGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGcc      <  1:1022269/150‑1 (MQ=255)
                                                                                                                                                cGGCGTTGGCGCAGAATACACTGATGCGCCGGCTTTGCAACGGATTGATGCTCGCGGCGAAACGGCGGTGCCAGGGTTTATTGATGCTCCCCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTGACCTTCcgcgc  <  2:848979/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGC  >  NC_000913/3843999‑3844293

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