Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,573,873 C→T 6.6% A444T (GCG→ACG)  pqqL ← putative periplasmic M16 family zinc metalloendopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,573,8730CT 6.6% 186.6 / 2.7 60A444T (GCG→ACG) pqqLputative periplasmic M16 family zinc metalloendopeptidase
Reads supporting (aligned to +/- strand):  ref base C (30/26);  new base T (2/2);  total (32/28)
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.90e-01

ACGTAAGCCGCCAGCTTTTTGTTGGCGTACTCCTTTTCCAGCGCCAGGATAGCCGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGC  >  NC_000913/1573725‑1574020
                                                                                                                                                    |                                                                                                                                                   
aCGTAAGCCGCCAGCTTTTTGTTGGCGTACTCCTTTTCCAGCGCCAGGATAGCCGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGc                                                                                                                                                     <  3:219772/149‑1 (MQ=255)
    aaGCCGCCAGCTTTTTGTTGGCGTACTCCTTTTCCAGCGCCAGGATAGCCGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGa                                                                                                                                                 >  4:32798/1‑149 (MQ=255)
        cgccAGCTTTTTGTTGGCGTACTCCTTTTCCAGCGCCAGGATAGCCGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACttt                                                                                                                                             >  4:10612/1‑149 (MQ=255)
        cgccAGCTTTTTGTTGGCGTACTCCTTTTCCAGCGCCAGGATAGCCGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACttt                                                                                                                                             <  4:300400/149‑1 (MQ=255)
                     ttgGCGTACTCCTTTTCCAGCGCCAGGATAGCCGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTtg                                                                                                                                >  1:126168/1‑149 (MQ=255)
                       ggCGTACTCCTTTTCCAGCGCCAGGATAGCCGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTtgct                                                                                                                              <  4:282012/149‑1 (MQ=255)
                              tacTTTTCCAGCGCCATGATAGCCGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTATTTACCATTTGCACAAAAAATGTATCCTGGTTCTTTATTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAg                                                                                                                       >  4:393914/3‑149 (MQ=255)
                              tCCTTTTCCAGCGCCAGGATAGCCGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTTCCACAg                                                                                                                       >  1:346767/1‑149 (MQ=255)
                              tCCTTTTCCAGCGCCAGGATAGCCGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAg                                                                                                                       >  1:371403/1‑149 (MQ=255)
                               ccTTTTCCAGCGCCAGGATAGCCGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGTCAGACTTTGCACGGTAATTTGCTGCCACAGa                                                                                                                      >  3:388295/1‑149 (MQ=255)
                               ccTTTTCCAGCGCCAGGATAGCCGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGTCTTTGCACGGTAATTTGCTGCCACAGa                                                                                                                      <  3:518800/149‑1 (MQ=255)
                                     ccAGCGCCAGGATAGCCGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTc                                                                                                                <  2:346767/149‑1 (MQ=255)
                                         cgcCAGGATAGCCGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGaaa                                                                                                            <  1:140468/149‑1 (MQ=255)
                                         cgcCAGGATAGCCGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGaaa                                                                                                            <  2:16553/149‑1 (MQ=255)
                                               gATAGCCGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCAATTTTTCCGCAAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGat                                                                                                      <  3:369520/149‑1 (MQ=255)
                                                    ccGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGTCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGtc                                                                                                 >  3:522459/1‑149 (MQ=255)
                                                             aCAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGc                                                                                        <  3:32798/149‑1 (MQ=255)
                                                               aaTGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGtt                                                                                                         >  4:117623/1‑130 (MQ=255)
                                                               aaTGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGtt                                                                                                         <  3:117623/130‑1 (MQ=255)
                                                                  gCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACaa                                                                                   >  1:272618/1‑149 (MQ=255)
                                                                     tttttGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGTCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGaa                                                                                <  4:522459/149‑1 (MQ=255)
                                                                            cggcAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGta                                                                         >  1:191003/1‑149 (MQ=255)
                                                                              gcAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGtatt                                                                       <  3:10612/149‑1 (MQ=255)
                                                                                aaCCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGtattat                                                                     >  3:109918/1‑149 (MQ=255)
                                                                                aaCCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGtattat                                                                     <  4:352571/149‑1 (MQ=255)
                                                                                 aCCTCGTTGTTTACCATTTGCTCCCAACATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTattatt                                                                    <  4:506259/149‑1 (MQ=255)
                                                                                    tCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTaa                                                                 <  4:586122/149‑1 (MQ=255)
                                                                                    tCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTaa                                                                 >  2:282194/1‑149 (MQ=255)
                                                                                     cATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAAt                                                                <  2:336102/149‑1 (MQ=255)
                                                                                              ccATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGAt                                                       >  4:518123/1‑149 (MQ=255)
                                                                                                      tCCCAAAATGCGTCCTGGTGCTTTCTTAACTGCTGCCATTTTGCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCg                                               >  3:344975/1‑149 (MQ=255)
                                                                                                        ccAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGAc                                             >  4:501316/1‑149 (MQ=255)
                                                                                                        ccAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGAc                                             <  3:501316/149‑1 (MQ=255)
                                                                                                        ccAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGAc                                             <  3:518123/149‑1 (MQ=255)
                                                                                                        ccAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTGCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGAc                                             >  1:173788/1‑149 (MQ=255)
                                                                                                         cAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTTGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACt                                            >  1:110959/1‑149 (MQ=255)
                                                                                                         cAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGta                                                                         <  3:269008/120‑1 (MQ=255)
                                                                                                         cAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGta                                                                         >  4:269008/1‑120 (MQ=255)
                                                                                                         cAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACt                                            >  1:172891/1‑149 (MQ=255)
                                                                                                          aaaaTGCGTCCTGGTTCTTTCTTACCTGCTGCCATTTTTCCGCCAGACTTTGCACGGTGATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCGCTTCCGGCGACAAGAACTGCGTATTATTTAATGAGCTGGATGCCAGGCGACTg                                           <  3:66180/149‑1 (MQ=255)
                                                                                                             aTGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTTTTATTTAATGa                                                              <  3:461210/127‑1 (MQ=255)
                                                                                                             aTGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGa                                                              >  4:461210/1‑127 (MQ=255)
                                                                                                                  tCCTGGTTCTTTCTTCACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCAt                                   <  2:173788/149‑1 (MQ=255)
                                                                                                                   ccTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATa                                  <  3:261006/149‑1 (MQ=255)
                                                                                                                    cTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGt                                         >  3:579758/1‑141 (MQ=255)
                                                                                                                    cTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGt                                         <  4:579758/141‑1 (MQ=255)
                                                                                                                    cTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATAc                                 >  3:355292/1‑149 (MQ=255)
                                                                                                                     tGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACg                                >  1:195346/1‑149 (MQ=255)
                                                                                                                        ttctttctTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTaa                             >  2:545840/1‑149 (MQ=255)
                                                                                                                         tctttcttAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTaaa                            <  1:41824/149‑1 (MQ=255)
                                                                                                                           tttcttAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATc                          >  2:505577/1‑149 (MQ=255)
                                                                                                                                tAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGtt                     >  3:567527/1‑149 (MQ=255)
                                                                                                                                    tgctgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGc                 <  1:505577/149‑1 (MQ=255)
                                                                                                                                        gcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGc             <  2:371403/149‑1 (MQ=255)
                                                                                                                                          cATTTTTCCGTCAGACTTTGCACGGTAATTTGCTGCCACAGACGTCTCGAAAGTTGATATGTCTCTTCCGGCGa                                                                                      <  2:426522/74‑1 (MQ=255)
                                                                                                                                          cATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGa                                                                                      >  1:426522/1‑74 (MQ=255)
                                                                                                                                          cATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTg           >  1:261141/1‑149 (MQ=255)
                                                                                                                                           aTTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGa          >  1:202232/1‑149 (MQ=255)
                                                                                                                                              tttCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCAAGCTGaacg       <  4:352479/149‑4 (MQ=255)
                                                                                                                                                   gCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGc  >  2:292862/1‑149 (MQ=255)
                                                                                                                                                   gCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGc  >  1:92456/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACGTAAGCCGCCAGCTTTTTGTTGGCGTACTCCTTTTCCAGCGCCAGGATAGCCGCAGGAGACAATGCTTTTTTGGCGGCAACCTCATTGTTTACCATTTGCTCCCAAAATGCGTCCTGGTTCTTTCTTAACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGC  >  NC_000913/1573725‑1574020

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑