Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,944,115 A→G noncoding (412/2904 nt) rrlC → 23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,944,1150AG98.9% 536.6 / ‑1.2 185noncoding (412/2904 nt)rrlC23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (100/83);  minor base . (1/1);  total (101/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 187 total aligned reads displayed.

TGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGG  >  NC_000913/3943979‑3944198
                                                                                                                                         |                                                                                   
tGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCa                                                                                  <  1:1679452/140‑1 (MQ=37)
 gAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCAt                                                                                 >  1:310337/1‑140 (MQ=37)
 gAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCAt                                                                                 >  1:98747/1‑140 (MQ=37)
 gAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCAt                                                                                 >  2:1810770/1‑140 (MQ=37)
 gAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCAt                                                                                 >  1:1791936/1‑140 (MQ=38)
      aGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCAtcctcc                                                                            >  2:3468054/1‑140 (MQ=38)
      aGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCAtcctcc                                                                            >  2:2669379/1‑140 (MQ=38)
      aGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCAtcctcc                                                                            >  2:3355092/1‑140 (MQ=38)
       gTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCa                                                                           <  1:2737371/140‑1 (MQ=37)
        tATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCaa                                                                          >  1:2878415/1‑140 (MQ=37)
          tgtgtgTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGCCATCCTCCAAg                                                                         <  2:1858212/140‑1 (MQ=255)
          tgtgtgTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAgg                                                                        <  1:1524475/140‑1 (MQ=21)
          tgtgtgTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAgg                                                                        <  2:1120136/140‑1 (MQ=21)
          tgtgtgTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAgg                                                                        <  1:3468054/140‑1 (MQ=21)
            tgtgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATTCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCt                                                                      >  1:2597009/1‑140 (MQ=21)
            tgtgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCt                                                                      <  2:340386/140‑1 (MQ=21)
            tgtgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCt                                                                      <  2:1394215/140‑1 (MQ=21)
            tgtgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCt                                                                      >  1:635710/1‑140 (MQ=21)
            tgtgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCt                                                                      >  2:843795/1‑140 (MQ=21)
            tgtgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCt                                                                      >  2:1822161/1‑140 (MQ=21)
             gtgtTAGTGGAAGCGTCTGGCAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTa                                                                     >  1:2307217/1‑140 (MQ=21)
              tgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTaa                                                                    <  1:1822161/140‑1 (MQ=21)
              tgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTaa                                                                    >  2:691499/1‑140 (MQ=21)
              tgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTaa                                                                    <  1:166002/140‑1 (MQ=21)
              tgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTaa                                                                    <  1:1617074/140‑1 (MQ=21)
               gtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTaaa                                                                   <  2:2071417/140‑1 (MQ=21)
                  aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCAtcctc                                                                             >  2:168535/1‑127 (MQ=25)
                  aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAg                                                                         >  2:1597322/1‑131 (MQ=21)
                  aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTaa                                                                    >  1:1064295/1‑136 (MQ=21)
                  aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATAc                                                                >  2:1410859/1‑140 (MQ=21)
                  aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATAc                                                                >  2:1956986/1‑140 (MQ=21)
                  aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATAc                                                                >  1:1286107/1‑140 (MQ=21)
                   gTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACt                                                               <  2:3239987/140‑1 (MQ=21)
                    tGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTc                                                              >  1:1086245/1‑140 (MQ=21)
                      gAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTaa                                                                    >  2:1272845/1‑132 (MQ=21)
                      gAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCct                                                            >  1:858856/1‑140 (MQ=21)
                      gAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCct                                                            >  1:2705820/1‑140 (MQ=21)
                      gAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAACCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCct                                                            >  1:1703056/1‑140 (MQ=21)
                        aGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCctga                                                          <  2:2007154/140‑1 (MQ=21)
                        aGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCctga                                                          <  2:1286107/140‑1 (MQ=21)
                        aGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCctga                                                          >  1:633462/1‑140 (MQ=21)
                        aGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGCGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCct                                                            >  2:2155092/1‑138 (MQ=21)
                        aGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGCGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCct                                                            <  1:2155092/138‑1 (MQ=19)
                          cGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCctgact                                                        >  2:2997677/1‑140 (MQ=21)
                             ttGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCACGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATGCTCCtgactgac                                                     <  2:2307217/139‑1 (MQ=255)
                             cTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCtgactgac                                                     <  1:1961100/140‑1 (MQ=21)
                              tGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCctga                                                          <  1:2018162/134‑1 (MQ=21)
                              tGGAAAGGCGCGCGATACAGGGTAACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGAcc                                                    <  2:2554590/140‑1 (MQ=21)
                                gAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGa                                                  >  2:1677365/1‑140 (MQ=21)
                                 aaaGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGAt                                                 <  1:499636/140‑1 (MQ=21)
                                   aGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAg                                               <  2:277105/140‑1 (MQ=21)
                                   aGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAg                                               <  1:1998717/140‑1 (MQ=21)
                                   aGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAg                                               <  1:1715069/140‑1 (MQ=21)
                                   aGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAg                                               <  1:798001/140‑1 (MQ=21)
                                     gcgcgcGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACt                                                               >  1:3365284/1‑122 (MQ=25)
                                     gcgcgcGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAg                                               <  1:252661/138‑1 (MQ=21)
                                      cgcgcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCtgactgac                                                     <  2:722591/131‑1 (MQ=21)
                                      cgcgcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCtgactgac                                                     >  1:722591/1‑131 (MQ=21)
                                      cgcgcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTg                                             <  2:1363207/139‑1 (MQ=21)
                                      cgcgcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTg                                             >  1:1363207/1‑139 (MQ=21)
                                       gcgcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAg                                               >  2:1486905/1‑136 (MQ=21)
                                       gcgcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGaa                                           <  2:157636/140‑1 (MQ=21)
                                          cgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCa                                        >  2:3320734/1‑140 (MQ=21)
                                           gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAg                                       <  1:1576944/140‑1 (MQ=21)
                                               cAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAg                                               <  2:1876127/128‑1 (MQ=25)
                                                aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCg                                  >  1:3176529/1‑140 (MQ=21)
                                                  ggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTg                                <  1:1971508/140‑1 (MQ=21)
                                                     gACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAgg                             >  1:2836088/1‑140 (MQ=37)
                                                      aCAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAggg                            >  1:2631613/1‑140 (MQ=37)
                                                       cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTg                                             <  1:1098200/122‑1 (MQ=25)
                                                       cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTg                                             >  1:2480753/1‑122 (MQ=37)
                                                       cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTg                                             >  2:1098200/1‑122 (MQ=37)
                                                       cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGt                                      <  1:3287667/129‑1 (MQ=25)
                                                       cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGa                           >  2:1615950/1‑140 (MQ=37)
                                                       cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGa                           >  2:3217917/1‑140 (MQ=37)
                                                       cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGa                           >  2:2842916/1‑140 (MQ=37)
                                                         gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGaaa                         >  1:552210/1‑140 (MQ=37)
                                                          ccccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAg                        >  1:2047455/1‑140 (MQ=37)
                                                          ccccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAg                        >  1:1326701/1‑140 (MQ=37)
                                                            ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGt                                      >  1:3083959/1‑124 (MQ=37)
                                                            ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGc                      <  1:1482446/140‑1 (MQ=21)
                                                            ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGc                      >  1:1022285/1‑140 (MQ=37)
                                                            ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGc                      >  1:1750338/1‑140 (MQ=37)
                                                              gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGa                    >  2:245581/1‑140 (MQ=25)
                                                                 cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGaaaa                 <  1:2578245/140‑1 (MQ=21)
                                                                    aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGaa              <  2:3062502/140‑1 (MQ=21)
                                                                        aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAAcccc          <  2:1750338/140‑1 (MQ=21)
                                                                         tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCg         >  2:2659332/1‑140 (MQ=21)
                                                                         tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCg         >  2:622740/1‑140 (MQ=37)
                                                                           cacaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCAGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGc       <  2:2480821/140‑1 (MQ=21)
                                                                            acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCg      >  1:12598/1‑140 (MQ=21)
                                                                             caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAg    <  2:3092234/140‑1 (MQ=18)
                                                                              atatTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAGCCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAg    <  2:2789489/140‑1 (MQ=18)
                                                                              atatTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAg    >  1:2269610/1‑140 (MQ=18)
                                                                                atTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAggg  <  1:2574460/140‑1 (MQ=14)
                                                                                atTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAggg  <  1:1002481/140‑1 (MQ=14)
                                                                                atTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAggg  <  2:1080717/140‑1 (MQ=14)
                                                                                atTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAggg  >  2:1964802/1‑140 (MQ=14)
                                                                                atTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAggg  <  2:259840/140‑1 (MQ=14)
                                                                                atTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAggg  >  2:2725507/1‑140 (MQ=14)
                                                                                                                                         |                                                                                   
TGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGG  >  NC_000913/3943979‑3944198

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