Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,925,407 T→A 6.4% D20E (GAT→GAA queF → 7‑cyano‑7‑deazaguanine reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,925,4070TA 6.4% 197.4 / 3.5 62D20E (GAT→GAAqueF7‑cyano‑7‑deazaguanine reductase
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base A (2/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCTTGTGCGGTCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCTGGACGCTGTATGAACTTTCCTGGCTGAATGCGAAAG  >  NC_000913/2925258‑2925552
                                                                                                                                                     |                                                                                                                                                 
gCTTGTGCGGTCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGAt                                                                                                                                                   >  1:95157/1‑150 (MQ=255)
     tgCGGTCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCg                                                                                                                                              >  2:47698/1‑150 (MQ=255)
     tgCGGTCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCg                                                                                                                                              >  2:1421744/1‑150 (MQ=255)
     tgCGGTCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCg                                                                                                                                              >  2:1151349/1‑150 (MQ=255)
      gCGGTCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCgg                                                                                                                                             <  1:240332/150‑1 (MQ=255)
         gTCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGAt                                                                                                                                          <  2:945554/150‑1 (MQ=255)
          tCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATa                                                                                                                                         >  1:841990/1‑150 (MQ=255)
              aTCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTa                                                                                                                                     <  1:1421744/150‑1 (MQ=255)
                    cACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGc                                                                                                                               >  2:754806/1‑150 (MQ=255)
                     aCTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTTGTTTATGACGGTGAAACATGTATTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGcc                                                                                                                              <  2:401803/150‑1 (MQ=255)
                       tAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAg                                                                                                                            <  2:432519/150‑1 (MQ=255)
                         aTTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGcc                                                                                                                          <  2:1216697/150‑1 (MQ=255)
                          tttCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGAATACCGGGATACCTATGACGCCAGCCt                                                                                                                         >  2:752677/1‑150 (MQ=255)
                          tttCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGAATACCGGGATACCTATGACGCCAGCCt                                                                                                                         >  2:757730/1‑150 (MQ=255)
                          tttCTCTCTTCACTTAAACCAGATACACTTGTCGTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTa                                                                                                                        >  2:272933/1‑150 (MQ=255)
                              tctctTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTg                                                                                                                     <  1:308725/150‑1 (MQ=255)
                              tctctTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTg                                                                                                                     <  2:912059/150‑1 (MQ=255)
                               ctctTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGc                                                                                                                    <  1:1155681/150‑1 (MQ=255)
                                      cTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGt                                                                                                             <  1:1170246/150‑1 (MQ=255)
                                      cTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGt                                                                                                             >  2:729999/1‑150 (MQ=255)
                                          aaCCAGATACACTTGTCGTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCAc                                                                                                        <  1:272933/150‑1 (MQ=255)
                                           aCCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGAATACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCAc                                                                                                        <  1:757730/150‑1 (MQ=255)
                                           aCCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGAATACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCAc                                                                                                        <  1:752677/150‑1 (MQ=255)
                                            ccAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACg                                                                                                       >  1:1108111/1‑150 (MQ=255)
                                            ccAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGAGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACg                                                                                                       >  1:1330540/1‑150 (MQ=255)
                                              aGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCa                                                                                                     <  1:291146/150‑1 (MQ=255)
                                               gATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAg                                                                                                    >  2:527689/1‑150 (MQ=255)
                                                 tACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGcc                                                                                                  >  1:38324/1‑150 (MQ=255)
                                                   cacTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTg                                                                                                >  1:240409/1‑150 (MQ=255)
                                                          cGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATcgcg                                                                                         >  1:386121/1‑150 (MQ=255)
                                                             ttagtttaTCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGAcc                                                                                      >  2:662732/1‑150 (MQ=255)
                                                              tagtttaTCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCCGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGAccc                                                                                     >  2:23688/1‑150 (MQ=255)
                                                                gtttaTCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGc                                                                                   >  1:325365/1‑150 (MQ=255)
                                                                    aTCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTggg                                                                               >  1:748249/1‑150 (MQ=255)
                                                                    aTCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTggg                                                                               <  1:1360035/150‑1 (MQ=255)
                                                                       tGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCt                                                                            >  2:615083/1‑150 (MQ=255)
                                                                             aTGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCCCGACCCGCTGGGTCTGAAAGc                                                                      <  1:754806/150‑1 (MQ=255)
                                                                              tGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCg                                                                     <  2:1108111/150‑1 (MQ=255)
                                                                              tGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCg                                                                     >  1:361947/1‑150 (MQ=255)
                                                                                aCGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGa                                                                   <  2:748249/150‑1 (MQ=255)
                                                                                 cGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCATCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGAt                                                                  <  2:1052675/150‑1 (MQ=255)
                                                                                       aaCATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTg                                                            <  2:841990/150‑1 (MQ=255)
                                                                                           tGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCtt                                                        <  2:718960/150‑1 (MQ=255)
                                                                                                tcttATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCAc                                                   >  2:635613/1‑150 (MQ=255)
                                                                                                tcttATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCAc                                                   <  1:1151349/150‑1 (MQ=255)
                                                                                                tcttATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCAc                                                   <  1:729999/150‑1 (MQ=255)
                                                                                                    aTGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCCGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTa                                               >  1:928237/1‑150 (MQ=255)
                                                                                                      gCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACg                                             >  1:645444/1‑150 (MQ=255)
                                                                                                      gCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACg                                             >  2:393835/1‑150 (MQ=255)
                                                                                                        aaaCCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGa                                           >  2:1123699/1‑150 (MQ=255)
                                                                                                                   cACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCTGGACGc                                >  1:914012/1‑150 (MQ=255)
                                                                                                                    agttGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCTGGACGCt                               <  2:438501/148‑1 (MQ=255)
                                                                                                                          gggCCTGACTCTTGGAAAATCAAGCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCTGGACGCTGTATGa                         <  2:65364/150‑1 (MQ=255)
                                                                                                                              cTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCTGGACGCTGTATGAACtt                     <  1:662732/150‑1 (MQ=255)
                                                                                                                                    ctTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCTGGACGCTGTATGAACTTTCCTgg               >  1:396256/1‑150 (MQ=255)
                                                                                                                                      tGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGAAATCTGGACGCTGTATGAACTTTCCAGGCt             >  1:563497/1‑150 (MQ=255)
                                                                                                                                       ggAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCTGGACGCTGTATGAACTTTCCTGGCTg            <  2:914012/150‑1 (MQ=255)
                                                                                                                                            aTCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCTGGACGCTGTATGAACTTTCCTGGCTGAATGc       >  1:348143/1‑150 (MQ=255)
                                                                                                                                             tCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCTGGACGCTGTATGAACTTTCCTGGCTGAATGCg      <  2:38324/150‑1 (MQ=255)
                                                                                                                                              cAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCTGGACGCTGTATGAACTTTCCTGGCTGAATGc       >  2:301372/1‑148 (MQ=255)
                                                                                                                                              cAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCTGGACGCTGTATGAACTTTCCTGGCTGAATGc       <  1:301372/148‑1 (MQ=255)
                                                                                                                                              cAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCTGGACGCTGTATGAACTTTCCTGGCTGAATGCGa     >  1:372039/1‑150 (MQ=255)
                                                                                                                                                 ccGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCTGGACGCTGTATGAACTTTCCTGGCTGAATGCGAAAg  <  1:47698/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GCTTGTGCGGTCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCTGGACGCTGTATGAACTTTCCTGGCTGAATGCGAAAG  >  NC_000913/2925258‑2925552

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