Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,018,810 (T)7→6 coding (109/810 nt) mgt ← transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,018,8040T.96.9% 284.6 / 1.3 64coding (115/810 nt)mgttransglycosylase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base . (25/37);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

AGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAA  >  CP000730/2018664‑2018939
                                                                                                                                            |                                                                                                                                       
aGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACt                                                                                                                                         <  2:132603/141‑1 (MQ=255)
   tGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACtttt                                                                                                                                      <  1:780250/141‑1 (MQ=255)
     acacAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACtttttt                                                                                                                                    <  2:842827/141‑1 (MQ=255)
       acaAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACttttttct                                                                                                                                  >  1:155519/1‑139 (MQ=255)
       acaAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACttttttct                                                                                                                                  <  1:177651/141‑3 (MQ=255)
              tactTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTcgg                                                                                                                          <  1:888627/141‑1 (MQ=255)
              tactTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTcgg                                                                                                                          >  1:508007/1‑141 (MQ=255)
                 ttttATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTtcgtcgt                                                                                                                       >  2:129194/1‑140 (MQ=255)
                  tttATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCTGTGGt                                                                                                                      >  2:863591/1‑141 (MQ=255)
                  tttATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGt                                                                                                                      >  1:365639/1‑141 (MQ=255)
                       tttCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGcc                                                                                                                 >  1:817565/1‑141 (MQ=255)
                       tttCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGcc                                                                                                                 <  1:721623/141‑1 (MQ=255)
                         tCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCaa                                                                                                               <  1:469841/141‑1 (MQ=255)
                         tCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCaa                                                                                                               <  2:261615/141‑1 (MQ=255)
                          cAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAAc                                                                                                              >  1:855274/1‑141 (MQ=255)
                          cAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAACGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAAc                                                                                                              >  1:292712/1‑141 (MQ=255)
                            aTTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTg                                                                                                            <  1:88675/141‑1 (MQ=255)
                             tttttCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTgg                                                                                                           <  2:316761/141‑1 (MQ=255)
                                ttCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTg                                                                                                        >  1:490566/1‑141 (MQ=255)
                                 tCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGt                                                                                                          <  2:385039/138‑1 (MQ=255)
                                 tCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGt                                                                                                          >  1:385039/1‑138 (MQ=255)
                                        tCATCCACATAATCGCGAATAGATAAAAAATACATAATACCAATAAGCAATGCTATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACg                                                                                                <  2:648449/141‑1 (MQ=255)
                                          aTCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGta                                                                                              <  1:259249/141‑1 (MQ=255)
                                           tCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGtat                                                                                             >  1:616095/1‑141 (MQ=255)
                                               caTTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATtata                                                                                         <  1:532025/141‑1 (MQ=255)
                                                 ttATCGCGTGTATATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGt                                                                                       <  2:673476/141‑1 (MQ=255)
                                                      gcgTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGt                                                                                  <  1:554817/141‑1 (MQ=255)
                                                        gtgtAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTc                                                                                <  1:497448/141‑1 (MQ=255)
                                                        gtgtAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTc                                                                                <  1:812147/141‑1 (MQ=255)
                                                           tAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCgtg                                                                             <  2:866279/141‑1 (MQ=255)
                                                            aGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGTATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGTTTGATACTACGTATTATAGTGAGGTTCgtga                                                                            >  2:102357/1‑140 (MQ=255)
                                                            aGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCgtgt                                                                            <  1:893630/141‑1 (MQ=255)
                                                              aTAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGttt                                                                          <  1:72624/141‑1 (MQ=255)
                                                                   aaaTACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATAGTTTTAATCGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCctat                                                                     >  1:536190/1‑141 (MQ=255)
                                                                     aTACATGATACCACTAAACTATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGc                                                                   <  1:741136/141‑1 (MQ=255)
                                                                     aTACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGTTAATACGTATTATAGTGAGGTTCGTGTTTCATATGc                                                                   <  2:276738/141‑1 (MQ=255)
                                                                        cATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCa                                                                >  1:768946/1‑141 (MQ=255)
                                                                          tGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCaaa                                                              <  2:549000/141‑1 (MQ=255)
                                                                             tACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATg                                                           <  2:566860/141‑1 (MQ=255)
                                                                                  aTAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCAt                                                      >  2:300600/1‑141 (MQ=255)
                                                                                       cAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCAtttgag                                                 <  1:520712/141‑1 (MQ=255)
                                                                                           gCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATtata                                                                                         >  2:483027/1‑97 (MQ=255)
                                                                                           gCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATtata                                                                                         <  1:483027/97‑1 (MQ=255)
                                                                                                 aTGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCTGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTAcc                                       >  2:738887/1‑141 (MQ=255)
                                                                                                 aTGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTAcc                                       >  1:224678/1‑141 (MQ=255)
                                                                                                      tAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTc                                                                                <  2:371215/95‑1 (MQ=255)
                                                                                                      tAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTc                                                                                >  1:371215/1‑95 (MQ=255)
                                                                                                       aGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGc                                 <  1:608818/141‑1 (MQ=255)
                                                                                                        gAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGct                                <  1:641993/141‑1 (MQ=255)
                                                                                                        gAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGct                                <  1:908426/141‑1 (MQ=255)
                                                                                                          aTGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGcttc                              <  1:864061/141‑1 (MQ=255)
                                                                                                               taataaTATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTc                         <  1:126537/141‑1 (MQ=255)
                                                                                                                     tatTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGt                   <  1:387237/141‑1 (MQ=255)
                                                                                                                     tatTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGt                   <  2:422714/141‑1 (MQ=255)
                                                                                                                     tatTTTTAATAGTATTCGTTTAC‑GTTTTTCTTTTTCGGAGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTGATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGt                   <  2:422730/141‑1 (MQ=255)
                                                                                                                      atTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTcggcgg                                                                                                                       >  1:763720/1‑40 (MQ=255)
                                                                                                                      atTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTcggcgg                                                                                                                       <  2:763720/40‑1 (MQ=255)
                                                                                                                        ttttAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGCGGTTCGAGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTg                >  1:270342/1‑141 (MQ=255)
                                                                                                                         tttAATAGTATTCGTTTACTTTTTTT‑TTTTTTTTTTTTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTTAGTTTGTGTACTTATCTTTTCTTTTCATTCGttttt               >  2:611844/1‑139 (MQ=255)
                                                                                                                         tttAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGc               <  2:621747/141‑1 (MQ=255)
                                                                                                                         tttAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGc               >  1:801537/1‑141 (MQ=255)
                                                                                                                          ttAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCt              >  1:661535/1‑141 (MQ=255)
                                                                                                                          ttAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCt              <  2:176751/141‑1 (MQ=255)
                                                                                                                          ttAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGCTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCt              >  1:803173/1‑141 (MQ=255)
                                                                                                                           tAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCttt                           >  1:236892/1‑127 (MQ=255)
                                                                                                                           tAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCttt                           <  2:236892/127‑1 (MQ=255)
                                                                                                                            aaTAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTcc            >  1:916566/1‑141 (MQ=255)
                                                                                                                              tAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCctt          <  2:485744/141‑1 (MQ=255)
                                                                                                                                      gTTTACTTTTTTT‑TTTTTTTTCGGATTGCCAACTGGTTGATAATACGTATTATTGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCTTTTGAGTTTGAGTACATATTGCTTTTTTTTATTCGTTTGCTACTTTTTTTtaa  >  2:270143/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                       
AGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAA  >  CP000730/2018664‑2018939

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