Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,146,736 A→G 5.3% noncoding (148/169 nt) psrD → novel Novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,146,7360AG 5.3% 339.7 / 6.5 112noncoding (148/169 nt)psrDnovel Novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base A (55/51);  new base G (2/4);  total (57/55)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 113 total aligned reads displayed.

TTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCTGATCAGGT  >  NC_000913/1146588‑1146885
                                                                                                                                                    |                                                                                                                                                     
ttAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTGTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAAc                                                                                                                                                      <  1:1467749/150‑1 (MQ=255)
 tAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAAcc                                                                                                                                                     <  2:872193/150‑1 (MQ=255)
  aGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACca                                                                                                                                                    >  2:774192/1‑150 (MQ=255)
   ggCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACcac                                                                                                                                                   >  2:1265474/1‑150 (MQ=255)
    gCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACcaca                                                                                                                                                  <  1:766636/150‑1 (MQ=255)
     cATATTTTTTTCCATCAGATATAGCGTATTGATCATAGCCATTTTAAACTATGCGATTCGTTTTGCAGGTAGATTATTGTTATCAGCACTGAACGAAAATAAAGCAGTAGCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAg                                                                                                                                                 >  1:88578/1‑150 (MQ=255)
      atatTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAgg                                                                                                                                                >  2:24168/1‑150 (MQ=255)
      atatTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAgg                                                                                                                                                >  2:1554079/1‑150 (MQ=255)
      atatTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAgg                                                                                                                                                >  1:1392394/1‑150 (MQ=255)
       tatTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAGcctgtggt                                                                                                                                               <  1:624436/150‑7 (MQ=255)
       tatTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGc                                                                                                                                               >  1:1651099/1‑150 (MQ=255)
           tttttCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGcc                                                                                                                                           <  1:1453718/150‑1 (MQ=255)
              ttCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCttt                                                                                                                                        <  2:198887/150‑1 (MQ=255)
                ccATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCttttt                                                                                                                                      >  2:1587011/1‑150 (MQ=255)
                ccATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAAATAAAATAAAGCAGTAACCCGCAATGTGTGCTAATTATTGGCAAAAGGCAGCCTGTGGTTGCCttttg                                                                                                                                      <  2:778205/150‑2 (MQ=255)
                 cATCAGATATAGCGTATTGATGAAAGCCATTTTAAACTATGGGCTTCGTTTTGCAGGTTGATGTTTGTTATAAGCACTGAACGAAAATTAAGCAGAAACCCGCAATGTGTGCGAATAATTGGCAAAAGGCAACCACAGGCTGCCTTATTc                                                                                                                                     >  2:1173251/1‑150 (MQ=255)
                   tCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCtt                                                                                                                                   <  1:484728/150‑1 (MQ=255)
                    cAGATATATCGTATTGATGATAGACATTTTATACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTTAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAGCCTGTGGTTGCCttttgccaa                                                                                                                                  <  2:846295/150‑6 (MQ=255)
                     aGATCTACCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTg                                                                                                                                 <  2:681241/150‑1 (MQ=255)
                     aGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTg                                                                                                                                 <  1:593388/150‑1 (MQ=255)
                     aGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTg                                                                                                                                 <  2:882341/150‑1 (MQ=255)
                      gATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGa                                                                                                                                <  2:1373302/150‑1 (MQ=255)
                       atatAGCGTATTGATGATAGCCATTTTAAACTATGCGCGTCGTTTTGCAGGGTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTATCCCGCAATGTGTGCGAATTATTTGCAAAAGGCAACCACAGTCTGCCTTTTTCTTTTAc                                                                                                                               >  1:266460/1‑150 (MQ=255)
                         ataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACATGCTGCCTTTTTCTATGActc                                                                                                                             >  2:305948/1‑150 (MQ=255)
                          taGCGTATTGATGATAGCCATTTTACACTATGCGCTTCGTTTTGCAGGTTGATGCTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCTCAGGCTGCCTTTTTCTTTGActct                                                                                                                            <  2:1196171/150‑1 (MQ=255)
                          taGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGActct                                                                                                                            <  2:650917/150‑1 (MQ=255)
                          taGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAAACACAGGCTGCCTTTTTCTTTGActct                                                                                                                            >  1:1247340/1‑150 (MQ=255)
                           aGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTa                                                                                                                           >  1:719772/1‑150 (MQ=255)
                             cGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATg                                                                                                                         <  1:1623918/150‑1 (MQ=255)
                                       atAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACaaa                                                                                                               <  1:313046/150‑1 (MQ=255)
                                        tAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAg                                                                                                              >  2:639711/1‑150 (MQ=255)
                                         aGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGt                                                                                                             >  2:505292/1‑150 (MQ=255)
                                          gCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGtt                                                                                                            >  1:1327610/1‑150 (MQ=255)
                                          gCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGtt                                                                                                            >  2:1110984/1‑150 (MQ=255)
                                          gCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGtt                                                                                                            >  2:364741/1‑150 (MQ=255)
                                           ccATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACTAAAATAAAGCAGTAACACGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTa                                                                                                           <  1:715582/150‑1 (MQ=255)
                                           ccATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTa                                                                                                           <  1:1074484/150‑1 (MQ=255)
                                             aTTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAat                                                                                                         <  1:364741/150‑1 (MQ=255)
                                              ttttAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAata                                                                                                        <  1:910806/150‑1 (MQ=255)
                                                 tAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgc                                                                                                     >  1:600506/1‑150 (MQ=255)
                                                  aaaCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCTCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgcg                                                                                                    <  2:1265981/150‑1 (MQ=255)
                                                  aaaCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACAACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgcg                                                                                                    >  2:85453/1‑150 (MQ=255)
                                                    aCTACGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgcgcg                                                                                                  <  2:1465807/150‑1 (MQ=255)
                                                       atccgcATCGTTTTGCAGGTAGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGGGAATTAATGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGccc                                                                                               <  2:127013/147‑1 (MQ=255)
                                                       aTGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGccc                                                                                               <  2:1470704/150‑1 (MQ=255)
                                                               cGTTTTGCAGGTTGATGTTTGTTATCAGCAGTGAACGAAAATAAAGCAGTAACCAGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCAATGACGTTACAAAGTTAATATGCGCGCCCTATGCaaa                                                                                       <  1:301898/150‑1 (MQ=255)
                                                                 ttttGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAg                                                                                     <  1:1368298/150‑1 (MQ=255)
                                                                        ggTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCTAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAAt                                                                              <  1:864093/150‑1 (MQ=255)
                                                                        ggTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAAt                                                                              >  1:1224382/1‑150 (MQ=255)
                                                                        ggTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAAt                                                                              >  2:1337049/1‑150 (MQ=255)
                                                                        ggTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAAt                                                                              <  1:408364/150‑1 (MQ=255)
                                                                         gTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAAtt                                                                             >  1:144709/1‑150 (MQ=255)
                                                                            gATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTAcc                                                                          <  2:1651099/150‑1 (MQ=255)
                                                                             aTGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTACCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTAccc                                                                         >  1:1273363/1‑150 (MQ=255)
                                                                              tgtttgttATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTAcccc                                                                        >  2:241186/1‑150 (MQ=255)
                                                                               gtttgttATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCt                                                                       >  2:671104/1‑150 (MQ=255)
                                                                                      aTCAGCACTGAACGAAAATAAAGCAGTAACCCGCATTGTTTTCGAATTAATGGCAAAAGGCAACAACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACtctc                                                                <  2:1395904/150‑1 (MQ=255)
                                                                                      aTCAGCACTGAACGAAAATAAAGCAGTAACACGCAATTGGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGAAGTTACAAACTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACtctc                                                                <  2:286855/150‑1 (MQ=255)
                                                                                       tCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCg                                                               <  2:197017/150‑1 (MQ=255)
                                                                                          gCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCAGTTTTCTTTGACTCTATTACGTTAAAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTAAAACTGAATCTCGagc                                                            >  1:1642072/1‑148 (MQ=255)
                                                                                             cTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCgg                                                         >  1:542942/1‑150 (MQ=255)
                                                                                              tGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGt                                                        <  1:505292/150‑1 (MQ=255)
                                                                                              tGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGAAAATATGCTCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCTATCCGGt                                                        >  1:876824/1‑150 (MQ=255)
                                                                                                aaCGAAAATAAAGCAGTAACCCGATATGTTTGATAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTc                                                      <  1:106658/150‑1 (MQ=255)
                                                                                                  cgatattAAAGCAGTAACCCGCATTGTGTGCGAATTATTGGCAAAAGGCAACGACAGGCTGCGTTTTTCTTTGACTCTATGACGTTAGAAAGTTTATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGCCTCTCGGTCCGGTTCGt                                                    >  2:34705/7‑150 (MQ=255)
                                                                                                  cGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGt                                                    <  1:495819/150‑1 (MQ=255)
                                                                                                        tAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCt                                              >  2:1454743/1‑150 (MQ=255)
                                                                                                         aaaGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTc                                             <  1:1644009/150‑1 (MQ=255)
                                                                                                             cAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCaaaa                                         <  1:1265474/150‑1 (MQ=255)
                                                                                                              aGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCaaaaa                                        >  2:158519/1‑150 (MQ=255)
                                                                                                              aGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCaaaaa                                        >  2:15836/1‑150 (MQ=255)
                                                                                                               gTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAGTTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAAc                                       >  2:781230/1‑150 (MQ=255)
                                                                                                               gTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAAc                                       <  1:879365/150‑1 (MQ=255)
                                                                                                                 aaCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGc                                     >  1:513125/1‑150 (MQ=255)
                                                                                                                  aCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACTTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGcc                                    <  2:1273363/150‑1 (MQ=255)
                                                                                                                   cccGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCt                                   >  2:1319334/1‑150 (MQ=255)
                                                                                                                   cccGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCt                                   >  1:537047/1‑150 (MQ=255)
                                                                                                                   cccGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCt                                   <  2:1392394/150‑1 (MQ=255)
                                                                                                                    ccGCAATGTGTGCGAATTATTGGCAAAAGTCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCtt                                  <  1:881749/150‑1 (MQ=255)
                                                                                                                     cGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTg                                 <  1:1587011/150‑1 (MQ=255)
                                                                                                                            tgtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGTCTCAAAAACGCCTTGATTACCa                          <  2:600506/150‑1 (MQ=255)
                                                                                                                             gtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAg                         <  1:639711/150‑1 (MQ=255)
                                                                                                                             gtgCGAATTAATGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAg                         <  2:124581/150‑1 (MQ=255)
                                                                                                                               gCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAggg                       <  1:1225977/150‑1 (MQ=255)
                                                                                                                                cGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCATGACTCACGATCCGGTTCGTACGGCTCAAAAAGGCCTTGAGTACCAGGGt                      >  1:1554530/1‑150 (MQ=255)
                                                                                                                                   attattTGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACTAAGTTAATATGAGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATc                   <  1:85453/150‑1 (MQ=255)
                                                                                                                                   attattGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATc                   <  1:1500231/150‑1 (MQ=255)
                                                                                                                                   ataattGGCAAAAGGCAACCACAGGCTGGCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATc                   <  1:305948/147‑1 (MQ=255)
                                                                                                                                    ttattGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCt                  <  1:1454743/150‑1 (MQ=255)
                                                                                                                                      caaagaaaaaGGCAGCCTGTGGTTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAAAGCCTTGATTACCAGGGTATCtat                >  2:439821/7‑150 (MQ=255)
                                                                                                                                      attGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCtat                >  1:739719/1‑150 (MQ=255)
                                                                                                                                        tGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATAc              >  1:404193/1‑150 (MQ=255)
                                                                                                                                           cAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATAcccc           >  2:328176/1‑150 (MQ=255)
                                                                                                                                           cAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACAAGGGTATCTATAcccc           >  2:551384/1‑150 (MQ=255)
                                                                                                                                            aaaaGGCAGCCTGTGGTTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCt          >  1:1337432/1‑150 (MQ=255)
                                                                                                                                            aaaaGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCt          >  1:805820/1‑150 (MQ=255)
                                                                                                                                             aaaGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCTg         <  2:513125/150‑1 (MQ=255)
                                                                                                                                                  cAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATTCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCTGATCAg    <  1:1110984/150‑1 (MQ=255)
                                                                                                                                                  cAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCTGATCAg    >  1:627316/1‑150 (MQ=255)
                                                                                                                                                    aCCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAACTTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCTGATCAGGt  <  1:1554079/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCTGATCAGGT  >  NC_000913/1146588‑1146885

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: