Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,844,348 G→T 8.3% A129S (GCC→TCC)  adeD → cryptic adenine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,844,3480GT 8.3% 334.3 / 17.2 96A129S (GCC→TCC) adeDcryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base G (45/43);  new base T (4/4);  total (49/47)
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.95e-01

AACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCAT  >  NC_000913/3844203‑3844494
                                                                                                                                                 |                                                                                                                                                  
aaCGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGcc                                                                                                                                                  >  4:847316/1‑148 (MQ=255)
                tttATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCgg                                                                                                                                                <  4:811058/134‑1 (MQ=255)
                tttATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCgg                                                                                                                                                >  3:811058/1‑134 (MQ=255)
                  tATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGcc                                                                                                                                <  3:847316/148‑1 (MQ=255)
                       aTGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCACTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCtg                                                                                                                           >  3:56504/1‑148 (MQ=255)
                        tGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGc                                                                                                                             <  2:249097/145‑1 (MQ=255)
                        tGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGc                                                                                                                             >  1:249097/1‑145 (MQ=255)
                            cACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTGGCCTGGTTTGCCCGCtg                                                                                                                           >  4:56212/1‑143 (MQ=255)
                            cACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCtg                                                                                                                           <  3:56212/143‑1 (MQ=255)
                             aCCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCACTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGa                                                                                                                     <  4:56504/148‑1 (MQ=255)
                                  cATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACagg                                                                                                                >  4:419309/1‑148 (MQ=255)
                                    tatTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACaggca                                                                                                              >  2:566641/1‑148 (MQ=255)
                                    tatTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACaggca                                                                                                              <  3:419309/148‑1 (MQ=255)
                                       tGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGcc                                                                                                                                <  3:425867/127‑1 (MQ=255)
                                       tGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGcc                                                                                                                                >  4:425867/1‑127 (MQ=255)
                                       tGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACaggcaagg                                                                                                           >  2:443543/1‑148 (MQ=255)
                                       tGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACaggcaagg                                                                                                           >  4:366980/1‑148 (MQ=255)
                                         aaTCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGcc                                                                                                                                         <  1:902910/116‑1 (MQ=255)
                                         aaTCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGcc                                                                                                                                         >  2:902910/1‑116 (MQ=255)
                                           tCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCaaa                                                                                                       <  3:366980/148‑1 (MQ=255)
                                            ccAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGttt                                                                                                                                   <  1:943678/119‑1 (MQ=255)
                                            ccAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGttt                                                                                                                                   >  2:943678/1‑119 (MQ=255)
                                                cATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAATCCGGCTTCGcc                                                                                                                                         >  2:468244/1‑109 (MQ=255)
                                                cATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAATCCGGCTTCGcc                                                                                                                                         <  1:468244/109‑1 (MQ=255)
                                                    atgaCGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGAGGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACt                                                                                              >  4:20116/1‑148 (MQ=255)
                                                     tgaCGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACtt                                                                                             <  1:443543/148‑1 (MQ=255)
                                                       aCGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTgg                                                                                                                                      <  1:1084303/105‑1 (MQ=255)
                                                       aCGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTgg                                                                                                                                      >  2:1084303/1‑105 (MQ=255)
                                                          ccGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCg                                                                                                                                           >  4:633859/1‑97 (MQ=255)
                                                          ccGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCg                                                                                                                                           <  3:633859/97‑1 (MQ=255)
                                                            ggTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTc                                                                                      <  1:566641/148‑1 (MQ=255)
                                                                   tttGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTgg                                                                                                                                      >  1:623715/1‑93 (MQ=255)
                                                                   tttGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTgg                                                                                                                                      <  2:623715/93‑1 (MQ=255)
                                                                   tttGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGcc                                                                                                                                >  1:684147/1‑99 (MQ=255)
                                                                   tttGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGcc                                                                                                                                <  2:684147/99‑1 (MQ=255)
                                                                   tttGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACag                                                                                                                 >  2:108501/1‑114 (MQ=255)
                                                                   tttGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACag                                                                                                                 <  1:108501/114‑1 (MQ=255)
                                                                   tttGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAg                                                                                                  <  1:138161/129‑1 (MQ=255)
                                                                   tttGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAg                                                                                                  >  2:138161/1‑129 (MQ=255)
                                                                    ttGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGc                                                                                                                                          <  3:534356/88‑1 (MQ=255)
                                                                    ttGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGc                                                                                                                                          >  4:534356/1‑88 (MQ=255)
                                                                        aaCCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACaggcaa                                                                                                             >  4:1073245/1‑113 (MQ=255)
                                                                        aaCCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACaggcaa                                                                                                             <  3:1073245/113‑1 (MQ=255)
                                                                         aCCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTAc                                                                         >  4:1006447/1‑148 (MQ=255)
                                                                         aCCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTAc                                                                         >  2:804262/1‑148 (MQ=255)
                                                                          ccGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTc                                                                                      <  3:498066/134‑1 (MQ=255)
                                                                          ccGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTc                                                                                      >  4:498066/1‑134 (MQ=255)
                                                                          ccGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTAcc                                                                        >  1:169049/1‑148 (MQ=255)
                                                                            gCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAATCCGGCTTCGccca                                                                                                                                       <  2:811638/83‑3 (MQ=255)
                                                                            gCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAATCCGGCTTCGccca                                                                                                                                       >  1:811638/1‑81 (MQ=255)
                                                                              tACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGcc                                                                    <  1:804262/148‑1 (MQ=255)
                                                                                   tGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACagg                                                                                                                <  4:500102/99‑1 (MQ=255)
                                                                                   tGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACagg                                                                                                                >  3:500102/1‑99 (MQ=255)
                                                                                   tGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGa                                                               <  3:20116/148‑1 (MQ=255)
                                                                                    gCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACt                                                                                              <  4:989194/116‑1 (MQ=255)
                                                                                    gCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACt                                                                                              >  3:989194/1‑116 (MQ=255)
                                                                                    gCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACt                                                                                              <  2:562252/116‑1 (MQ=255)
                                                                                    gCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACt                                                                                              >  1:562252/1‑116 (MQ=255)
                                                                                           ggacTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCg                                                       <  3:1006447/145‑1 (MQ=255)
                                                                                             ccTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTg                                                                                                                                  >  3:15150/1‑71 (MQ=255)
                                                                                             ccTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTg                                                                                                                                  <  4:15150/71‑1 (MQ=255)
                                                                                                gacgacCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAg                                                             <  4:80416/137‑1 (MQ=255)
                                                                                                gacgacCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAg                                                             >  3:80416/1‑137 (MQ=255)
                                                                                                 acgacCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACaggca                                                                                                              >  1:825371/1‑87 (MQ=255)
                                                                                                 acgacCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACaggca                                                                                                              <  2:825371/87‑1 (MQ=255)
                                                                                                  cgacCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTaa                                                >  2:198762/1‑148 (MQ=255)
                                                                                                    acCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGTCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTa                                                 <  3:985548/145‑1 (MQ=255)
                                                                                                    acCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTa                                                 >  4:985548/1‑145 (MQ=255)
                                                                                                     ccGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGa                                                               <  1:997588/130‑1 (MQ=255)
                                                                                                     ccGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGa                                                               >  2:997588/1‑130 (MQ=255)
                                                                                                          aTTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTGGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGAt                                                     >  2:916864/1‑135 (MQ=255)
                                                                                                          aTTTGCGACCCTCATGAAATCGTCAACGTCATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGGAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGAt                                                     <  1:916864/135‑1 (MQ=255)
                                                                                                                          aaaTCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTa                                                                                            >  3:19822/1‑80 (MQ=255)
                                                                                                                          aaaTCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTa                                                                                            <  4:19822/80‑1 (MQ=255)
                                                                                                                          aaaTCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACa                                                                                          >  4:309174/1‑82 (MQ=255)
                                                                                                                          aaaTCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACa                                                                                          <  3:309174/82‑1 (MQ=255)
                                                                                                                            aTCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGa                      <  2:169049/148‑1 (MQ=255)
                                                                                                                            aTCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGa                      >  4:886120/1‑148 (MQ=255)
                                                                                                                             tCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGctc                                                                                 >  1:867288/1‑88 (MQ=255)
                                                                                                                             tCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGctc                                                                                 <  2:867288/88‑1 (MQ=255)
                                                                                                                                  aaCGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCg                <  3:886120/148‑1 (MQ=255)
                                                                                                                                   aCGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGcc                                                                    >  4:1116569/1‑95 (MQ=255)
                                                                                                                                   aCGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGcc                                                                    <  3:1116569/95‑1 (MQ=255)
                                                                                                                                      tGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCtt                             <  3:816460/131‑1 (MQ=255)
                                                                                                                                      tGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCtt                             >  4:816460/1‑131 (MQ=255)
                                                                                                                                        ccaggCGAATCCGGCTTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATg                                                    >  3:657561/4‑106 (MQ=255)
                                                                                                                                        ccaggCGAATCCGGCTTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATg                                                    <  4:657561/103‑1 (MQ=255)
                                                                                                                                        aTGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGc          >  4:13540/1‑148 (MQ=255)
                                                                                                                                         tGGGCGAAGCCGGATTCGCCTGGCTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCt                                                          >  3:963664/1‑99 (MQ=255)
                                                                                                                                         cgggCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTCACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCt                                                          <  4:963664/98‑1 (MQ=255)
                                                                                                                                         caggCGAATCCGGCTTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTa                                                 >  1:1061671/3‑108 (MQ=255)
                                                                                                                                         caggCGAATCCGGCTTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTa                                                 <  2:1061671/106‑1 (MQ=255)
                                                                                                                                             cGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAAc                                                                                                                   >  4:993831/1‑38 (MQ=255)
                                                                                                                                             cGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAAc                                                                                                                   <  3:993831/38‑1 (MQ=255)
                                                                                                                                               aaGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCa   >  2:466425/1‑148 (MQ=255)
                                                                                                                                                aGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCAt  <  1:198762/148‑1 (MQ=255)
                                                                                                                                                 gCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTa                                  <  2:554771/115‑1 (MQ=255)
                                                                                                                                                 gCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTa                                  >  1:554771/1‑115 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCAT  >  NC_000913/3844203‑3844494

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