Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,102,340 (C)6→7 coding (513/720 nt) trmI ← tRNA m(7)G46 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,102,3341.C98.7% 277.7 / ‑4.3 78coding (519/720 nt)trmItRNA m(7)G46 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (35/42);  total (35/43)
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

TTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGC  >  NC_000913/3102190‑3102478
                                                                                                                                                 |                                                                                                                                                 
ttgttCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑‑cccccc                                                                                                                                            >  2:313959/1‑151 (MQ=255)
   ttCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑‑ccccccc                                                                                                                                           >  2:1476330/1‑148 (MQ=255)
      aaaTTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTg                                                                                                                                       <  1:1577732/151‑1 (MQ=255)
      aaaTTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTg                                                                                                                                       >  2:1572185/1‑151 (MQ=255)
          ttCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGt                                                                                                                                   >  1:840741/1‑151 (MQ=255)
            cGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGttt                                                                                                                                 <  2:1347480/151‑1 (MQ=255)
            cGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGtt                                                                                                                                  <  1:1685957/150‑1 (MQ=255)
                aCCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCtt                                                                                                                             >  1:1402542/1‑151 (MQ=255)
                aCCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCtt                                                                                                                             <  2:717766/151‑1 (MQ=255)
                  cGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCtttt                                                                                                                           <  1:1547027/151‑1 (MQ=255)
                       gTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCa                                                                                                                      <  2:202568/151‑1 (MQ=255)
                         gATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCTGTCGCCATATGGAATACG‑‑CCCCCCAGCTGCAGTTTGCTTTTTGCCAGt                                                                                                                    <  1:850151/150‑1 (MQ=255)
                         gATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAg                                                                                                                     <  1:1379346/150‑1 (MQ=255)
                         gATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGt                                                                                                                    >  1:1365514/1‑151 (MQ=255)
                         gATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGt                                                                                                                    >  2:1074935/1‑151 (MQ=255)
                                 aCGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCa                                                                                                              >  1:1257798/1‑150 (MQ=255)
                                  cgcgGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAAc                                                                                                           <  2:792581/151‑1 (MQ=255)
                                      gTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCa                                                                                                       <  2:707729/151‑1 (MQ=255)
                                             aTCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAAc                                                                                                >  2:1001715/1‑151 (MQ=255)
                                               cATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACg                                                                                               >  2:1570772/1‑150 (MQ=255)
                                                aTTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGAt                                                                                             >  2:395872/1‑151 (MQ=255)
                                                  tGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATa                                                                                            <  1:1654486/150‑1 (MQ=255)
                                                  tGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATa                                                                                            <  1:1521344/150‑1 (MQ=255)
                                                      ctctGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCg                                                                                       >  2:1290810/1‑151 (MQ=255)
                                                      ctctGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCCGAACGATACGGCg                                                                                       >  2:358927/1‑151 (MQ=255)
                                                       tctGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGCTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGt                                                                                      >  2:77001/1‑151 (MQ=255)
                                                        ctGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGtt                                                                                     >  2:829578/1‑151 (MQ=255)
                                                          gACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGttt                                                                                    <  2:69953/150‑1 (MQ=255)
                                                                 ttttATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATgg                                                                             <  1:1196460/150‑1 (MQ=255)
                                                                 ttttATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGgc                                                                            >  2:1425318/1‑151 (MQ=255)
                                                                 ttttATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGgc                                                                            >  1:856264/1‑151 (MQ=255)
                                                                     ataACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGgcgcgc                                                                        >  1:176624/1‑151 (MQ=255)
                                                                     ataACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGgcgcg                                                                         <  1:316638/150‑1 (MQ=255)
                                                                      taACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCt                                                                       >  1:1482300/1‑151 (MQ=255)
                                                                          cGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCttt                                                                     <  2:1284008/149‑1 (MQ=255)
                                                                          cGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTtgt                                                                   <  2:1545957/151‑1 (MQ=255)
                                                                                  gaagaCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACggg                                                           <  1:1212726/151‑1 (MQ=255)
                                                                                     gaCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCa                                                        <  2:898531/151‑1 (MQ=255)
                                                                                      aCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAg                                                       <  2:1528453/151‑1 (MQ=255)
                                                                                          cACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGa                                                    >  2:1534770/1‑150 (MQ=255)
                                                                                           aCTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGaaa                                                  <  1:879442/151‑1 (MQ=255)
                                                                                            cTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGaaa                                                  >  1:850653/1‑150 (MQ=255)
                                                                                            cTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGaaa                                                  >  2:18249/1‑150 (MQ=255)
                                                                                            cTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGaaa                                                  <  2:77169/150‑1 (MQ=255)
                                                                                              tCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGaaaa                                                 >  2:1619089/1‑149 (MQ=255)
                                                                                               cAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGaaaa                                                 >  2:1289544/1‑148 (MQ=255)
                                                                                               cAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAag                                               <  1:53785/150‑1 (MQ=255)
                                                                                                aaGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAagag                                             <  2:1653056/151‑1 (MQ=255)
                                                                                                 aGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGc                                            >  2:1507822/1‑151 (MQ=255)
                                                                                                    atatGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGc                                         <  1:1572088/151‑1 (MQ=255)
                                                                                                     tatGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGc                                         >  2:928195/1‑150 (MQ=255)
                                                                                                       tGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCAcc                                      >  1:1047999/1‑151 (MQ=255)
                                                                                                       tGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCAcc                                      >  2:727595/1‑151 (MQ=255)
                                                                                                          tCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCAt                                    <  1:274359/150‑1 (MQ=255)
                                                                                                          tCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATg                                   >  2:1404285/1‑151 (MQ=255)
                                                                                                             gCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATgcgc                                <  1:313932/151‑1 (MQ=255)
                                                                                                              cATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATgcgc                                >  2:721708/1‑150 (MQ=255)
                                                                                                              cATAAGGTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCa                               <  2:274059/151‑1 (MQ=255)
                                                                                                                 aaGGGTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATg                            >  1:1484210/1‑151 (MQ=255)
                                                                                                                   ggTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGaa                          <  1:1570675/151‑1 (MQ=255)
                                                                                                                    gTTCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAAt                         <  2:1585506/151‑1 (MQ=255)
                                                                                                                     ttCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAAtt                        >  2:319365/1‑151 (MQ=255)
                                                                                                                      tCCCAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAAtt                        <  1:335115/150‑1 (MQ=255)
                                                                                                                        ccAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTc                     <  1:1658464/151‑1 (MQ=255)
                                                                                                                         cAGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCa                    <  2:1485223/151‑1 (MQ=255)
                                                                                                                          aGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTc                     >  1:1369517/1‑149 (MQ=255)
                                                                                                                          aGTCGGTCGCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAg                   <  1:865310/151‑1 (MQ=255)
                                                                                                                            tcggtcgCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAgg                  >  1:1570253/1‑150 (MQ=255)
                                                                                                                              ggtcgCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAAt               <  1:198993/151‑1 (MQ=255)
                                                                                                                              ggtcgCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAAt               <  1:1611077/151‑1 (MQ=255)
                                                                                                                               gtcgCCATATGGAAGACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAAAGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGaa                <  2:1277071/149‑1 (MQ=255)
                                                                                                                                tcgCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGaa                <  1:94983/148‑1 (MQ=255)
                                                                                                                                tcgCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATc              <  2:486366/150‑1 (MQ=255)
                                                                                                                                 cgCCATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCAt            <  1:604260/151‑1 (MQ=255)
                                                                                                                                    cATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCAtttt         <  2:560889/151‑1 (MQ=255)
                                                                                                                                    cATATGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCAtttt         >  2:1690587/1‑151 (MQ=255)
                                                                                                                                     atatGGAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTa        <  1:869299/151‑1 (MQ=255)
                                                                                                                                         ggAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATgc     >  1:1635558/1‑150 (MQ=255)
                                                                                                                                          gAATACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATgc     >  2:45661/1‑149 (MQ=255)
                                                                                                                                           aaTACGC‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATgcagc  <  1:1183797/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑‑CCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGC  >  NC_000913/3102190‑3102478

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑