Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,815,819 C→G 5.5% intergenic (‑51/+44) pyrE ← / ← rph orotate phosphoribosyltransferase/truncated RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,815,8190CG 5.5% 304.4 / 7.0 91intergenic (‑51/+44)pyrE/rphorotate phosphoribosyltransferase/truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base C (35/49);  new base G (3/2);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

AGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCC  >  NC_000913/3815672‑3815967
                                                                                                                                                    |                                                                                                                                                    
aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAg                                                                                                                                                    >  2:434253/1‑150 (MQ=255)
    gCTTTTGCGCCCGTATTTAAGCGTAAACTCGCAAAACTTTAACACCTGCTTGCTAAGCGAAAATTCAATAAACTGGCGCTTATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAAAAAAAAAAGGCGACTCATCAGTCGc                                                                                                                                                <  2:429965/150‑1 (MQ=255)
     cTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGcc                                                                                                                                               <  2:189272/150‑1 (MQ=255)
     cTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGcc                                                                                                                                               <  1:643782/150‑1 (MQ=255)
       tttGCGCCGGTTTTTCTGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCtt                                                                                                                                             <  1:343517/150‑1 (MQ=255)
          gcgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaa                                                                                                                                          <  1:1057744/150‑1 (MQ=255)
          gcgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAGGCGACTCATCAGTCGCCTTaaaa                                                                                                                                         <  2:1421504/150‑1 (MQ=255)
           cgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaa                                                                                                                                         >  2:864857/1‑150 (MQ=255)
           cgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaa                                                                                                                                         >  1:922256/1‑150 (MQ=255)
             cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCtttttttt                                                                                                                                       <  2:477663/150‑7 (MQ=255)
             cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAt                                                                                                                                       >  1:1202372/1‑150 (MQ=255)
             cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAt                                                                                                                                       >  1:1604079/1‑150 (MQ=255)
              ccGGATTTCAGCGTAAACTCGCCAAACTTTAACCCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttttg                                                                                                                                      <  1:1148074/150‑8 (MQ=255)
               cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCa                                                                                                                                     >  2:991096/1‑150 (MQ=255)
                ggATTTCAGCGTAAACTCGGCAATCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGGCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAg                                                                                                                                    <  2:838323/150‑1 (MQ=255)
                ggATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAg                                                                                                                                    <  1:425821/150‑1 (MQ=255)
                  aTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGtt                                                                                                                                  >  1:1234994/1‑150 (MQ=255)
                      cAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                                                                              <  2:1202372/150‑1 (MQ=255)
                               cTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCC                                                                                                                     <  2:1652824‑M1/150‑1 (MQ=255)
                                 cGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTT                                                                                                                   <  2:801288‑M1/150‑1 (MQ=255)
                                     aaaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                                                                               <  1:850023‑M1/150‑1 (MQ=255)
                                     aaaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                                                                               <  2:1204854‑M1/150‑1 (MQ=255)
                                     aaaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                                                                               <  2:1207525‑M1/150‑1 (MQ=255)
                                     aaaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                                                                               <  2:1234994‑M1/150‑1 (MQ=255)
                                       aCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcg                                                                                                             <  2:659799‑M1/150‑3 (MQ=255)
                                       aCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcg                                                                                                             <  2:744130‑M1/150‑3 (MQ=255)
                                        cTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgt                                                                                                            >  2:146929‑M1/1‑147 (MQ=255)
                                        cTTTAACACCTGATTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCATTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGAGCCTCATAAGTCGCCTTAAAAATAAGTTTGCCACCGCACCCTTCTGCctt                                                                                                            >  2:447880‑M1/1‑147 (MQ=255)
                                          ttAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTATACAGACAAAAAAAAGGCGAATCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcc                                                                                                          >  2:1541045‑M1/1‑145 (MQ=255)
                                             aaacatgcttgctaagcgCAAATTCAATAAACTGGCGCTGATATGGTTTCATGACTTCGCTCCTCATCTTAATTTTATAAACACAAAAAAAAGGCGAATCATAAGTAGCCTTAAAAATAAGTTTGCAAGCGCCGCCTTCTGCcgtcccct                                                                                                       >  2:431912‑M1/6‑142 (MQ=255)
                                              cacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctg                                                                                                      >  1:1469970‑M1/1‑141 (MQ=255)
                                                ccTGCTTGCTAAGCGCAAATTCAATAATCTGGCGCTGATTTGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCGAGCGCCGCCTTCTGCcgtccccttca                                                                                                    <  2:238930‑M1/150‑12 (MQ=255)
                                                ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                                                                    <  2:1396764‑M1/150‑12 (MQ=255)
                                                       gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                                             >  1:319379‑M1/1‑132 (MQ=255)
                                                       gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                                             >  1:860433‑M1/1‑132 (MQ=255)
                                                          aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatga                                                                                           <  2:1296321‑M1/150‑21 (MQ=255)
                                                           aGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                                                                         <  1:146929‑M1/150‑23 (MQ=255)
                                                           aGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                                                                         >  1:1368070‑M1/1‑128 (MQ=255)
                                                               cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                                     >  1:32823‑M1/1‑124 (MQ=255)
                                                               cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                                     <  1:864857‑M1/150‑27 (MQ=255)
                                                                     cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCCGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtctt                                                                               >  1:1572889‑M1/1‑118 (MQ=255)
                                                                      aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                                                                              <  2:104585‑M1/150‑34 (MQ=255)
                                                                      aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                                                                              <  1:113761‑M1/150‑34 (MQ=255)
                                                                      aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                                                                              <  1:1359369‑M1/150‑34 (MQ=255)
                                                                        tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcgg                                                                            <  2:922256‑M1/150‑36 (MQ=255)
                                                                        tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcgg                                                                            <  1:1495164‑M1/150‑36 (MQ=255)
                                                                        tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcgg                                                                            <  2:635492‑M1/150‑36 (MQ=255)
                                                                         aaaCTGGGGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTAGTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggt                                                                           <  1:371377‑M1/150‑37 (MQ=255)
                                                                         aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGTGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggt                                                                           <  1:1130833‑M1/150‑37 (MQ=255)
                                                                          aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                                                                          <  2:319379‑M1/150‑38 (MQ=255)
                                                                          aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCCCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtacactgcacttcaatgatgcgcccgtcttcggta                                                                          >  2:1245646‑M1/1‑113 (MQ=255)
                                                                           aCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtca                                                                         <  1:743109‑M1/150‑39 (MQ=255)
                                                                                 gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcacta                                                                   >  2:180468‑M1/1‑106 (MQ=255)
                                                                                   tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacg                                                                 <  1:596046‑M1/150‑47 (MQ=255)
                                                                                     atatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt                                                               >  2:1180165‑M1/1‑102 (MQ=255)
                                                                                      tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc                                                              >  2:1247927‑M1/1‑101 (MQ=255)
                                                                                      tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc                                                              <  1:1215834‑M1/150‑50 (MQ=255)
                                                                                      tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc                                                              <  2:1469970‑M1/150‑50 (MQ=255)
                                                                                       atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactaggttca                                                             >  2:556827‑M1/1‑100 (MQ=255)
                                                                                         ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatg                                                           >  1:1020893‑M1/1‑98 (MQ=255)
                                                                                         ggTTTCATGCCTTCGCTCCTCATCCTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGt                                                                                                                                                   <  1:390348/62‑1 (MQ=255)
                                                                                         ggTTTCATGCCTTCGCTCCTCATCCTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGt                                                                                                                                                   >  2:390348/1‑62 (MQ=255)
                                                                                              cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCCGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt                                                      >  2:28348‑M1/1‑93 (MQ=255)
                                                                                                tGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtct                                                    >  2:1516547‑M1/1‑91 (MQ=255)
                                                                                                 gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctc                                                   <  1:1180165‑M1/150‑61 (MQ=255)
                                                                                                       gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggc                                             <  2:170938‑M1/150‑67 (MQ=255)
                                                                                                       gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggc                                             <  1:441909‑M1/150‑67 (MQ=255)
                                                                                                       gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggc                                             <  2:1604079‑M1/150‑67 (MQ=255)
                                                                                                        ctcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggca                                            <  2:705610‑M1/150‑68 (MQ=255)
                                                                                                        ctcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggca                                            <  1:502978‑M1/150‑68 (MQ=255)
                                                                                                           ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagag                                         <  1:434253‑M1/150‑71 (MQ=255)
                                                                                                           ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagag                                         <  1:1636844‑M1/150‑71 (MQ=255)
                                                                                                            tcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagt                                        >  1:1096038‑M1/1‑79 (MQ=255)
                                                                                                             cATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtc                                       >  2:412961‑M1/1‑78 (MQ=255)
                                                                                                                 ttACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttca                                   >  1:1088879‑M1/1‑74 (MQ=255)
                                                                                                                    cTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacg                                <  2:899024‑M1/150‑80 (MQ=255)
                                                                                                                    cTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacg                                >  2:29347‑M1/1‑71 (MQ=255)
                                                                                                                    cTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacg                                >  2:608074‑M1/1‑71 (MQ=255)
                                                                                                                    cTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagag                                         <  1:736037‑M1/141‑71 (MQ=255)
                                                                                                                    cTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagag                                         >  2:736037‑M1/1‑71 (MQ=255)
                                                                                                                          tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcc                          >  2:56423‑M1/1‑65 (MQ=255)
                                                                                                                          tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcc                          <  1:56423‑M1/150‑86 (MQ=255)
                                                                                                                            cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccag                        >  2:633666‑M1/1‑63 (MQ=255)
                                                                                                                             aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccaga                       <  1:217473‑M1/150‑89 (MQ=255)
                                                                                                                              gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagat                      >  1:75987‑M1/1‑61 (MQ=255)
                                                                                                                              gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagat                      >  1:1309826‑M1/1‑61 (MQ=255)
                                                                                                                               aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatc                     <  1:1491692‑M1/150‑91 (MQ=255)
                                                                                                                                 gatttttaaggcgactGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcg                    >  1:1261561‑M1/8‑59 (MQ=255)
                                                                                                                                   aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgca                  >  2:1582653‑M1/1‑57 (MQ=255)
                                                                                                                                        aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccg             >  2:1644927‑M1/1‑52 (MQ=255)
                                                                                                                                         aGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgc            >  2:1262610‑M1/1‑51 (MQ=255)
                                                                                                                                                tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccg     >  1:431212‑M1/1‑44 (MQ=255)
                                                                                                                                                gcatcagTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTCCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccg     <  2:1467568‑M1/149‑107 (MQ=255)
                                                                                                                                                   tcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttc  <  1:180468‑M1/150‑110 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCC  >  NC_000913/3815672‑3815967

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: