Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,887,671 A→G 5.5% A107A (GCT→GCC rnd ← ribonuclease D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,887,6710AG 5.5% 262.3 / 3.6 89A107A (GCT→GCCrndribonuclease D
Reads supporting (aligned to +/- strand):  ref base A (44/40);  new base G (3/2);  total (47/42)
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.00e+00

ACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAG  >  NC_000913/1887525‑1887820
                                                                                                                                                  |                                                                                                                                                     
aCATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGgcagca                                                                                                                                                    <  1:858689/150‑1 (MQ=255)
aCATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGgcagca                                                                                                                                                    <  1:1667293/150‑1 (MQ=255)
aCATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGgcagca                                                                                                                                                    <  1:714792/150‑1 (MQ=255)
  aTCCGCCGCTGCGTATTCACACTGACGTTCGGGCAGCGGTCGGGCCAGCCCGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCTGGCGTCCGCAGAAGGCAGCAAg                                                                                                                                                  <  2:249870/150‑1 (MQ=255)
   tCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAgg                                                                                                                                                 <  1:6097/150‑1 (MQ=255)
   tCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAgg                                                                                                                                                 >  2:1422454/1‑150 (MQ=255)
   tCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAgg                                                                                                                                                 >  2:594651/1‑150 (MQ=255)
     cgccgcTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGAt                                                                                                                                               >  1:1484624/1‑150 (MQ=255)
         gcTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGc                                                                                                                                           >  2:804041/1‑150 (MQ=255)
          cTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCg                                                                                                                                          <  1:587798/150‑1 (MQ=255)
          cTGCGTATTCACACTGACGTTCGGTCAGCGGGCTGGCCAGCCAGTCGGGGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCg                                                                                                                                          <  2:948738/150‑1 (MQ=255)
             cGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCgtgt                                                                                                                                       >  1:211367/1‑150 (MQ=255)
               tATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtca                                                                                                                                     <  2:389332/150‑1 (MQ=255)
                aTTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaa                                                                                                                                    <  2:1398741/150‑1 (MQ=255)
                  tCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCTAGGAAGCGCTACCCCATGACTTCGGGCGTCCGCAGCAGGCAGCAAGGAGTTGCGTTtcaatc                                                                                                                                  >  1:285069/1‑150 (MQ=255)
                  tCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaatc                                                                                                                                  >  1:189557/1‑150 (MQ=255)
                   cacaCTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaatca                                                                                                                                 >  1:1212247/1‑150 (MQ=255)
                        tGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGt                                                                                                                            <  2:191599/150‑1 (MQ=255)
                        tGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGt                                                                                                                            >  2:1428181/1‑150 (MQ=255)
                         gACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCAGTCTTGTCCGGCGTAACGCCGGAATACTCTGCCACCAGGGAGGCGACACGCGATGCGATCGGGCGTCCGCTGCTGGGGGCAAGGAGTTGCGGGTCTATCAAggggg                                                                                                                           >  2:283530/1‑148 (MQ=255)
                          aCGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTtg                                                                                                                          <  2:1718594/150‑1 (MQ=255)
                             ttCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTgg                                                                                                                       <  2:1653412/150‑1 (MQ=255)
                                   ctgcgggcggGCCAGCCTGGCGGTGCGCGATTCACTCTTGTCCAGCGTAAGGCCGGAATTCTCTTCCACCATCGAAGCGTAACCGCATGACCTCGGGCGTCCGCAGAAGGCAGCAATGATTTGCGTGTCAATCAAGGGTTGTGGTAATTc                                                                                                                 <  2:289374/148‑1 (MQ=255)
                                    aGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCg                                                                                                                <  2:1151749/150‑1 (MQ=255)
                                      cGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGcc                                                                                                              <  2:589397/150‑1 (MQ=255)
                                       ggTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCa                                                                                                             <  2:1484624/150‑1 (MQ=255)
                                           tGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGa                                                                                                         <  1:641032/150‑1 (MQ=255)
                                            ggCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGAc                                                                                                        <  2:196546/150‑1 (MQ=255)
                                             gccagccaGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACa                                                                                                       >  1:1077288/1‑150 (MQ=255)
                                                    agTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGa                                                                                                >  2:381716/1‑150 (MQ=255)
                                                       cagTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCGGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAaca                                                                                             >  1:1710282/3‑150 (MQ=255)
                                                       cGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCGGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAaca                                                                                             >  1:1710280/1‑150 (MQ=255)
                                                        ggTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAacac                                                                                            >  1:767522/1‑150 (MQ=255)
                                                            cgcgATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTcc                                                                                        <  1:1627482/150‑1 (MQ=255)
                                                             gcgATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCa                                                                                       >  1:532337/1‑150 (MQ=255)
                                                             gcgATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCa                                                                                       <  2:43048/150‑1 (MQ=255)
                                                              cgATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCACATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACGCTTCCAg                                                                                      >  1:294439/1‑150 (MQ=255)
                                                               gATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGa                                                                                     >  1:1318729/1‑150 (MQ=255)
                                                                 ttCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATc                                                                                   <  1:365561/150‑1 (MQ=255)
                                                                     ctctTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCa                                                                               >  2:572673/1‑150 (MQ=255)
                                                                     ctctTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCa                                                                               <  2:767522/150‑1 (MQ=255)
                                                                      tctTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCAc                                                                              <  1:572673/150‑1 (MQ=255)
                                                                      tctTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCAc                                                                              >  2:948359/1‑150 (MQ=255)
                                                                      tctTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCAc                                                                              <  2:658707/150‑1 (MQ=255)
                                                                       ctTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCAct                                                                             <  2:1212247/150‑1 (MQ=255)
                                                                           tCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCActgcct                                                                         <  2:189557/150‑1 (MQ=255)
                                                                             cAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCActgcctgc                                                                       <  1:1430474/150‑1 (MQ=255)
                                                                                cGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATg                                                                    <  1:594651/150‑1 (MQ=255)
                                                                                 gTAACTCCGTAATACTCATCCACCATTGAAGCGAAACCCCATGACATCTGGCGTCCGCATAAGTCAGCAAGGATTTGCGTGTCAATCAAGGGTTGGGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATgg                                                                   <  1:1493717/150‑1 (MQ=255)
                                                                                     cGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGaa                                                               >  2:1521496/1‑150 (MQ=255)
                                                                                     cGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGaa                                                               >  1:1030603/1‑150 (MQ=255)
                                                                                      gCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGaaa                                                              >  2:374079/1‑150 (MQ=255)
                                                                                      gCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGaaa                                                              >  2:1139677/1‑150 (MQ=255)
                                                                                         ggAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAAttt                                                           >  1:548772/1‑150 (MQ=255)
                                                                                         ggAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAAttt                                                           >  2:1294735/1‑150 (MQ=255)
                                                                                         ggAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAAttt                                                           <  1:1521496/150‑1 (MQ=255)
                                                                                         ggAATACCCTTCCTCCATGGAATGGAAAACCCATGACCTCGGGCGTCCGCCGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACAATGAGGAACACTTCCAGATCTTCACTGCCTGCATGTAGAAAttt                                                           <  2:285069/150‑1 (MQ=255)
                                                                                           aaTACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTtg                                                         <  1:381716/150‑1 (MQ=255)
                                                                                           aaTACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTtg                                                         <  1:384304/150‑1 (MQ=255)
                                                                                             tACTCTTCCACCATGGAAGCGACACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTGGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTtgtg                                                       <  1:919893/150‑1 (MQ=255)
                                                                                               ctctTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGAt                                                     >  1:302930/1‑150 (MQ=255)
                                                                                                  ttCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGa                                                  >  1:469729/1‑150 (MQ=255)
                                                                                                     caccaTGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACgg                                               >  2:1665679/1‑150 (MQ=255)
                                                                                                       ccaTGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGAt                                             >  2:981545/1‑150 (MQ=255)
                                                                                                       ccaTGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCAAGGAGAAATTTTGTGATGGACGGaa                                             >  2:242925/1‑149 (MQ=255)
                                                                                                         aTGGAAGTGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATcg                                           >  1:107827/1‑150 (MQ=255)
                                                                                                         aTGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATcg                                           >  2:1308123/1‑150 (MQ=255)
                                                                                                          tGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATcgc                                          <  1:778969/150‑1 (MQ=255)
                                                                                                           ggAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATcgcg                                         >  2:381635/1‑150 (MQ=255)
                                                                                                            gAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCgcgc                                        >  2:1669745/1‑150 (MQ=255)
                                                                                                            gAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCgcgc                                        >  1:1288474/1‑150 (MQ=255)
                                                                                                            gAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCgcgc                                        >  2:1468227/1‑150 (MQ=255)
                                                                                                             aaGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCa                                       <  1:1665679/150‑1 (MQ=255)
                                                                                                              aGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAg                                      >  1:1159270/1‑150 (MQ=255)
                                                                                                                cGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGa                                    >  2:557533/1‑150 (MQ=255)
                                                                                                                 gAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACAATGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAAATATGTGATGGACGGATCGCGCAGGAt                                   >  1:256992/1‑150 (MQ=255)
                                                                                                                       ccATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCttt                             >  2:1674018/1‑150 (MQ=255)
                                                                                                                           gACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGc                         <  2:1077288/150‑1 (MQ=255)
                                                                                                                           gACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAAAACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGc                         >  2:1186767/1‑150 (MQ=255)
                                                                                                                                 gggCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGAc                   <  1:1674018/150‑1 (MQ=255)
                                                                                                                                 gggCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGAc                   <  2:1159270/150‑1 (MQ=255)
                                                                                                                                  ggCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGAcc                  <  1:1139677/150‑1 (MQ=255)
                                                                                                                                  ggCGTCCGCAGAAGGCAGAAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCATAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACTGAACGCGCACGATCGCTTTCAGCGGTGAcc                  >  1:1101707/1‑150 (MQ=255)
                                                                                                                                  ggCCTCCGCAGAAGGCGGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGAcc                  <  2:1710282/150‑1 (MQ=255)
                                                                                                                                  ggCCTCCGCAGAAGGCGGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGAcc                  <  2:1710280/150‑1 (MQ=255)
                                                                                                                                   gCGTCCGCAGAAGGCAGCAACGATTTGCGTTTCAATCAAGGGTTGTGTTATTTCGCCAAAGACATTCAGGAACACTGCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCa                 <  1:283530/150‑1 (MQ=255)
                                                                                                                                      tCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTc              >  2:1114928/1‑150 (MQ=255)
                                                                                                                                         gCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGt           >  1:534202/1‑150 (MQ=255)
                                                                                                                                         gCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGt           <  1:1422454/150‑1 (MQ=255)
                                                                                                                                          cAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTg          >  1:201755/1‑150 (MQ=255)
                                                                                                                                            gAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGAt        >  2:1522216/1‑150 (MQ=255)
                                                                                                                                              aGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCACTTTCAGCGGTGACCAGTCGGTGATcc      <  2:532337/150‑1 (MQ=255)
                                                                                                                                                 cagcaAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGa   >  1:719194/1‑150 (MQ=255)
                                                                                                                                                  agcaAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAg  >  2:710666/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAG  >  NC_000913/1887525‑1887820

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