Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,844,353 A→C 8.2% G130G (GGA→GGC adeD → cryptic adenine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,844,3530AC 8.2% 330.4 / 21.2 97G130G (GGA→GGCadeDcryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base A (45/44);  new base C (4/4);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.21e-01

TATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCAT  >  NC_000913/3844221‑3844494
                                                                                                                                    |                                                                                                                                             
tATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGcc                                                                                                                                <  3:847316/148‑1 (MQ=255)
     aTGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCACTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCtg                                                                                                                           >  3:56504/1‑148 (MQ=255)
      tGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGc                                                                                                                             >  1:249097/1‑145 (MQ=255)
      tGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGc                                                                                                                             <  2:249097/145‑1 (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                                                                                                                                >  4:425867/1‑127 (MQ=255)
                     tGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGcc                                                                                                                                <  3:425867/127‑1 (MQ=255)
                     tGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACaggcaagg                                                                                                           >  2:443543/1‑148 (MQ=255)
                     tGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACaggcaagg                                                                                                           >  4:366980/1‑148 (MQ=255)
                       aaTCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGcc                                                                                                                                         >  2:902910/1‑116 (MQ=255)
                       aaTCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGcc                                                                                                                                         <  1:902910/116‑1 (MQ=255)
                         tCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCaaa                                                                                                       <  3:366980/148‑1 (MQ=255)
                          ccAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGttt                                                                                                                                   >  2:943678/1‑119 (MQ=255)
                          ccAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGttt                                                                                                                                   <  1:943678/119‑1 (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                                                                                                                                <  2:684147/99‑1 (MQ=255)
                                                 tttGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGcc                                                                                                                                >  1:684147/1‑99 (MQ=255)
                                                 tttGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACag                                                                                                                 >  2:108501/1‑114 (MQ=255)
                                                 tttGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACag                                                                                                                 <  1:108501/114‑1 (MQ=255)
                                                 tttGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAg                                                                                                  >  2:138161/1‑129 (MQ=255)
                                                 tttGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAg                                                                                                  <  1:138161/129‑1 (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                                                                                      >  4:498066/1‑134 (MQ=255)
                                                        ccGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTc                                                                                      <  3:498066/134‑1 (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                                                                                              <  2:562252/116‑1 (MQ=255)
                                                                  gCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACt                                                                                              >  3:989194/1‑116 (MQ=255)
                                                                  gCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACt                                                                                              >  1:562252/1‑116 (MQ=255)
                                                                         ggacTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCg                                                       <  3:1006447/145‑1 (MQ=255)
                                                                           ccTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTg                                                                                                                                  <  4:15150/71‑1 (MQ=255)
                                                                           ccTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTg                                                                                                                                  >  3:15150/1‑71 (MQ=255)
                                                                              gacgacCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAg                                                             >  3:80416/1‑137 (MQ=255)
                                                                              gacgacCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAg                                                             <  4:80416/137‑1 (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                                                                                          <  3:309174/82‑1 (MQ=255)
                                                                                                        aaaTCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACa                                                                                          >  4:309174/1‑82 (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                                                                    <  3:1116569/95‑1 (MQ=255)
                                                                                                                 aCGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGcc                                                                    >  4:1116569/1‑95 (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)
                                                                                                                                   gATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGa                                                      <  2:339909/91‑1 (MQ=255)
                                                                                                                                   gATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGa                                                      >  1:339909/1‑91 (MQ=255)
                                                                                                                                    aTTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCg         <  3:203415/135‑1 (MQ=255)
                                                                                                                                    aTTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCg         >  4:203415/1‑135 (MQ=255)
                                                                                                                                    |                                                                                                                                             
TATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCAT  >  NC_000913/3844221‑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: