Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,438,346 C→T 6.2% Q204* (CAG→TAG)  lysA2 → diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,438,3460CT 6.2% 197.1 / 3.8 65Q204* (CAG→TAG) lysA2diaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (34/26);  new base T (2/2);  total (37/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 = 6.26e-01

CACACGAATTTATTCAAACTGGGCAAGAAGATAGTAAGTTTGGATTATCAATTCAATATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTAAATATGTTGAAGGTGACGAA  >  CP000730/1438206‑1438477
                                                                                                                                            |                                                                                                                                   
cacacGAATTTATTCAAACTGGGCAAGAAGATAGTAAGTTTGGATTATCAATTCAATATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTcac                                                                                                                                     >  1:781271/1‑141 (MQ=255)
      aaTTTATTCAAACTGGGCAAGAAGATAGTAAGTTTGGATTATCAATTCAATATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTg                                                                                                                               >  1:353249/1‑141 (MQ=255)
          tATTCAAACTGGGCAAGAAGATAGTAAGTTTGGATTATCAATTCAATATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAgg                                                                                                                           >  1:312316/1‑141 (MQ=255)
                     ggCAAGAAGATAGTAAGTTTGGATTATCAATTCAATATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCAt                                                                                                                <  2:329073/141‑1 (MQ=255)
                     ggCAAGAAGATAGTAAGTGTGGATTATCAATTCAATATGGCATAGCAAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCAt                                                                                                                <  2:49935/141‑1 (MQ=255)
                             gaTAGTAAGTTTGGATTATCAATTCAATATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGaa                                                                                                        >  2:640747/1‑141 (MQ=255)
                               tAGTAAGTTTGGATTATCAATTCAATATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAAc                                                                                                      >  1:415619/1‑141 (MQ=255)
                               tAGTAAGTTTGGATTATCAAGTCAATATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTGAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAAc                                                                                                      >  2:274634/1‑141 (MQ=255)
                                  tAAGTGTGGATTATCAATTCAATATTGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGc                                                                                                   <  2:777116/141‑1 (MQ=255)
                                       ttGGATTATCAATTCAATATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTaaaa                                                                                              <  2:312316/141‑1 (MQ=255)
                                          gATTATCAATTCAATATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTg                                                                                           >  2:776757/1‑141 (MQ=255)
                                             tATCAATTCAATATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGttt                                                                                        >  1:213436/1‑141 (MQ=255)
                                             tATCAATTCAATATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGttt                                                                                        >  1:480232/1‑141 (MQ=255)
                                              aTCAATTCAATATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGtttt                                                                                       <  1:85014/141‑1 (MQ=255)
                                              aTCAATTCAATATGGATTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAAGTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGCAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGtttt                                                                                       <  1:576350/141‑1 (MQ=255)
                                                 aattcaatATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACg                                                                                    >  2:745728/1‑141 (MQ=255)
                                                  attcaatATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGt                                                                                   >  1:722933/1‑141 (MQ=255)
                                                  attcaatATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGt                                                                                   >  2:633264/1‑141 (MQ=255)
                                                  attcaatATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGt                                                                                   >  2:652236/1‑141 (MQ=255)
                                                  attcaatATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATAAAAAGGTGTACATTGTCATATTGGATAACAGATTGAAGGTACAAAAGCATTTATAGAAACTGATAAAATTGTTGTACGt                                                                                   >  2:21728/1‑141 (MQ=255)
                                                    tcaatATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTg                                                                                 <  2:197021/141‑1 (MQ=255)
                                                    tcaatATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTg                                                                                 <  2:680714/141‑1 (MQ=255)
                                                            gCTTAGCTAAAAAAGCAATTGACCAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAACTTGTTTTACGTTGGGTTAAag                                                                         >  2:273156/1‑141 (MQ=255)
                                                            gCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAag                                                                         >  1:608408/1‑141 (MQ=255)
                                                              ttAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAagag                                                                       <  1:661343/141‑1 (MQ=255)
                                                                aGCTAAAAAAGCTATTGACAACGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGAGTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCa                                                                     >  1:419612/1‑141 (MQ=255)
                                                                aGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGCGCa                                                                     >  2:81506/1‑141 (MQ=255)
                                                                aGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTAAGTTGGCTTAAAGAGCa                                                                     >  2:361894/1‑141 (MQ=255)
                                                                  cTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAg                                                                   <  1:176263/141‑1 (MQ=255)
                                                                     aaaaaGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTGGGCTTAAAGAGCAAGGCa                                                                >  1:238544/1‑141 (MQ=255)
                                                                        aaGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTc                                                             <  2:781271/141‑1 (MQ=255)
                                                                            aaTTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGt                                                         <  2:799269/141‑1 (MQ=255)
                                                                                 aCAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAGGTTGAat                                                    >  1:670363/1‑141 (MQ=255)
                                                                                          aacaaTCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCATAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAAcc                                           >  1:275770/1‑141 (MQ=255)
                                                                                          aacaaTCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCATAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAAcc                                           >  1:275992/1‑141 (MQ=255)
                                                                                          aacaaTCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAAcc                                           >  1:553349/1‑141 (MQ=255)
                                                                                                cTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTtggtg                                     >  2:391473/1‑141 (MQ=255)
                                                                                                cTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTGGAATTATTAAACCTgggcg                                     >  1:654899/1‑139 (MQ=255)
                                                                                                 tAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTtggtgg                                    >  2:404434/1‑141 (MQ=255)
                                                                                                 tAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTtggtgg                                    <  1:271936/141‑1 (MQ=255)
                                                                                                 tAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTtggtgg                                    <  1:289232/141‑1 (MQ=255)
                                                                                                 caaaCACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTtggtgg                                    <  1:295876/140‑1 (MQ=255)
                                                                                                       cTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACATTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCtt                              <  1:539919/141‑1 (MQ=255)
                                                                                                           aaattaaaaGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGt                          >  2:242182/1‑141 (MQ=255)
                                                                                                           aaattaaaaGGTGTACATTGTCATATTGGATCAAAGAATGAAGGTACAGAAGCATGTATTGAAACGGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATACAAGTTGAATGATTTAACCTTGGTGGTGGCTTTGGt                          >  2:385474/1‑141 (MQ=255)
                                                                                                              ttaaaaGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGTCTTTGGTAtt                       <  1:585894/141‑1 (MQ=255)
                                                                                                                aaaaGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTaa                     >  1:233400/1‑141 (MQ=255)
                                                                                                                aaaaGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGCCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGCTGGTGGCTTTGGTATTaa                     >  1:43428/1‑141 (MQ=255)
                                                                                                                 aaaGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTaaa                    >  1:123344/1‑141 (MQ=255)
                                                                                                                 aaaGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGGATTaaa                    >  2:805132/1‑141 (MQ=255)
                                                                                                                      tttACATTGTCTTATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATCATTAAACCTTGGTGGTGGCTTTGGTATTAAATATGt               <  1:408717/139‑1 (MQ=255)
                                                                                                                      tgtACATTTTAAGATTGGCTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTAAATATGt               <  1:247988/141‑1 (MQ=255)
                                                                                                                      tgtACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTAAATATGt               <  1:443823/141‑1 (MQ=255)
                                                                                                                      tgtACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTAAATATGt               <  1:372568/141‑1 (MQ=255)
                                                                                                                          cATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAATGTGCTAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTAAATATGTTGaa           >  1:187994/1‑141 (MQ=255)
                                                                                                                             tGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTAAATATGTTGAAGGt        <  1:274634/141‑1 (MQ=255)
                                                                                                                             tGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTAAATATGTTGAAGGt        <  2:415619/141‑1 (MQ=255)
                                                                                                                               tCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTAAATATGTTGAAGGTGa      >  1:313013/1‑141 (MQ=255)
                                                                                                                               tCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTAAATATGTTGAAGGTGa      <  1:408584/141‑1 (MQ=255)
                                                                                                                                 atatTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGAGGCTTTGGTATTAAATATGTTGAAGGTGACg    >  2:566863/1‑141 (MQ=255)
                                                                                                                                  tatTGGTTCATAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTAAATATGTTGAAGGTGACGa   <  2:275992/141‑1 (MQ=255)
                                                                                                                                  tatTGGTTCATAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTAAATATGTTGAAGGTGACGa   <  2:275770/141‑1 (MQ=255)
                                                                                                                                  tatTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTAAATATGTTGAAGGTGACGa   <  1:655856/141‑1 (MQ=255)
                                                                                                                                  tatTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTAAATATGTTGAAGGTGACGa   <  2:553349/141‑1 (MQ=255)
                                                                                                                                   atTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTAAATATGTTGAAGGTGACGaa  >  2:86452/1‑141 (MQ=255)
                                                                                                                                   atTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTAAAAATGTTGAAGGTGAGGaa  >  2:290516/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                   
CACACGAATTTATTCAAACTGGGCAAGAAGATAGTAAGTTTGGATTATCAATTCAATATGGCTTAGCTAAAAAAGCAATTGACAAAGTCCAACAATCTAAACACTTAAAATTAAAAGGTGTACATTGTCATATTGGTTCACAGATTGAAGGTACAGAAGCATTTATTGAAACTGCTAAAATTGTTTTACGTTGGCTTAAAGAGCAAGGCATTCAAGTTGAATTATTAAACCTTGGTGGTGGCTTTGGTATTAAATATGTTGAAGGTGACGAA  >  CP000730/1438206‑1438477

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: