Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,620,134 C→T 8.1% A163A (GCG→GCA rnz ← ribonuclease Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,620,1340CT 8.1% 199.8 / 4.1 69A163A (GCG→GCArnzribonuclease Z
Reads supporting (aligned to +/- strand):  ref base C (22/41);  new base T (3/3);  total (25/44)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

AAATTGCTACAACTGGACCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGATCATCGATTTCAATAT  >  CP000730/1619998‑1620270
                                                                                                                                        |                                                                                                                                        
aaaTTGCTACAACTGGACCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTATTGATATTTAGGACCTGGTTCATGACCAATGTTTTTTAACGCtt                                                                                                                                      >  2:627340/1‑141 (MQ=255)
      cTACAACTGGACCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAATGCTTCAACat                                                                                                                                >  2:725935/1‑141 (MQ=255)
      cTACAACTGGACCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAATGCTTCAACat                                                                                                                                >  2:725781/1‑141 (MQ=255)
      cTACAACTGGACCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACat                                                                                                                                <  1:301623/141‑1 (MQ=255)
         cAACTGGCCCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACattga                                                                                                                             <  1:756707/141‑1 (MQ=255)
                aCCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTAcc                                                                                                                      <  2:110038/141‑1 (MQ=255)
                  ccTGTTTGGACTCATCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGGATTGTAACGCTTCAACATTGATTGTACCTg                                                                                                                    <  2:602886/141‑1 (MQ=255)
                       ttGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATAAATCCAGAGAGTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGtt                                                                                                               <  1:210629/141‑1 (MQ=255)
                         ggACTGAACTCTAAATTCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGACATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTCCCTGTAGTTTc                                                                                                             <  1:797500/141‑1 (MQ=255)
                           aCTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAg                                                                                                           >  2:574908/1‑141 (MQ=255)
                           aCTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAg                                                                                                           >  2:633564/1‑141 (MQ=255)
                             tCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGt                                                                                                         <  1:764700/141‑1 (MQ=255)
                              cACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTg                                                                                                        <  2:334080/141‑1 (MQ=255)
                              cACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTg                                                                                                        <  1:730154/141‑1 (MQ=255)
                              cACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACACTGATTCTACCTGTAGTTTCAGGTg                                                                                                        <  2:670487/141‑1 (MQ=255)
                              cACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATCAAAATTAAGTTCTTGATATTTAGGACCTCGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTg                                                                                                        <  1:898050/141‑1 (MQ=255)
                               aCCTCTAAAATCTTTTGATTGATAAACTTGGCCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGc                                                                                                       <  2:131472/141‑1 (MQ=255)
                               aCCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGc                                                                                                       <  1:453154/141‑1 (MQ=255)
                                    tAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTGGTTTCAGGTGCCTTTa                                                                                                  >  2:103043/1‑141 (MQ=255)
                                    tAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTa                                                                                                  >  1:719336/1‑141 (MQ=255)
                                    tAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTa                                                                                                  >  2:487692/1‑141 (MQ=255)
                                           ttttGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAATGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCtata                                                                                           <  1:725935/141‑1 (MQ=255)
                                           ttttGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAATGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCtata                                                                                           <  1:725781/141‑1 (MQ=255)
                                            tttGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTatat                                                                                          <  1:82851/141‑1 (MQ=255)
                                            tttGATTGATAAACTTGACCAATATGTTCAAATGCATCACGAGATGTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTatat                                                                                          <  2:399806/141‑1 (MQ=255)
                                              tgattgatAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATcc                                                                                        >  2:800362/1‑141 (MQ=255)
                                               gattgatAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCa                                                                                       >  1:297462/1‑141 (MQ=255)
                                               gattgatAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCGTGAAATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCa                                                                                       >  1:604803/1‑141 (MQ=255)
                                                attgatAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCat                                                                                      <  1:879139/141‑1 (MQ=255)
                                                    atAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATatga                                                                                  >  2:275131/1‑141 (MQ=255)
                                                    aaaaaCTTGACCATTATGTTGTGATGTAGCATGAGATTTTACTTCTTAATAGTAAGGGTCTGGTTCCAGACCAATGTTTATTAACGCTTCAACTTTGATTGGACCTGTAGTTTCAGGCGCCATTACTCTATATCCATatga                                                                                  <  2:383932/139‑1 (MQ=255)
                                                     tAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATatgat                                                                                 >  2:143567/1‑141 (MQ=255)
                                                             aCCCTTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATAcc                                                                         <  1:512612/141‑1 (MQ=255)
                                                              ccATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCa                                                                        <  1:4032/141‑1 (MQ=255)
                                                                 ttATGTTCAAATGTATCATGAGATTTTACTTCTTGATATATAGGACATGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGa                                                                     >  2:373541/1‑141 (MQ=255)
                                                                         agaTGTATTATGAGATTTTAATTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTaa                                                             <  1:397729/139‑1 (MQ=255)
                                                                            tGTTTCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATg                                                          <  2:578098/141‑1 (MQ=255)
                                                                            tGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATTTGATGGTATACCATGATTAAGTATATg                                                          <  2:260730/141‑1 (MQ=255)
                                                                              tATCATGAGATTTTACTTCTTGCTATTTAGGACCTGGTTCAAGACCAATGTTTTTGAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATgcg                                                        >  2:664641/1‑141 (MQ=255)
                                                                              tATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATgcg                                                        >  2:703315/1‑141 (MQ=255)
                                                                               aTCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATgcgc                                                       >  1:778913/1‑141 (MQ=255)
                                                                                  aTGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCctc                                                    >  1:159375/1‑141 (MQ=255)
                                                                                       aTTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAg                                               <  1:137733/141‑1 (MQ=255)
                                                                                       aTTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAg                                               <  1:231696/141‑1 (MQ=255)
                                                                                       aTTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAg                                               <  1:330940/141‑1 (MQ=255)
                                                                                        ttttACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAATGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGt                                              <  1:120932/141‑1 (MQ=255)
                                                                                           tACATCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTATATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTaaacaaaa                                           >  1:654629/1‑134 (MQ=255)
                                                                                            aCTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCCGCACCTTGGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTaaaa                                          <  1:121177/141‑1 (MQ=255)
                                                                                               tcttGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGAACGTGCCAGTGGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACca                                       <  1:662726/141‑1 (MQ=255)
                                                                                                tttGATATTTAGGACCTTGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTTTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACcat                                      <  2:146389/140‑1 (MQ=255)
                                                                                                tttGATATGTAGGAGCTTGTTCAACGCCATTGTCTGTTAACGCTTCAACTTTTATTCTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGACTAAGTAAATGCGCCTCTACAGTAAAACcat                                      <  2:829280/140‑1 (MQ=255)
                                                                                                cttGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACcat                                      <  2:332383/141‑1 (MQ=255)
                                                                                                     tatTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCatga                                 <  1:69171/141‑1 (MQ=255)
                                                                                                       tttAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGCTGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCatgatg                               >  1:895073/1‑141 (MQ=255)
                                                                                                          aGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGata                            >  2:708661/1‑141 (MQ=255)
                                                                                                          aGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGata                            >  1:748135/1‑141 (MQ=255)
                                                                                                          aGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACAGTGAGTATGGCTGGAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGata                            >  1:436911/1‑141 (MQ=255)
                                                                                                                tGGTTCAAGACCAATGTTTTTTTTTTTTTTTACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTTTATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATTATGATATTTCaa                      >  2:471876/1‑140 (MQ=255)
                                                                                                                  gTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGCTTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAgat                    <  1:861326/141‑1 (MQ=255)
                                                                                                                   ttCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGTGCCTCTACAGTAAAACCATCATGATGATATGTCAgatg                   <  1:824932/141‑1 (MQ=255)
                                                                                                                      aaGCCCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGatc                <  2:12157/141‑1 (MQ=255)
                                                                                                                        gACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGatcat              <  1:275131/141‑1 (MQ=255)
                                                                                                                         gccAATGTTTTTTAACGCTTCAACATGGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCAGCATGATGATATGTCAGATGatcatc             <  2:736275/140‑1 (MQ=255)
                                                                                                                         aCCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGatcatc             <  1:143567/141‑1 (MQ=255)
                                                                                                                         aCCAATGCTTTTTAACGCTTAAACATTGATTGTACCTGTAGTTTCAGGGTCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGatcctc             <  2:298960/141‑1 (MQ=255)
                                                                                                                              tGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGATCATCGAttt        <  1:543406/141‑1 (MQ=255)
                                                                                                                                 tttttAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGATCATCGATTTCaa     <  1:184096/141‑1 (MQ=255)
                                                                                                                                  ttttAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAATAAAACCATCATGAAGATAAGTCAGATGATCATCAATTTCaaa    >  2:375885/1‑140 (MQ=255)
                                                                                                                                    ttAACGCTTCGACATTGACTGTACCTGTAGTTTCACGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGATCATCGATTTCAatat  <  2:106228/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
AAATTGCTACAACTGGACCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGATCATCGATTTCAATAT  >  CP000730/1619998‑1620270

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