Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,004,984 +CG coding (1516/1833 nt) recQ → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,004,9831.G99.4% 612.6 / ‑4.1 174coding (1515/1833 nt)recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (90/83);  total (90/84)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
*NC_0009134,004,9832.C99.4% 611.7 / ‑3.9 174coding (1515/1833 nt)recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (90/83);  total (90/84)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 192 total aligned reads displayed.

CCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTCGCG‑‑GAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  NC_000913/4004893‑4005060
                                                                                               ||                                                                             
 cGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCg‑‑gc                                                                             >  1:2253464/1‑92 (MQ=255)
 cGCCAGCTGATTCACCTCGGCCTGGTAACGCAAAATATTGCCCAGAATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCg‑‑gc                                                                             <  1:4748485/93‑2 (MQ=255)
  gCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCg‑‑gcg                                                                            >  2:2213893/1‑91 (MQ=255)
  gCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCg‑‑gcg                                                                            >  1:5045380/1‑91 (MQ=255)
   ccAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAAGTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCg‑‑gcga                                                                           >  1:5149951/1‑90 (MQ=255)
   ccAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCg‑‑gcga                                                                           >  1:1698713/1‑90 (MQ=255)
   ccAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCACTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCg‑‑gcga                                                                           >  2:1108254/1‑90 (MQ=255)
    cAGCTGATTCACCTCGGCCTGGTGATGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCg‑‑gcgaa                                                                          <  2:3151120/93‑5 (MQ=255)
          aTTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCctct                                                                      >  2:5530869/1‑93 (MQ=255)
          aTTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCctct                                                                      <  1:609470/93‑1 (MQ=255)
           ttCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCtt                                                                     <  2:2150175/93‑1 (MQ=255)
           ttCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCtt                                                                     >  1:5648801/1‑93 (MQ=255)
           ttCACCTCGGCCGGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCtt                                                                     <  1:3918640/93‑1 (MQ=255)
            tCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCttt                                                                    <  2:4258183/93‑1 (MQ=255)
              cccTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGc                                                                  <  2:1764047/92‑1 (MQ=255)
              aCCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGc                                                                  >  2:3015734/1‑93 (MQ=255)
                 tCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAAc                                                               >  2:5626980/1‑93 (MQ=255)
                 tCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAAc                                                               >  2:647907/1‑93 (MQ=255)
                  cGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACt                                                              <  1:485284/93‑1 (MQ=255)
                   ggCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACtt                                                             <  1:2437194/93‑1 (MQ=255)
                    gCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtg                                                            >  2:4415834/1‑93 (MQ=255)
                      cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgcc                                                          >  1:2557025/1‑93 (MQ=255)
                      cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgcc                                                          >  2:1945712/1‑93 (MQ=255)
                       tGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgccg                                                         <  1:1366513/93‑1 (MQ=255)
                       tGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgccg                                                         >  1:4225600/1‑93 (MQ=255)
                        ggTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgt                                                        <  2:3845482/93‑1 (MQ=255)
                        ggTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgt                                                        <  2:5968971/93‑1 (MQ=255)
                           gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgcc                                                     <  2:3048351/93‑1 (MQ=255)
                           gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgcc                                                     >  1:737909/1‑93 (MQ=255)
                              gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCcgcg                                                  >  2:4896244/1‑93 (MQ=255)
                               cAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGt                                                 >  1:2814729/1‑93 (MQ=255)
                                aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                >  2:4664803/1‑93 (MQ=255)
                                aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                >  1:478854/1‑93 (MQ=255)
                                aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                >  2:1683698/1‑93 (MQ=255)
                                aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                <  2:3710425/93‑1 (MQ=255)
                                 aaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTAt                                               <  2:3489286/93‑1 (MQ=255)
                                  aaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATc                                              <  2:4431861/93‑1 (MQ=255)
                                  aaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATc                                              >  2:3127768/1‑93 (MQ=255)
                                   atatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCg                                             <  2:2935356/93‑1 (MQ=255)
                                   atatTCCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCg                                             <  1:2727634/93‑1 (MQ=255)
                                    tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGt                                            >  1:5610065/1‑93 (MQ=255)
                                    tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGt                                            >  2:2648882/1‑93 (MQ=255)
                                     atTGCCCAGCATTGTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTg                                           <  2:498025/93‑1 (MQ=255)
                                     atTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTCTCGTg                                           >  2:6066578/1‑93 (MQ=255)
                                     atTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTg                                           >  1:4110116/1‑93 (MQ=255)
                                     atTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGAATCGTg                                           >  1:1030131/1‑93 (MQ=255)
                                       tGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgc                                         <  2:68579/93‑1 (MQ=255)
                                         cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcgc                                       >  1:4983243/1‑93 (MQ=255)
                                         cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcgc                                       >  1:1902402/1‑93 (MQ=255)
                                          ccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                                      >  1:4326206/1‑93 (MQ=255)
                                          ccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                                      >  2:5797120/1‑93 (MQ=255)
                                            aGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCa                                    >  2:1363504/1‑93 (MQ=255)
                                            aGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCa                                    >  2:42412/1‑93 (MQ=255)
                                              cATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCaaa                                  <  1:5977428/93‑1 (MQ=255)
                                              cATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCaaa                                  <  1:3342355/93‑1 (MQ=255)
                                               aTTCTGCCATACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 <  1:1477040/93‑1 (MQ=255)
                                                ttCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAcc                                <  2:2582152/93‑1 (MQ=255)
                                                 tCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCg                               <  1:4527807/93‑1 (MQ=255)
                                                  cTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGa                              >  2:4790687/1‑93 (MQ=255)
                                                   tGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaa                             >  2:3797509/1‑93 (MQ=255)
                                                   tGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaa                             <  2:4809197/93‑1 (MQ=255)
                                                     cccTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAg                           <  2:3400904/93‑1 (MQ=255)
                                                     cccTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAg                           <  2:141115/93‑1 (MQ=255)
                                                     cccTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAg                           <  1:5771225/93‑1 (MQ=255)
                                                     cccTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAg                           <  2:141302/93‑1 (MQ=255)
                                                      ccTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          <  1:4644784/93‑1 (MQ=255)
                                                      ccTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          >  1:892257/1‑93 (MQ=255)
                                                      ccTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          >  2:368662/1‑93 (MQ=255)
                                                      ccTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          >  2:5234938/1‑93 (MQ=255)
                                                       cTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGTCGCTCAAACCGAAAGCg                         >  2:201437/1‑93 (MQ=255)
                                                       cTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCg                         >  1:3724428/1‑93 (MQ=255)
                                                       cTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCg                         <  1:6040366/93‑1 (MQ=255)
                                                        tACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        >  1:1510820/1‑93 (MQ=255)
                                                          cAACTGGCAGAGGCCGCGCGCCCGATGCTGCTCGGC‑‑GAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATg                      <  2:2111542/93‑1 (MQ=255)
                                                          cAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATg                      >  1:2987551/1‑93 (MQ=255)
                                                           aaCTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGc                     >  2:5389368/1‑93 (MQ=255)
                                                           aaCTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGc                     <  2:1693106/93‑1 (MQ=255)
                                                            aCTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    <  2:3896197/93‑1 (MQ=255)
                                                            aCTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    <  1:1233682/93‑1 (MQ=255)
                                                              tGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGa                  >  2:3787711/1‑93 (MQ=255)
                                                               gACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAATGCGATGCAGaa                 <  1:3685651/93‑1 (MQ=255)
                                                               gACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGaa                 <  2:2396815/93‑1 (MQ=255)
                                                                aCAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGaaa                <  1:3131271/93‑1 (MQ=255)
                                                                 cAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAAt               <  1:3622664/93‑1 (MQ=255)
                                                                   gagGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCg             >  1:3996234/1‑93 (MQ=255)
                                                                     ggCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           <  2:2665705/93‑1 (MQ=255)
                                                                     ggCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           <  1:5082222/93‑1 (MQ=255)
                                                                     ggCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           >  2:3173523/1‑93 (MQ=255)
                                                                      gCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          <  1:6071625/93‑1 (MQ=255)
                                                                      gCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          <  2:5869253/93‑1 (MQ=255)
                                                                        cgcgcgCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcgg        >  2:4038532/1‑93 (MQ=255)
                                                                        cgcgcgCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcgg        <  1:1631922/93‑1 (MQ=255)
                                                                         gcgcgcCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       <  2:1830231/93‑1 (MQ=255)
                                                                          cgcgcCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggcg      <  1:6141102/93‑1 (MQ=255)
                                                                           gcgcCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggcgg     >  1:3410942/1‑93 (MQ=255)
                                                                            cgcCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggcgga    >  1:2500519/1‑92 (MQ=255)
                                                                              cccG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCaa  >  1:1046045/1‑93 (MQ=255)
                                                                              cccG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCaa  <  1:5406918/93‑1 (MQ=255)
                                                                              cccG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCaa  >  2:1215068/1‑93 (MQ=255)
                                                                              cccG‑‑GTGCTGCGCGGCGAATCCTATTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCaa  >  2:5425263/1‑93 (MQ=255)
                                                                                               ||                                                                             
CCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTCGCG‑‑GAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  NC_000913/4004893‑4005060

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