Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 404,848 A→G 14.9% L266P (CTC→CCC)  proC ← pyrroline‑5‑carboxylate reductase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913404,8480AG14.9% 167.5 / 4.7 79L266P (CTC→CCC) proCpyrroline‑5‑carboxylate reductase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (25/40);  new base G (11/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 4.52e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCATGTTACGGGTGCATGAAGGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGTGGCA  >  NC_000913/404624‑405061
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                       
gCCATGTTACGGGTTCATGAAGGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAg                                                                                                                                                                                                                        <  1:220788/226‑1 (MQ=255)
   aTGTTACGGGTGCATGAAGGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGGGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAACGGGCGGCCGCAGGTCCGGCGAAAGTCGACAGGGTTGTCTGGGgtttt                                                                                                                                                                                                                   >  1:675303/1‑228 (MQ=255)
         cGGTTGAATTAAGGGCAAAATTAGTTCCGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGGCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGAtt                                                                                                                                                                                                             <  2:189783/228‑1 (MQ=255)
         cGGGTGCATTAAGGGCTAAATACATTACGTTTGCGGAATGCGCCACAGTTAACTTTGCGGATTAGTGTGGGGGTTGCGTCGCTGAACCCACAAATGAGTCACTAGCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGAtt                                                                                                                                                                                                             <  2:540951/228‑1 (MQ=255)
              gCATGAAGGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAGGACATACGGAGTTTCTGGGGTTTTCTGATTTTTc                                                                                                                                                                                                         >  1:326376/1‑227 (MQ=255)
               gaTGAAGGGTTAAATACATTACGTTTGCCGAATAACCCAAAGGTAACTTTACGGATTAGGTTGGGGGTTGCGCCGCTGAACCCAAAAATGAGGCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCa                                                                                                                                                                                                       <  1:580631/227‑1 (MQ=255)
               cttGAAGTGTTAAATTCATTAGGTTTGCCGACTACGCCAACGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTAGCGGGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCa                                                                                                                                                                                                       <  2:298886/226‑1 (MQ=255)
               cATGAAGGGCTAAATTCATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCa                                                                                                                                                                                                       <  1:141006/228‑1 (MQ=255)
               cATGAAGGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCCGCTCAAGTCATCAGGGTTTGCTGGGTTTTTCTGATTTTTCCa                                                                                                                                                                                                       >  1:631873/1‑228 (MQ=255)
               cATGAAGGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCCGACGACACGCGAAAGACCACAAGAGTTGCTTGGGTTTTCTGTTTTTTCCa                                                                                                                                                                                                       >  2:32509/1‑228 (MQ=255)
               aaTGAAGGGCTAAATACATTGCGTTTGCCGAATGCGCCACAGGGAATGTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCGTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCa                                                                                                                                                                                                       <  2:530188/227‑1 (MQ=255)
                      ggCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGCCGTCCGGCGAAAGTCATCAGGATTTGCTGAGGTTTTCTGATTTTTCCATaacat                                                                                                                                                                                                 >  1:86097/1‑223 (MQ=255)
                                aTTACGTTTGCCGAATATGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGGGTTTTTCTGATTTTTCCATTAACTTCGTCATCGCt                                                                                                                                                                                      >  1:147519/1‑228 (MQ=255)
                                aTTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTTTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCt                                                                                                                                                                                      <  1:178569/228‑1 (MQ=255)
                                  tgcGTTTGCCGAAATCGCCACAGGTAACTGTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCAAAAATGAGTCACTATCGTGAGTGGTTGAAATCGGAAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCt                                                                                                                                                                                      <  2:86097/224‑1 (MQ=255)
                                  tACGTTTGTCCAATCCGCCACATGTAACGTTACGGATTAGTGTGGGGGTTGCGCCCCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCTCCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCt                                                                                                                                                                                      <  1:290418/226‑1 (MQ=255)
                                  tACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTTGCGGCCTGACGTCCGGCGAAAGACATCAGGATTTGCTAAGTTATTCTGCTTTTTCCCTACCCTTCGTAATCGCtc                                                                                                                                                                                     >  1:94490/1‑226 (MQ=255)
                                  tACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTTACGACCGGCGAAAGTCATCAGGCTTTGCTGGGTTTTTTTGATTTTTCCTGAAACTTCGTCATCCCtt                                                                                                                                                                                     >  1:637776/1‑227 (MQ=255)
                                     gTTTGACGAAAACGCCACAGGTAACTTTCCGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGTGCTTTACAACGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATTCACTTCGTCATCGCTTCGAt                                                                                                                                                                                 <  1:173184/228‑1 (MQ=255)
                                       ttGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAACTGAGTAACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGGAAGTCATCAGGATTTGCTTAGTTTTTCTGATTTTTCCATCAACCTCGTAATCGCTTCGTTCa                                                                                                                                                                               >  1:79714/1‑228 (MQ=255)
                                        tGCCGAATACGCCACCGGTAATTTTTCGGATTAGTTGGGGGTTTTGGCCGGTGAACCCAAAAATGAGGCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCAc                                                                                                                                                                              <  2:79714/228‑1 (MQ=255)
                                        tGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCCGACGTGCGGGGAAAGACATAAGGATTTGCGTAGGTTTTCGTGTTTTTCCCTAACATTCGCCAACGCTTCGATAAc                                                                                                                                                                              >  2:178569/1‑228 (MQ=255)
                                          ccGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTTATGAGGTTTTCTGATTTTTCCATAAACTTCCTCATGCCTTCGATTACt                                                                                                                                                                             >  1:397617/1‑227 (MQ=255)
                                                    cacaGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATAAAATTCGTCATCGCTTCGCACACTGCCGCACGGa                                                                                                                                                                 >  1:564140/1‑227 (MQ=255)
                                                          taagtgtgcGGATTAGTGTGGGGGTGGCGCAGGTGAACCCACAAAGGAGTCACTATCGTGAGTGGTTGAAATCGGCAGAGTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCttt                                                                                                                                                          <  1:629149/220‑1 (MQ=255)
                                                                  cGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCCTCAGGATTTGCTGAGTTTTTCTGATTTTTCCAAACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTcc                                                                                                                                                   >  1:384053/1‑227 (MQ=255)
                                                                  cGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCCCTTCGTTAACTGGAGCACGGAAGCCTTTCTCTTcc                                                                                                                                                   >  1:30975/1‑227 (MQ=255)
                                                                   ggATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGTCGTAAACGCACCGAAGTGGCGGCCTGCACGCCGGGGACAAGAACACAGATTTGCTAGGTTTTTCTGATTTTTCCATACCATTCGTCATACCGTCGTTCCCTCAGCCACGGAATCCTTTTTCTTCCa                                                                                                                                                  >  2:660281/1‑227 (MQ=255)
                                                                   ggATAATTGTGGGGGTTGCGCCGATGAACGCACAAATTAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGAt                                                                                                                                                                                 <  1:373538/198‑1 (MQ=255)
                                                                       tAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGGGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGGTCACTGCAGCACGGAAGCCTTTCTCTTCCAGTAcc                                                                                                                                             >  1:359598/1‑227 (MQ=255)
                                                                          tgtgGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAACGCAAAGAAAGCCGGGCGTACCCGCACCGAAGTGCCGGCCTGCCGTCCGGCGAAAGTCATCCGGCTTTGCTGCGTTTTTCTGATCTTTTCACACACTTCGCCCTAGCCTTGACCCCTGCACCACGCGAGCCCTTCTCTTCCCGTTcg                                                                                                                                             >  2:116597/1‑225 (MQ=255)
                                                                            tgGGGTTTGCACCGCTGAACCCAAAAATGTGCAACCAGCGTGAAGGGGTGAACTCGGCCGAGTTGGTGCCTTACAATGCAAAGAAAACCGGACGTAACCGCGACGAATGGGCGTCCTGGCGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTTATTTTTCCATACCCTTCGTCATCTCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTCCGCGTAc                                                                                                                                        <  2:675303/228‑1 (MQ=255)
                                                                            tgGGGTTTGACCCGCTGAACCCAAAAATGATTCCATATCGTGAGTGGTTGAAATCGGCAGTTTTGGGGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTTACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACCCTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTAc                                                                                                                                        <  2:637776/228‑1 (MQ=255)
                                                                            tgGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTAc                                                                                                                                        <  1:471218/228‑1 (MQ=255)
                                                                               gggTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGGCCGGCGAAAGTCATCAGGATTTTCGGAGGTTTTCTGATTTTTCCATACACTTCGTCCTCCCCTTCATTACTTCGCCACGGAAACCTTTTCCTTCCAGTAccctgccccc                                                                                                                                     >  2:524250/1‑219 (MQ=255)
                                                                                    gcgcCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAATAACGCCGGACGTAACCGCACCGAAGTGGCGGCCTCGAGTCCGGCGAACGGCATTCGGCGTTGCTGGGGTTTGCTGGTTTTTCCATACACCTTCTCATCGCGTAGATCCATGTAGCCCCGAAGCCTTTCTCTTCCCAGACCCGGACAGCTTCa                                                                                                                                 >  2:362755/1‑227 (MQ=255)
                                                                                            gAACCCACAAGTGGGTCCCTGTCGTGAGTGGTTGAAATCGGCAGATTTGGCCCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATggtggt                                                                                                                         <  1:281896/227‑1 (MQ=255)
                                                                                                       tGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGGGGTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCCGCACGGAAGCCTTTCTCTTCTACTACGCGTACCGCTTCAATGCGGGTGCCTCCCGGTGa                                                                                                             >  1:596971/1‑228 (MQ=255)
                                                                                                         aGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGc                                                                                                           <  1:440858/228‑1 (MQ=255)
                                                                                                         aGTCACTATCGTGAGTGGTTGAAAGCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGGAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAg                                                                                                            <  1:58653/227‑1 (MQ=255)
                                                                                                          gTCACTATCGCGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGGAAGGCATCCGGAATTGCTGAGTCTTTCTGGTTTTTCCATCCACTTCCTCCTCGCTTCGCTCCCTCCACCCGCGGAACCCTTTCTTTCCAGGACCCCTACCGCTTCAATTGTGGGGTCTCCCggggcgc                                                                                                           >  2:9840/1‑222 (MQ=255)
                                                                                                           tCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCCTCCGGGTTTGCTGAGTTTTTCTGATTTTTCCATACACCTCGTCATACCTTCGATCACTGCAGCACGGAATCCTTTTTCTTCCAGTAACCCTACCCCTTCAAATGTGGTGTCTCCCCGTggataa                                                                                                         >  2:440858/1‑223 (MQ=255)
                                                                                                            tccTATCGTTAGTGGTTGAAAGCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACCTAACCGCACCGAAGTGGCTGCCTGACGGCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGa                                                                                                        <  2:385098/226‑1 (MQ=255)
                                                                                                              cTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGGCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACCCGTACCGCTTCTATGGGGGGGGCCCCCGGTggagcgacc                                                                                                      >  2:93163/1‑220 (MQ=255)
                                                                                                               tagcGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCa                                                                                                          <  1:93163/220‑1 (MQ=255)
                                                                                                                           ttGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGGGAAAGTTATCAGGATTTGCTGAGGTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCATCACGGTAACCTTTCTCTTCCATGACCAGTACCCCGTCACGTGTGGTGCCTCCCCGGGAGCAGACCACATCTTTCAGTg                                                                                         >  1:218626/1‑228 (MQ=255)
                                                                                                                           ttGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGGGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGACTTTTTTTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCCCCCGGTGAGCAGTCCATATCTTTTAGTg                                                                                         >  1:46430/1‑228 (MQ=255)
                                                                                                                           gtGAAATCGGGTGTTTTGGCCATTGCAAAAGCAAAGAACGCCGGACGTAACCGCGCCGAAGTTGCGGCCTGACGTCCGGCGAAAGTCAGCAGGATTTGCTGAGTTTTTCTGATTTTTCCATCCACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTg                                                                                         <  2:631873/227‑1 (MQ=255)
                                                                                                                                   gggaGATTTGGGCCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCCAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCgg                                                                                  <  2:147519/224‑1 (MQ=255)
                                                                                                                                   ggCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGGAAGTCATCCGGGTTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGGTCCCTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCCCTTCAATGGTGGTGCCTCCCGGGGAGCAGACCCACTCTTTCAGTGCCCCCggc                                                                                 >  1:398563/1‑227 (MQ=255)
                                                                                                                                    gCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGGCCGGCGAAAGTCATCAGGATTTGCGGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTAGATCACTGCAGCACGCAAGTCTTTCTTTTCCCTACCCCTTCCCCCTTCACTAGTGGTGCCCCCCCGTCCTCCAACCATATCTTTCACTGCCCCCggtt                                                                                >  2:676086/1‑226 (MQ=255)
                                                                                                                                     cAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGGCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTTCCGCACGGAAGCCTTTCTCTTCCCGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGGGAGCGGACCACATCTTTCATGGCCCCCGGaa                                                                                >  1:93548/1‑226 (MQ=255)
                                                                                                                                     cAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTACGGCGAAAGTAATCAAGATTTGCTGAGGTTTTCTGATTTTTCCCTACACTTCGTCATCGCTTCGGTAACTTCACCCCCGTAGCCCTTCTCTTCCAGATCCCGGTACGCTTCAATGTGGGGGCCTCCCCGTGAGCCGACCCTCACTTTCCCTGCCCCCCGaac                                                                               >  2:690177/1‑226 (MQ=255)
                                                                                                                                      aGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCGGTCGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATg                                                                               <  1:9840/227‑1 (MQ=255)
                                                                                                                                                       aagcaaagAAAGCCGGACGTAACCGCACCTAAGTGGCGGCCCGAGGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTGTCCATACACTTCGTCATCGCGTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCAc                                                             <  1:406528/228‑1 (MQ=255)
                                                                                                                                                             agaaagCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGGTTTGCTGAGTTTTTCTGATTTTTCCATAAACTTCGTCATCGCTTCGATCACTTCACCACGGAAGCCTTTTTCTTCCACTACCCGTCCCGCTTCAATGGCGGGGGCTCCCGGTGAGTAGCCCATATCTTTCCGTGCCCCCGCGGGTCTTCCCCTTTCCACCCCCCttt                                                        >  2:629041/1‑227 (MQ=255)
                                                                                                                                                             agaaagCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATAAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAATCCTTTTTCTTCCAGTACGCGTACCGCTTCAATGGGGGTGACTCCCGGTGAGCAGGCCCTATCTTTCCGTTACCCCGGATGTTTTCCCGTTTCCAGCAcccattt                                                       >  1:294238/1‑223 (MQ=255)
                                                                                                                                                                 agcccgccttAACCCCCCCGGAGTTGCGGCCGTGCTTCCGGCGAAAGTCATCAGGTTTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCg                                                   <  2:93548/219‑1 (MQ=255)
                                                                                                                                                                 agCCGGGCGTTACCTGCCCGGAGGTGCGGCCTTGCGTCCGTCGAAAGTCATCAGGATTTGCTGAGTTTTTCTTATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCGGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCg                                                   <  2:46430/228‑1 (MQ=255)
                                                                                                                                                                 agCCGGGAGTAACCGCACCGAGGTGGCGGCCCTTCGTCGCGCTAAAGCCATCAGGATTTGCCGAGTTTGTCTGATTTTTCCATACACTGCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGTGTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCg                                                   <  2:94102/228‑1 (MQ=255)
                                                                                                                                                                 agCCGGACGTAACCGCACCGAAGTGGCTGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCg                                                   <  1:1972/228‑1 (MQ=255)
                                                                                                                                                                 agCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGGCAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTCCCATACCCTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCg                                                   <  2:398563/228‑1 (MQ=255)
                                                                                                                                                                 agCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCg                                                   <  1:629041/228‑1 (MQ=255)
                                                                                                                                                                  gCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCCATGCCCCCGGATGTGCTCCCGTTTCCAGCACCCTTTTTTCg                                                   >  1:459147/1‑227 (MQ=255)
                                                                                                                                                                   ccGGGCGTAACCGCACCGAAGTGGCGGCCTTACGTCCGGCGAAAGTCACCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGa                                                 <  2:359598/228‑1 (MQ=255)
                                                                                                                                                                   ccGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAg                                                                                                            <  2:282173/169‑1 (MQ=255)
                                                                                                                                                                   ccGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACCCGTACCGCTTCAATGGTGGGGCCTCCCGGTGAg                                                                                                            >  1:282173/1‑169 (MQ=255)
                                                                                                                                                                    cGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCCGATGTTCTCCCGTTTCCAGCCCCATTTTTGc                                                    >  1:94102/1‑224 (MQ=255)
                                                                                                                                                                     ggACGTAACCGTATCGAAGTGTCCGCCTGACCTTCTGCCAAAGTCACCAGGATTTGCTGGTTTTTTCTGATTTTTCCATCCACTTCGGCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAAc                                               <  2:185837/228‑1 (MQ=255)
                                                                                                                                                                      gtcgctaCCCCACCGAATTTGTGGCCTGACGTCCGGCGAAAGGCACCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACCCTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAAc                                               <  2:459147/221‑1 (MQ=255)
                                                                                                                                                                                cACCGAAGTGGCGGCCTGACGTCCGGCGAAAGACATCAGGATTTGCTGAGTTTTTCTGCGTTTTCCATACACTTCGTCATCGCTTCGCTCACTGCAGCCCGGAAGTCCTTCTCTTCCAGTACGCGTACCGCCTCACTGGGTGGCACTCCCGGTGAGCCGACCAAATCTTTCACGGCCCCCCGAGGTTCTCCCGCTTCCCGAACCAATTTTCCCGAAACCaacacctgt                                    >  2:233257/1‑220 (MQ=255)
                                                                                                                                                                                 aCCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCg                                                   <  1:233257/212‑1 (MQ=255)
                                                                                                                                                                                 aCCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTGTTTCTGATTTTTCCATACACTTCGTCATCGCTTTGATCACTGCAGCCCGGAAGCCTTTCTCTTCCGTTACCCGTACCGCTTCCATGGCGGTGCCTCCCGGTGAGCCGACCACTACGTTCAGTGCCCCCGGAAGTTGCCCCGTTTCCCGCAACCCTTTTTCGGAAACCACTACAGccc                                   >  2:296469/1‑227 (MQ=255)
                                                                                                                                                                                     aaGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCCTTACCGCCTGAGc                               >  1:385766/1‑228 (MQ=255)
                                                                                                                                                                                       gTGGCGCCCTGGCGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTGTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCgg                             <  1:676086/228‑1 (MQ=255)
                                                                                                                                                                                               cTGACGTCCGGCGAAAGTCATCAGGCTTTGCTGAGTTTTTCTGATTTTTCCATTCCCTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTta                     <  1:32509/228‑1 (MQ=255)
                                                                                                                                                                                               cTGACGTCCGGCGAAAGGCATCAGGTTTTTCTGAGTTTTTCTGATTTTTCCATACCCTTCGGCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTta                     <  2:596971/228‑1 (MQ=255)
                                                                                                                                                                                                 ttccGTCCGGCGAAAGTCCCTAGGATTTGCGCAGTTTTTCTGATTTTCCCTCTCCCCTCGTCATCGCTTCGATCCCTGCAGTACGGACGCCTTTCTCTTCCAGTACGCGTACCGCTGCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGGGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTtat                    <  2:294238/225‑1 (MQ=255)
                                                                                                                                                                                                                    cgctgggtttGGTGGGTTTTTCTGATTTTTCCATCCCCTCCGTCAGCGCTTCGCTCCCTGCAGCCCGGAAGCCTGTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGTGGCa  <  2:431936/221‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                       
GCCATGTTACGGGTGCATGAAGGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGTGGCA  >  NC_000913/404624‑405061

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