Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,392,453 C→A 11.9% intergenic (+18/‑19) glyV → / → glyX tRNA‑Gly/tRNA‑Gly

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,392,4530CA11.9% 237.9 / 2.7 94intergenic (+18/‑19)glyV/glyXtRNA‑Gly/tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base C (21/59);  new base A (3/10);  total (25/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGG  >  NC_000913/4392304‑4392601
                                                                                                                                                     |                                                                                                                                                    
cgcCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGtgctgc                                                                                                                                                      <  2:431282/150‑1 (MQ=255)
  ccTCCCGTAACGACGCAGAAATGCGAAAATAACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGGCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCaa                                                                                                                                                    <  2:1627443/150‑1 (MQ=255)
    tCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAg                                                                                                                                                  <  2:861489/150‑1 (MQ=255)
     cccGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGc                                                                                                                                                 <  2:1844758/150‑1 (MQ=255)
     cccGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGc                                                                                                                                                 <  1:621672/150‑1 (MQ=255)
     cccGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGc                                                                                                                                                 <  2:827149/150‑1 (MQ=255)
       cGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGcac                                                                                                                                               <  2:906705/150‑1 (MQ=255)
         tAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAg                                                                                                                                             >  1:119180/1‑150 (MQ=255)
         tAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAg                                                                                                                                             <  2:171057/150‑1 (MQ=255)
           acgacgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGac                                                                                                                                           <  2:1545570/150‑1 (MQ=255)
           acgacgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGa                                                                                                                                            <  2:420798/149‑1 (MQ=255)
           acgacgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGa                                                                                                                                            >  1:420798/1‑149 (MQ=255)
            cgacgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGacc                                                                                                                                          <  1:99376/150‑1 (MQ=255)
              acgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGaccac                                                                                                                                        <  2:603959/150‑1 (MQ=255)
               cgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGaccacc                                                                                                                                       <  2:1158915/150‑1 (MQ=255)
                gCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTAGCTAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCAccc                                                                                                                                      <  2:258842/150‑1 (MQ=255)
                 cAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAGTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGTGGTCGCGACTTCGAGTCTCGGTTCACGCGCCAAAAGTAGAAACCTGGTGGAAAGCACTGagcaccaa                                                                                                                                     >  1:840704/1‑148 (MQ=18)
                  aGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGTTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAACGTGCTGCAGGGCCCAGCCCCCCCac                                                                                                                                    >  1:1497849/1‑149 (MQ=18)
                   gAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAATTCATACAAAGCAAACACAAACCaat                                                                                                                                   >  1:937304/1‑149 (MQ=18)
                   gAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTAGAGTCTAGTTTCCCGCTCCAAAATTTGAAAAGTGCTGAAAAGCACAGAAAACCCAAg                                                                                                                                   >  1:1581129/1‑150 (MQ=21)
                       tGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCggg                                                                                                                               <  2:1698819/150‑1 (MQ=255)
                       tGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCggg                                                                                                                               <  2:1790412/150‑1 (MQ=255)
                       tGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCggg                                                                                                                               <  1:239325/150‑1 (MQ=255)
                        gCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCTGTTGGTAGAGCACGACCTTGCCAAGTTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGATCACCCAAGCGTGa                                                                                                                              <  2:870167/150‑1 (MQ=255)
                        gCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGa                                                                                                                              <  1:173780/150‑1 (MQ=255)
                        gCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGa                                                                                                                              <  2:1984869/150‑1 (MQ=255)
                        gCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGa                                                                                                                              <  2:1512943/150‑1 (MQ=255)
                        gCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGa                                                                                                                              <  2:1344154/150‑1 (MQ=255)
                         cGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGGCTCGTTTCCCTCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGaa                                                                                                                             >  1:392648/1‑150 (MQ=255)
                          gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAAt                                                                                                                            <  2:119180/150‑1 (MQ=255)
                          gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAAt                                                                                                                            <  2:1375544/150‑1 (MQ=255)
                            aaaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTACAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAg                                                                                                                          <  1:293994/150‑1 (MQ=255)
                             aaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGATGCAAAGCACAGACAAACCAAGCGGGAATAGc                                                                                                                         >  1:1593/1‑150 (MQ=255)
                             aaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGGCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGc                                                                                                                         <  1:451766/150‑1 (MQ=255)
                               ttACGAAAGCAAACTTAAGAAGTACGAGGGAATATCTCAGTTGGTTGAGCACGTCCTTGCCAAGGGCGGGGTCGCGAGTTCGAGTCTCGTTTCCGGCTCCAAAATTTGAAAAGTGCTGCAAAGCACTGACCACCCAAGCGGGAATAGCTc                                                                                                                       <  2:1301634/150‑1 (MQ=17)
                                 aCGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCACTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAg                                                                                                                     <  2:547007/150‑1 (MQ=255)
                                  cGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGac                                                                                                                                           >  2:1706452/1‑127 (MQ=255)
                                  cGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGac                                                                                                                                           <  1:1706452/127‑1 (MQ=255)
                                   gAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAATTGATTAAAGGCAAAGACCACCAAACCGGGATTAGcaaaggt                                                                                                                   >  2:1902897/1‑144 (MQ=17)
                                    aaaGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTg                                                                                                                  <  1:165884/150‑1 (MQ=255)
                                     atggaaaaTTATGTAGTTCGCGGGAATATCTCAGTTGGTAGAGCACGACCTTTCCTAGGTCTTGGTCGCGAGTTCGTGTCTCGTTTCCCTCTCCAAAATTTGAAAAGTGCTCCACAGCACAGACCACCCAAGCGGGAATAGCTCAGTTgg                                                                                                                 <  2:1043858/146‑1 (MQ=14)
                                     aaGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTgg                                                                                                                 <  2:1964236/150‑1 (MQ=255)
                                      aGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGt                                                                                                                <  1:702996/150‑1 (MQ=255)
                                       gCAAAATTAAGTAGTACTCGGGAATAGCTCAGTTTGTAGAGCTCGACCTTTCCCACTTCGGGGTCGCCTGTTGGAGTCTCGTTTCCCGCGCCGAAATTTGAAAAGTGCTGCAAACCACTGAACACCCAAGCGGGAATAGCTCAGTTGGTa                                                                                                               <  2:1256407/150‑1 (MQ=17)
                                        cAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTag                                                                                                              <  1:245841/150‑1 (MQ=255)
                                        cAAAATTAAGTAGTAAGCTGGAATAGCTCAGTTGTTAGTGCACGTCCTTCCGAAGGTCGTTGTCGCTATTTGGAGGCTCCGTTACCTCTCCAAAATTTGAAAATTTCGTCATATCACAGACGACCCAAGCGGGAATAGCTCAGTTGGTag                                                                                                              <  2:1196251/150‑1 (MQ=9)
                                        cAAAAATACGTAGAACGCGGTAATACCTCAGTTGTTAGAGCACGACCTTGCCAAGCTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTCCTGCAAACCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTag                                                                                                              <  1:1708964/150‑1 (MQ=14)
                                         aaaaTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGAGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGATGCAAAGCACAGACCACCCAAGCGGGAATAGCTAAGTTGGTaga                                                                                                             >  2:923688/1‑150 (MQ=255)
                                           aaTTAAGTAGTACGCGCGAATAGCTCAGTTTGTAGAGCACGACCTTTCCAAGGTCGTGGTCGCTATTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTCCTGAAAATCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGc                                                                                                           <  2:346187/150‑1 (MQ=14)
                                            aTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCa                                                                                                          <  1:1331672/150‑1 (MQ=255)
                                                     taCGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGc                                                                                                 <  1:564267/150‑1 (MQ=255)
                                                               tAGCTCAGTTGCTCGAGAACTACCTTGACAAGTTCGGTGTCGCGAGTTCTAGTCTCGTTTCCCTCTCATAAATTTGAAAAGTGCTCCAAACCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                       <  1:1408972/150‑1 (MQ=14)
                                                                aGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgtgc                                                                                      <  1:496340/150‑4 (MQ=255)
                                                                   taaGTTGGTAGTGGAGGAGCTTTCCTAGGGCGTGGTCGTTATTTCGAGGCTCGGTTCCCTCTCAAAAATTTGAAAAGTTCTGAACATAACAGAACACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcg                                                                                   <  2:802949/148‑1 (MQ=2)
                                                                   tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCgg                                                                                                                                <  1:1305075/105‑1 (MQ=255)
                                                                   tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTAGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCgg                                                                                                                                >  2:1305075/1‑105 (MQ=255)
                                                                     agtttttAGAGCACGACCTTGCCAAGTTCGGTGTCGCGATTTCGACTCTCGTTTCCCTCTCATAAATTTGAAAATTCCTTCACATAACATAACACCCAAGCGGGAATAGCTCAGTTGGTAGATCACGACCTTGCCAAGGTCGGGGTcgcg                                                                                 <  2:431640/144‑1 (MQ=14)
                                                                       ttGGTAGAGCTCGTCCTTTCCCATGGCGGGGTCGCGATTTCGCGTCTCGGTTCCCTCGCCAAAATTTGAAAAGTTCTCCAAAACACATACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAg                                                                               <  2:516864/150‑1 (MQ=14)
                                                                           tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                       <  2:1412473/138‑1 (MQ=40)
                                                                           tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                       >  1:1412473/1‑138 (MQ=40)
                                                                           tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCg                                                                           <  2:1517229/150‑1 (MQ=39)
                                                                           tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCg                                                                           >  1:1517229/1‑150 (MQ=39)
                                                                            agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGtt                                                                                                                   <  2:262465/109‑1 (MQ=255)
                                                                            agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGtt                                                                                                                   >  1:262465/1‑109 (MQ=255)
                                                                             gagCCCGACCTGGCCAAGGGCGGGGTCGCGAGTTCGAGGCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCTCCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAg                                                                         <  2:156932/150‑1 (MQ=14)
                                                                             gagCACGACCTTGCCAAGGTCGGGGTCGCGATTTCGCGTCTCGGTTCCCGCTCCAAAATTTGAAAAGTGCTTCAAATCACAGAACACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAg                                                                         <  2:287569/150‑1 (MQ=14)
                                                                              accccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTTCAAAACACAGACCACAAAAGAGGGAATAGCGCAGTTGGTAGAGCACGACCTTGCCAAGGTCAGGGTAGAGAGTTCGAGt                                                                        >  2:428147/5‑150 (MQ=14)
                                                                                             aGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAACCACAGACCACCCAAGCGGGAATAGCTCAGTCGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTcc                                                         <  1:713215/150‑1 (MQ=255)
                                                                                                        gcgATTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGaaa                                              <  1:1564184/150‑1 (MQ=255)
                                                                                                         cgATTTCGAGTCTCGTTTCCCTCTCGAAAATTTTATAAGTTCTTAAAATAACAGACCAAGCAAGCTGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAg                                             <  2:1593/150‑1 (MQ=14)
                                                                                                           cgtgcGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTCCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTg                                           <  1:428147/146‑1 (MQ=255)
                                                                                                            tatCGAGTCCCGTTACCGCCTCAAAAATTAGAAAAGTCCTTCACATCAATGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGc                                          <  2:147115/148‑1 (MQ=17)
                                                                                                             ttCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGAAAAGCACAGCCCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTACGAGTCTCGTTACCCGCTCCAAAATTTAAAAGAGCt                                         >  2:5255/1‑150 (MQ=14)
                                                                                                               ggtgtgtcGGTTCCGTCTCAAAAATTGGAAAAGTGCTGAAAATCACAGACGTAGCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGt                                       <  1:1973419/144‑1 (MQ=17)
                                                                                                               cGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTGGAGTCTCGATGCGCGCTCCAAAATTAGCTAGTGCTGt                                       >  1:1355572/1‑150 (MQ=14)
                                                                                                                gcgcctcGGTTCCCGCTCCTAAATTTGAAAATTTCGGAATAGCACTTAAGTCCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTa                                      <  1:662877/146‑1 (MQ=17)
                                                                                                                 agcctcGTTTCCCCCTCATAAATTTGAAAAGTGCTGCAAATCACATAACACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTaa                                     <  1:1210926/147‑1 (MQ=17)
                                                                                                                      ccgttCCCGCTCCTAAATTTGAAAAGTTCTTCATATCACATACCTACAAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGcac                                <  2:597248/147‑1 (MQ=17)
                                                                                                                        tttCCGTCTCCTAAATTTGAAAATTTCTGAACATTAAATACGACGGTAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGTCACAg                              <  2:913110/150‑1 (MQ=14)
                                                                                                                           cccGCTCAAAAATTTGAAAAGTCCTCCACAGCACATACCACGCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGacc                           <  2:326648/150‑1 (MQ=18)
                                                                                                                            ccGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGACGCGAGATCGGCTCTCGTGTGCAGAGCAGAAAGTTGGAAGTGGAGTAAGGCAGAGacgc                          >  1:637636/1‑148 (MQ=14)
                                                                                                                              cctcataaaTTTGAAAAGTCCGCAAAATAAAATACCACGCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGaccacc                        <  2:346079/144‑1 (MQ=14)
                                                                                                                               cTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATTCAAGAGATCGAAGGCGCAcaccacac                       >  1:631606/1‑148 (MQ=14)
                                                                                                                                tCCAAAATTTGAAAAGTGCTCAAAATCACATACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCa                      <  2:966920/150‑1 (MQ=255)
                                                                                                                                 cctaaaTTTGAAAATTGCTGAACATAAATGACGACGCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCaa                     <  1:1257437/147‑1 (MQ=14)
                                                                                                                                  cAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTTTGAAAGTGATATAATGCCCACACCACCCAAg                    >  1:519639/1‑150 (MQ=17)
                                                                                                                                  cAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATTAAAGTGCTGTAAGGGACAGAGCcccccac                    >  1:445644/1‑147 (MQ=17)
                                                                                                                                   aaaaTTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGATGTAAGGACCACACCACCCAAGc                   >  2:446443/1‑150 (MQ=255)
                                                                                                                                      atgtGAAAAGTCCTGAACAACACATACGTCCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCggg                <  1:623678/147‑1 (MQ=14)
                                                                                                                                         tGAAAATTGCTTCAAATAACAGACGTCGTACGGGGGAATAGCTCAGTTGGTAGAGCAGGACCTTGCAAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAAt             <  2:257665/150‑1 (MQ=11)
                                                                                                                                         tGAAAATTGCTCCACATAACATACGACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAAt             <  1:196358/150‑1 (MQ=14)
                                                                                                                                          gAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGc                                                                                                 <  2:1732619/65‑1 (MQ=255)
                                                                                                                                          gAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGc                                                                                                 >  1:1732619/1‑65 (MQ=255)
                                                                                                                                             aaGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgggg                                                                                      >  2:581764/1‑73 (MQ=39)
                                                                                                                                             aaGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgggg                                                                                      <  1:581764/73‑1 (MQ=39)
                                                                                                                                                    gcAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATTAAAGAAATTTAATGAACACACCACCCAAGAGCGAATATCGAAGTTgt  >  2:547651/1‑149 (MQ=14)
                                                                                                                                                     |                                                                                                                                                    
CGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGG  >  NC_000913/4392304‑4392601

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