Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,727,000 A→G 5.7% Q114Q (CAA→CAG nemA → N‑ethylmaleimide reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,727,0000AG 5.7% 370.2 / 9.0 122Q114Q (CAA→CAGnemAN‑ethylmaleimide reductase
Reads supporting (aligned to +/- strand):  ref base A (54/61);  new base G (4/3);  total (58/64)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

Only 100 of 124 total aligned reads displayed.

CAGGCAAAAGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTCAATGATT  >  NC_000913/1726851‑1727148
                                                                                                                                                     |                                                                                                                                                    
cAGGCAAAAGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCaa                                                                                                                                                      >  2:214958/1‑150 (MQ=255)
 aGGCAAAAGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCCTGGAAAAAAATCAACGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAAc                                                                                                                                                     <  1:1958346/150‑1 (MQ=255)
 aGGCAAAAGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAAc                                                                                                                                                     <  1:1958350/150‑1 (MQ=255)
  ggCAAAAGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAAcc                                                                                                                                                    >  2:1863750/1‑150 (MQ=255)
   gCAAAAGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCt                                                                                                                                                   >  2:828253/1‑150 (MQ=255)
      aaaGGATATGCCGGGGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGc                                                                                                                                                <  2:1693066/150‑1 (MQ=255)
        aGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCgg                                                                                                                                              <  1:1123669/150‑1 (MQ=255)
        aGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCgg                                                                                                                                              <  2:1676717/150‑1 (MQ=255)
        aGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCgg                                                                                                                                              >  1:2069848/1‑150 (MQ=255)
        aGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCgg                                                                                                                                              >  1:1993409/1‑150 (MQ=255)
         ggATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGGTTGCAGGCTGGCgtg                                                                                                                                             <  1:2408962/150‑3 (MQ=255)
          gATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGGTTGCAGGCTGGCgtga                                                                                                                                            <  2:447770/150‑4 (MQ=255)
            tatGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAg                                                                                                                                          <  1:46022/150‑1 (MQ=255)
             atGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGGTTGCAGGCTGGCgtgagaa                                                                                                                                         <  1:10172/150‑7 (MQ=255)
              tGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACTCCAGCCTGCAACCTGGCGGTCAGGc                                                                                                                                        >  1:532036/1‑150 (MQ=255)
              tGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGc                                                                                                                                        >  1:530573/1‑150 (MQ=255)
              tGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCAGCATGCAAAAAAATCACCGCAGGCGTTCATGCTGAAAAAGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGc                                                                                                                                        >  2:719119/1‑150 (MQ=255)
               gCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTTCAGCTGTGGCACAGCGTACTCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCa                                                                                                                                       <  2:1551499/150‑1 (MQ=255)
                 atgTGCGACTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCAcc                                                                                                                                     <  2:2324191/148‑1 (MQ=255)
                 aGGTGCGCCTGGCATCCATAGTCCGGAGCTTATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCAcc                                                                                                                                     >  2:1531302/1‑150 (MQ=255)
                  ggTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCg                                                                                                                                    <  2:1244492/150‑1 (MQ=255)
                  ggTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCg                                                                                                                                    <  1:2475814/150‑1 (MQ=255)
                   gTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCgg                                                                                                                                   >  1:1238431/1‑150 (MQ=255)
                      cgcCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAg                                                                                                                                <  1:1257052/150‑1 (MQ=255)
                        ccTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAgcg                                                                                                                              <  2:1684884/150‑1 (MQ=255)
                         cTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAgcgc                                                                                                                             >  1:743126/1‑150 (MQ=255)
                           ggCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCt                                                                                                                           >  1:1851021/1‑150 (MQ=255)
                           ggCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCt                                                                                                                           >  1:1365340/1‑150 (MQ=255)
                               tCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAg                                                                                                                       <  2:1305916/150‑1 (MQ=255)
                               tCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAg                                                                                                                       <  1:2238772/150‑1 (MQ=255)
                                ccATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGc                                                                                                                      <  1:803174/150‑1 (MQ=255)
                                ccATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGc                                                                                                                      <  1:1942613/150‑1 (MQ=255)
                                    aGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACtt                                                                                                                  <  1:331499/150‑1 (MQ=255)
                                       ccGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAgc                                                                                                               >  2:1234741/1‑150 (MQ=255)
                                       ccGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAgc                                                                                                               <  1:639637/150‑1 (MQ=255)
                                         ggAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAgcgc                                                                                                             <  2:2482494/150‑1 (MQ=255)
                                         ggAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAgcgc                                                                                                             >  2:1554676/1‑150 (MQ=255)
                                         ggAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAgcgc                                                                                                             >  2:2049708/1‑150 (MQ=255)
                                            gCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCggg                                                                                                          <  1:1253578/150‑1 (MQ=255)
                                                  tGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACg                                                                                                    <  2:1118115/150‑1 (MQ=255)
                                                   gCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGt                                                                                                   <  1:1863750/150‑1 (MQ=255)
                                                   gCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGt                                                                                                   <  2:530573/150‑1 (MQ=255)
                                                    ccGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTa                                                                                                  <  1:1164564/150‑1 (MQ=255)
                                                    ccGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTa                                                                                                  <  1:2473505/150‑1 (MQ=255)
                                                     cGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTAc                                                                                                 <  1:2328040/150‑1 (MQ=255)
                                                     cGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTAc                                                                                                 <  1:24229/150‑1 (MQ=255)
                                                     cGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTAc                                                                                                 >  2:136397/1‑150 (MQ=255)
                                                     cGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGGCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTAc                                                                                                 <  1:2151618/150‑1 (MQ=255)
                                                           gAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTctct                                                                                           >  1:2475188/1‑150 (MQ=255)
                                                           gAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTctct                                                                                           <  1:1082436/150‑1 (MQ=255)
                                                           gAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTctct                                                                                           >  1:2303951/1‑150 (MQ=255)
                                                             aaaaaaTCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTgc                                                                                         <  1:645443/150‑1 (MQ=255)
                                                                       cGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGaaaa                                                                               >  1:662242/1‑150 (MQ=255)
                                                                        gCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCTGGAACACGTACTTCTCTGCGCGATGAAAAt                                                                              >  1:1723729/1‑150 (MQ=255)
                                                                                tCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGc                                                                      <  1:678976/150‑1 (MQ=255)
                                                                                 cATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCg                                                                     <  2:729426/150‑1 (MQ=255)
                                                                                     cTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATcc                                                                 >  2:53673/1‑150 (MQ=255)
                                                                                         aaaTGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCgtgt                                                             <  1:2476757/150‑1 (MQ=255)
                                                                                             ggTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGaa                                                         <  2:2321896/150‑1 (MQ=255)
                                                                                              gTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGaaa                                                        >  1:1549695/1‑150 (MQ=255)
                                                                                              gTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGaaa                                                        >  2:2066514/1‑150 (MQ=255)
                                                                                                  tatGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATc                                                    >  1:999006/1‑150 (MQ=255)
                                                                                                  tatGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATc                                                    >  1:1961006/1‑150 (MQ=255)
                                                                                                      gCCGTGCAGCTGTGGCACCCCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCCGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATg                                                >  1:743291/1‑150 (MQ=255)
                                                                                                       ccGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGc                                               <  1:322299/150‑1 (MQ=255)
                                                                                                           gCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCcgcg                                           >  2:1078646/1‑150 (MQ=255)
                                                                                                            cAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGt                                          >  1:535344/1‑150 (MQ=255)
                                                                                                                tgtgGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTgcgc                                      >  1:1883554/1‑150 (MQ=255)
                                                                                                                   ggCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTg                                   <  2:532036/150‑1 (MQ=255)
                                                                                                                    gCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGa                                  >  2:1385941/1‑150 (MQ=255)
                                                                                                                    gCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGa                                  >  2:1243588/1‑150 (MQ=255)
                                                                                                                    gCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGa                                  >  2:1918223/1‑150 (MQ=255)
                                                                                                                     cacaCCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGaa                                 >  1:15918/1‑150 (MQ=255)
                                                                                                                      acacCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAAc                                <  2:1883554/150‑1 (MQ=255)
                                                                                                                       cacCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACt                               <  2:2069848/150‑1 (MQ=255)
                                                                                                                       cacCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACt                               <  2:1098887/150‑1 (MQ=255)
                                                                                                                        acCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTg                              <  1:828253/150‑1 (MQ=255)
                                                                                                                         ccGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTgg                             >  1:479424/1‑150 (MQ=255)
                                                                                                                           ggACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGaa                           <  2:1549695/150‑1 (MQ=255)
                                                                                                                                cATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGAGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGAt                      >  1:1011526/1‑150 (MQ=255)
                                                                                                                                   cctgaccGCCAGGTTGCAGGCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTcc                   >  1:586079/7‑150 (MQ=255)
                                                                                                                                    tctcACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCATGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCa                  <  1:1234741/150‑1 (MQ=255)
                                                                                                                                    tctcACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCa                  >  1:1020577/1‑150 (MQ=255)
                                                                                                                                     ctcGCGCCAGCCTGCAACCTGGCGGTCAGGCACCGATAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAg                 >  2:1976115/1‑150 (MQ=255)
                                                                                                                                     ctcACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCGGTGTTTAAACATCCATGCCGCGTGCGCTTTAACTGGAAGAGATTCCAg                 <  2:252736/150‑1 (MQ=255)
                                                                                                                                     ctcACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAg                 >  2:2182300/1‑150 (MQ=255)
                                                                                                                                      tcACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAgg                <  2:1993409/150‑1 (MQ=255)
                                                                                                                                       cACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAgg                <  2:774316/149‑1 (MQ=255)
                                                                                                                                       cACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAgg                >  1:774316/1‑149 (MQ=255)
                                                                                                                                       cACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGt               >  2:2410629/1‑150 (MQ=255)
                                                                                                                                        aCGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTa              >  2:1752969/1‑150 (MQ=255)
                                                                                                                                         cGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTgcgc                                      <  1:1641281/125‑1 (MQ=255)
                                                                                                                                          gCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGACTCCAGGTATc            <  1:1554676/150‑1 (MQ=255)
                                                                                                                                           ccAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCg           <  2:2303951/150‑1 (MQ=255)
                                                                                                                                             aGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTc         <  2:1365340/150‑1 (MQ=255)
                                                                                                                                                cTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTCAAt      <  1:122294/150‑1 (MQ=255)
                                                                                                                                                 tGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTCAATg     >  1:2119612/1‑150 (MQ=255)
                                                                                                                                                  gCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTCAATGa    >  2:273051/1‑150 (MQ=255)
                                                                                                                                                   cAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTCAATGAt   <  1:1217127/150‑1 (MQ=255)
                                                                                                                                                    aaCCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGTGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTCAATGAtt  <  2:1238431/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGGCAAAAGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTCAATGATT  >  NC_000913/1726851‑1727148

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