Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 707,879 A→T L11I (TTA→ATA)  mntB ← ABC superfamily ATP binding cassette transporter, ABC protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730707,8790AT98.6% 241.7 / ‑3.7 69L11I (TTA→ATA) mntBABC superfamily ATP binding cassette transporter, ABC protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (41/27);  total (42/27)
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.97e-01

CCGGTAGCATTAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAAT  >  CP000730/707743‑708013
                                                                                                                                        |                                                                                                                                      
ccGGTAGCATTAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAAc                                                                                                                                    >  1:672782/1‑141 (MQ=255)
     aGCATTAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGAtt                                                                                                                               <  2:465684/141‑1 (MQ=255)
      gCATTAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGAttt                                                                                                                              >  2:25395/1‑141 (MQ=255)
      gCATTAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGAttt                                                                                                                              <  2:83363/141‑1 (MQ=255)
       cATTAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTa                                                                                                                             >  2:331127/1‑141 (MQ=255)
        aTTAAATTCACCAATTAAAGACTTGATAACGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTaa                                                                                                                            <  1:352911/141‑1 (MQ=255)
         ttAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTaaa                                                                                                                           <  2:322300/141‑1 (MQ=255)
         ttAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGACTTaaa                                                                                                                           >  2:620436/1‑141 (MQ=255)
                cACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTg                                                                                                                    <  2:397643/141‑1 (MQ=255)
                cACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTg                                                                                                                    >  1:27952/1‑141 (MQ=255)
                 aCCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGt                                                                                                                   >  2:591678/1‑141 (MQ=255)
                 aCCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGt                                                                                                                   >  2:268338/1‑141 (MQ=255)
                  ccAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGtt                                                                                                                  >  2:802701/1‑141 (MQ=255)
                     aTTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCt                                                                                                               >  2:245695/1‑141 (MQ=255)
                     aTTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCt                                                                                                               <  2:398841/141‑1 (MQ=255)
                      ttAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTa                                                                                                              <  2:637835/141‑1 (MQ=255)
                          aGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACaa                                                                                                          >  1:740103/1‑141 (MQ=255)
                          aGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACaa                                                                                                          >  1:721601/1‑141 (MQ=255)
                          aGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACaa                                                                                                          >  1:721487/1‑141 (MQ=255)
                           gACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACaaa                                                                                                         >  2:50556/1‑141 (MQ=255)
                            aCTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGCAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTACATCTTTTGTTTCTAACAAAc                                                                                                        >  1:79598/1‑141 (MQ=255)
                             cTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACg                                                                                                       >  1:669403/1‑141 (MQ=255)
                              ttGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGt                                                                                                      <  2:833512/141‑1 (MQ=255)
                                        aaGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATAcctcc                                                                                            <  2:540140/141‑1 (MQ=255)
                                            tttACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAAt                                                                                        >  2:196743/1‑141 (MQ=255)
                                            tttACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAAt                                                                                        >  1:795097/1‑141 (MQ=255)
                                            tttACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAAt                                                                                        >  1:290423/1‑141 (MQ=255)
                                            tttACCAGAACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAAt                                                                                        <  2:36111/141‑1 (MQ=255)
                                              tACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTa                                                                                      >  1:598413/1‑141 (MQ=255)
                                              tACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTa                                                                                      >  2:87320/1‑141 (MQ=255)
                                               aCCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTaa                                                                                     >  1:166360/1‑141 (MQ=255)
                                               aCCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTaa                                                                                     <  1:26708/141‑1 (MQ=255)
                                                ccAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTaaa                                                                                    >  1:322021/1‑141 (MQ=255)
                                                ccAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTaaa                                                                                    >  2:590937/1‑141 (MQ=255)
                                                  aGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAg                                                                                  >  1:412625/1‑141 (MQ=255)
                                                    cACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGtt                                                                                <  1:692315/141‑1 (MQ=255)
                                                          tCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCt                                                                          >  2:384360/1‑141 (MQ=255)
                                                          tCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCt                                                                          >  1:522903/1‑141 (MQ=255)
                                                            gggCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTaa                                                                        <  2:820849/141‑1 (MQ=255)
                                                             ggCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAAc                                                                       <  2:29239/141‑1 (MQ=255)
                                                             ggCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAAc                                                                       <  1:452721/141‑1 (MQ=255)
                                                               cccATGATCCCAATTATTTCGACGCGTACAGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCt                                                                     <  2:639671/141‑1 (MQ=255)
                                                               cccATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCt                                                                     <  2:125093/141‑1 (MQ=255)
                                                                           attattTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGtat                                                         >  2:448573/1‑141 (MQ=255)
                                                                            ttattTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTTAAAGTTTAGGCTAACCTAATTAATTGtata                                                        <  1:177701/141‑1 (MQ=255)
                                                                                  cgccgcGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAAc                                                  >  2:595601/1‑141 (MQ=255)
                                                                                   gccgcGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACt                                                 >  2:16119/1‑141 (MQ=255)
                                                                                    ccgcGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTg                                                <  2:68713/141‑1 (MQ=255)
                                                                                      gcgTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTgag                                              <  1:196743/141‑1 (MQ=255)
                                                                                        gTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACCAACGATTACACCACCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGaa                                            <  1:735275/141‑1 (MQ=255)
                                                                                                      aaGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCa                              <  1:70921/141‑1 (MQ=255)
                                                                                                       aGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCaa                             >  2:555067/1‑141 (MQ=255)
                                                                                                       aGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCaa                             <  1:754943/141‑1 (MQ=255)
                                                                                                          aaaTGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTa                          >  2:749758/1‑141 (MQ=255)
                                                                                                                  ttAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTg                  >  1:1900/1‑141 (MQ=255)
                                                                                                                   tAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGa                 >  2:587170/1‑141 (MQ=255)
                                                                                                                   tAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGa                 >  2:587137/1‑141 (MQ=255)
                                                                                                                    aaGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGat                >  1:382150/1‑141 (MQ=255)
                                                                                                                      gTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGatat              >  2:149050/1‑141 (MQ=255)
                                                                                                                      gTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGatat              <  2:705240/141‑1 (MQ=255)
                                                                                                                       tACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGatata             >  1:296169/1‑141 (MQ=255)
                                                                                                                        aCATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGatataa            >  1:472928/1‑141 (MQ=255)
                                                                                                                         cATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTAGACCTCCGAAGTAATAGTTTAGGCGAACCTAATTATGTGGAGAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGatataat           >  1:82018/1‑141 (MQ=255)
                                                                                                                           tGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAAAATATTTATCATGTCAAGTAAATTCGTGatataatat         <  1:749758/141‑1 (MQ=255)
                                                                                                                               tattaCCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGAc     <  2:245319/141‑1 (MQ=255)
                                                                                                                               tattaCCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGAc     <  1:595601/141‑1 (MQ=255)
                                                                                                                               tattaCCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGAc     >  1:818425/1‑141 (MQ=255)
                                                                                                                                 ttaCCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACaa   >  2:353313/1‑141 (MQ=255)
                                                                                                                                  taCCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAat  >  1:755589/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                      
CCGGTAGCATTAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAAT  >  CP000730/707743‑708013

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