Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,716,061 C→A 5.1% G388C (GGT→TGT)  tyrS ← tyrosine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,716,0610CA 5.1% 261.0 / 2.9 78G388C (GGT→TGT) tyrStyrosine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (32/42);  new base A (2/2);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

GGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAG  >  NC_000913/1715912‑1716206
                                                                                                                                                     |                                                                                                                                                 
ggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCAcc                                                                                                                                                   <  2:1333788/150‑1 (MQ=255)
 gAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCg                                                                                                                                                  >  2:1017922/1‑150 (MQ=255)
   gTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGtt                                                                                                                                                >  1:339184/1‑150 (MQ=255)
   gTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGtt                                                                                                                                                >  2:471988/1‑150 (MQ=255)
     gagaaatcACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACAGTTaa                                                                                                                                              >  1:954377/7‑150 (MQ=255)
     gagaaatcACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTAGTCAGGATCGGACTGTTTTTCACAGTTaa                                                                                                                                              >  1:948730/7‑150 (MQ=255)
     gATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTaa                                                                                                                                              >  2:881619/1‑150 (MQ=255)
     gATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTaa                                                                                                                                              >  2:649287/1‑150 (MQ=255)
      aTTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAAt                                                                                                                                             <  2:1274845/150‑1 (MQ=255)
       tttCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTATAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATg                                                                                                                                            <  1:1378981/150‑1 (MQ=255)
       tttCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATg                                                                                                                                            <  1:794156/150‑1 (MQ=255)
          ctcACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTg                                                                                                                                         >  2:82856/1‑150 (MQ=255)
           tcACTCCCCCTTTCCACTTAATGCATTATTTCTAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGa                                                                                                                                        <  1:365306/150‑1 (MQ=255)
            cACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCTCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGAt                                                                                                                                       <  1:571181/150‑1 (MQ=255)
               tCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGc                                                                                                                                    <  2:1427344/150‑1 (MQ=255)
               tCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGc                                                                                                                                    >  1:280611/1‑150 (MQ=255)
                  cccTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCATTAACGGTGaaaaacag                                                                                                                                 <  1:1535982/150‑8 (MQ=255)
                    cTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTgg                                                                                                                               <  1:553262/150‑1 (MQ=255)
                               aTGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGtt                                                                                                                    <  1:1361167/150‑1 (MQ=255)
                                tGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTg                                                                                                                                >  1:1162619/1‑137 (MQ=255)
                                tGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTg                                                                                                                                <  2:1162619/137‑1 (MQ=255)
                                 gCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGTCGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGtttt                                                                                                                  <  2:1285611/150‑1 (MQ=255)
                                  cATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTa                                                                                                                 <  2:222850/150‑1 (MQ=255)
                                    ttattTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCTAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACg                                                                                                               <  1:1222206/150‑1 (MQ=255)
                                    ttattTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACg                                                                                                               <  1:711825/150‑1 (MQ=255)
                                       tttCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGc                                                                                                            >  2:1102346/1‑150 (MQ=255)
                                          ccAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGGCTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTg                                                                                                         >  1:637117/1‑150 (MQ=255)
                                           cAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGa                                                                                                        >  2:1344625/1‑150 (MQ=255)
                                           cAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGa                                                                                                        >  2:1213606/1‑150 (MQ=255)
                                                 aTCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACgg                                                                                                  <  1:1421031/150‑1 (MQ=255)
                                                  tCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACggg                                                                                                 >  2:587889/1‑150 (MQ=255)
                                                     gACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAg                                                                                              <  1:768424/150‑1 (MQ=255)
                                                       cAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGt                                                                                            <  1:183710/150‑1 (MQ=255)
                                                        aGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGtt                                                                                           <  2:339184/150‑1 (MQ=255)
                                                         gTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTg                                                                                          <  1:350169/150‑1 (MQ=255)
                                                           aaTTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCACACAGACTATCTTGTTCTTTAAAGTAGTATTCAGGATCGGACTGTTTTTCAGCGTTAATGGTGATGGCATTTGAGGCGATAGTTTTAGGTGCCTGACCACGGGAAGGTTGCa                                                                                        <  2:988760/150‑1 (MQ=255)
                                                           aaTTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCTGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCa                                                                                        <  2:1076825/150‑1 (MQ=255)
                                                            aTTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAg                                                                                       <  1:1332652/150‑1 (MQ=255)
                                                            aTTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAg                                                                                       <  2:277403/150‑1 (MQ=255)
                                                            aTTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCAGTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAg                                                                                       <  1:82856/150‑1 (MQ=255)
                                                               cTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTg                                                                                                                                         <  1:1453433/97‑1 (MQ=255)
                                                               cTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTg                                                                                                                                         >  2:1453433/1‑97 (MQ=255)
                                                                 ttttACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAg                                                                                  <  2:1481466/150‑1 (MQ=255)
                                                                 ttttACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAg                                                                                  <  1:1017922/150‑1 (MQ=255)
                                                                     aCCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGTCTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATc                                                                              <  2:637117/150‑1 (MQ=255)
                                                                      ccGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCg                                                                             >  1:635157/1‑150 (MQ=255)
                                                                       cgcgACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGtttt                                                                                                                  >  1:1412476/1‑112 (MQ=255)
                                                                       cgcgACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGtttt                                                                                                                  <  2:1412476/112‑1 (MQ=255)
                                                                       cgcgACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGa                                                                            >  2:1301727/1‑150 (MQ=255)
                                                                        gcgACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGAc                                                                           <  1:1373596/150‑1 (MQ=255)
                                                                           aCGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAg                                                                        >  2:1500459/1‑150 (MQ=255)
                                                                            cGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGt                                                                       >  1:1053610/1‑150 (MQ=255)
                                                                                 tAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTg                                                                  <  2:100873/150‑1 (MQ=255)
                                                                                       aaaaCGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAg                                                            >  2:48701/1‑150 (MQ=255)
                                                                                       aaaaCGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTAAGGATCGGACTGTTTTTCACCGTTTATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTAAGAATCGACCAGTGCCTGCATCAg                                                            >  2:1314321/1‑150 (MQ=255)
                                                                                        aaaCGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAgg                                                           <  1:1500459/150‑1 (MQ=255)
                                                                                        aaaCGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAgg                                                           >  1:851740/1‑150 (MQ=255)
                                                                                            gACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTg                                                       <  2:5038/150‑1 (MQ=255)
                                                                                                         aTCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCAt                                          >  2:1401145/1‑150 (MQ=255)
                                                                                                          tcttcttctTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACAGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCAtc                                         <  2:948730/150‑1 (MQ=255)
                                                                                                          tcttcttctTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACAGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCAtc                                         <  2:954377/150‑1 (MQ=255)
                                                                                                           cttcttcttTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCAtct                                        <  1:881619/150‑1 (MQ=255)
                                                                                                             tcttcttTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTg                                                                                          <  2:234979/98‑1 (MQ=255)
                                                                                                             tcttcttTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTg                                                                                          >  1:234979/1‑98 (MQ=255)
                                                                                                                 cttTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGATAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCa                                  >  1:897646/1‑150 (MQ=255)
                                                                                                                     aaaGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCgg                              >  2:799565/1‑150 (MQ=255)
                                                                                                                      aagaagTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGt                             >  2:576134/1‑150 (MQ=255)
                                                                                                                               tCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTcc                    <  1:113361/150‑1 (MQ=255)
                                                                                                                               tCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTcc                    <  1:471988/150‑1 (MQ=255)
                                                                                                                                cAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCt                   >  1:808356/1‑150 (MQ=255)
                                                                                                                                  ggATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTgc                 >  1:377396/1‑150 (MQ=255)
                                                                                                                                   gATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTgcg                <  2:634541/150‑1 (MQ=255)
                                                                                                                                    aTCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGACCTTTTCCATCTCAACCATCGGTACGCCGTCCTgcgc               >  2:85805/1‑150 (MQ=255)
                                                                                                                                        gACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGc           <  1:799565/150‑1 (MQ=255)
                                                                                                                                             tttttCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTc      <  2:1064636/150‑1 (MQ=255)
                                                                                                                                             tttttCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTc      >  1:496036/1‑150 (MQ=255)
                                                                                                                                              ttttCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCg     >  2:463423/1‑150 (MQ=255)
                                                                                                                                               tttCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGa    <  1:719411/150‑1 (MQ=255)
                                                                                                                                                 tCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAg  <  1:463423/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAG  >  NC_000913/1715912‑1716206

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: ‑