Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,071,598 A→G R680G (AGG→GGG)  b2000 → outer membrane fluffing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,071,5980AG96.5% 343.0 / 1.1 114R680G (AGG→GGG) b2000outer membrane fluffing protein
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base G (52/58);  total (53/61)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 118 total aligned reads displayed.

GGGGCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAA  >  NC_000913/2071510‑2071687
                                                                                        |                                                                                         
ggggCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCt                                                                                       <  1:4455686/93‑1 (MQ=255)
  ggCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGc                                                                                     >  2:3478660/1‑93 (MQ=255)
  ggCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGc                                                                                     <  2:1715133/93‑1 (MQ=255)
   gCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGACGGGGCCTTTGTCCAGGGGAACGGGCTGCa                                                                                    <  2:2907550/93‑1 (MQ=255)
   gCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCa                                                                                    <  1:3659374/93‑1 (MQ=255)
   gCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCa                                                                                    <  1:1722406/93‑1 (MQ=255)
   gCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGAAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCa                                                                                    <  2:1425041/93‑1 (MQ=255)
     tGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAgg                                                                                  >  1:1957620/1‑93 (MQ=255)
      ggCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGc                                                                                 >  2:3119369/1‑93 (MQ=255)
      ggCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGc                                                                                 <  2:2160489/93‑1 (MQ=255)
       gCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGcc                                                                                <  1:2421143/93‑1 (MQ=255)
          aCCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGt                                                                             >  2:1631995/1‑93 (MQ=255)
            cAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGc                                                                           >  1:1062561/1‑93 (MQ=255)
            cAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGc                                                                           >  1:1062711/1‑93 (MQ=255)
             aGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGcc                                                                          <  2:1421148/93‑1 (MQ=255)
               cGGTAATGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCtt                                                                        <  1:979339/93‑1 (MQ=255)
                 gTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTa                                                                      >  2:1242042/1‑93 (MQ=255)
                 gTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTa                                                                      <  1:4732658/93‑1 (MQ=255)
                 gTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTa                                                                      >  2:2683136/1‑93 (MQ=255)
                  tAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTaa                                                                     >  1:4360132/1‑93 (MQ=255)
                   aaGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAac                                                                    >  2:3587995/1‑93 (MQ=255)
                    aGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAact                                                                   >  1:62680/1‑93 (MQ=255)
                     gggTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAacta                                                                  <  1:171505/93‑1 (MQ=255)
                        tATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTc                                                               >  2:580723/1‑93 (MQ=255)
                         aTTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTcc                                                              >  1:1712397/1‑93 (MQ=255)
                            cAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTc                                                           >  1:430877/1‑93 (MQ=255)
                              ggtggtGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCaa                                                         >  1:1974819/1‑93 (MQ=255)
                                tggtggAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAAcc                                                       <  2:2606642/93‑1 (MQ=255)
                                     gAAGCCATTAACGGTGCCACCACGGAGGAAGGGCCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGAc                                                  <  1:1190258/93‑1 (MQ=255)
                                      aaGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACa                                                 <  1:3871977/93‑1 (MQ=255)
                                       aGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAg                                                >  1:2772676/1‑93 (MQ=255)
                                       aGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAg                                                >  1:1745000/1‑93 (MQ=255)
                                       aGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAg                                                >  2:2195112/1‑93 (MQ=255)
                                        gCCATTAACGGTGCCAACACGGAGGAAGGGGCCTTTGTACAAGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGt                                               <  2:4204322/93‑1 (MQ=255)
                                         ccATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAGCGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGtg                                              >  1:1946813/1‑93 (MQ=255)
                                          cATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCGCTCAACCGGGACAGtgg                                             >  2:2060611/1‑92 (MQ=255)
                                          cATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGtga                                             <  1:4438396/93‑1 (MQ=255)
                                           aTTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGtgat                                            <  1:2859056/93‑1 (MQ=255)
                                           aTTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGtgat                                            <  1:2179427/93‑1 (MQ=255)
                                           aTTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGtgat                                            >  2:2864162/1‑93 (MQ=255)
                                           aTTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGtgat                                            >  1:2992028/1‑93 (MQ=255)
                                            ttAACGGTGCCCCCACGGAGGAAGGGGCCTTTGTCCAGGGGACCGGGCTGCCGGCCGGTGCCTTTACCCACTCCCTCAACCGGGACAGtgatg                                           >  1:3388595/1‑93 (MQ=255)
                                            ttAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGtgatg                                           <  1:4137678/93‑1 (MQ=255)
                                            ttAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGtgatg                                           <  2:4745764/93‑1 (MQ=255)
                                            ttAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGcgatg                                           >  1:3388596/1‑93 (MQ=255)
                                             tcaCGGTGCCCACACGGAGGAAGGGGCATTTGTCCAGCGGAACAGGCTGCAGGCAGGTGCCTTTAACTACTCCCTCAACCGGGACAGtgatga                                          <  1:4514695/91‑1 (MQ=255)
                                             tAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGtgatgg                                          >  1:2628603/1‑92 (MQ=255)
                                             tAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGtgatga                                          <  1:4179450/93‑1 (MQ=255)
                                             tAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGtgatga                                          <  2:2901589/93‑1 (MQ=255)
                                              aaCGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATgag                                         >  2:2094714/1‑93 (MQ=255)
                                              aaCGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATgag                                         >  1:1896893/1‑93 (MQ=255)
                                              aaCGGGGCCACCACGGAGGAAGTGGCCTTGGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATgag                                         <  2:912208/93‑1 (MQ=255)
                                                cGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGagag                                       >  2:1881554/1‑93 (MQ=255)
                                                 ggTGCCACCACGGAGGAAGGGGCCATTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGc                                      <  1:2865079/93‑1 (MQ=255)
                                                  gTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCt                                     <  1:4103412/93‑1 (MQ=255)
                                                  gTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCt                                     <  2:627487/93‑1 (MQ=255)
                                                  gTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCt                                     <  1:499192/93‑1 (MQ=255)
                                                   tGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTg                                    <  1:1311684/93‑1 (MQ=255)
                                                     ccaccaCGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGt                                  <  1:4121765/93‑1 (MQ=255)
                                                        ccacGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATc                               >  2:3737975/1‑93 (MQ=255)
                                                        ccacGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATc                               <  2:3593557/93‑1 (MQ=255)
                                                        ccacGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATc                               <  1:352732/93‑1 (MQ=255)
                                                        ccacGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATc                               <  1:3592690/93‑1 (MQ=255)
                                                          acGGAGGAAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTg                             >  1:3057205/1‑93 (MQ=255)
                                                            ggaggaAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTgcg                           <  2:751130/93‑1 (MQ=255)
                                                              aggaAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCa                         >  2:2274111/1‑93 (MQ=255)
                                                               ggaAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAg                        >  2:1379502/1‑93 (MQ=255)
                                                               ggaAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAg                        >  2:3228066/1‑93 (MQ=255)
                                                               ggaAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAg                        <  1:4227092/93‑1 (MQ=255)
                                                                gaAGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGt                       <  1:4149244/93‑1 (MQ=255)
                                                                  aGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGa                     >  1:3249067/1‑93 (MQ=255)
                                                                  aGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGa                     >  2:3978028/1‑93 (MQ=255)
                                                                  aGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGa                     >  2:3108480/1‑93 (MQ=255)
                                                                  aGGGGCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGa                     >  1:2404545/1‑93 (MQ=255)
                                                                    gtgCCTGTGTCCAGGGGAACGGGCTGCAGTCCGGGGCCTTTAACGACTCCCTCAGCCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGaaa                   <  1:2995233/91‑1 (MQ=255)
                                                                     ggCCTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGaaaa                  <  2:1606940/93‑1 (MQ=255)
                                                                     ggCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGaaaa                  <  1:3862273/93‑1 (MQ=255)
                                                                       ccTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATg                >  2:1514568/1‑93 (MQ=255)
                                                                       ccTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATg                >  2:1514566/1‑93 (MQ=255)
                                                                       ccTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATg                >  2:1403515/1‑93 (MQ=255)
                                                                        cTTTGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGc               <  2:654752/93‑1 (MQ=255)
                                                                         tttGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCt              <  2:4025209/93‑1 (MQ=255)
                                                                          ttGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCtt             <  1:2789089/93‑1 (MQ=255)
                                                                           tGTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTa            <  2:2514318/93‑1 (MQ=255)
                                                                            gTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTAt           <  2:1235254/93‑1 (MQ=255)
                                                                            gTCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTAt           <  2:4005101/93‑1 (MQ=255)
                                                                             tCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATc          <  1:2365711/93‑1 (MQ=255)
                                                                             tCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATc          <  1:2020730/93‑1 (MQ=255)
                                                                             tCCAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATc          <  1:3904955/93‑1 (MQ=255)
                                                                              ccAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCg         >  2:4679065/1‑93 (MQ=255)
                                                                              ccAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCg         >  2:4679066/1‑93 (MQ=255)
                                                                               cAGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGt        <  1:449902/93‑1 (MQ=255)
                                                                                aGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTg       >  2:2197289/1‑93 (MQ=255)
                                                                                aGGGGAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTg       >  2:2149694/1‑93 (MQ=255)
                                                                                 ggggAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGc      >  2:3846124/1‑93 (MQ=255)
                                                                                 ggggAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGc      <  1:1306421/93‑1 (MQ=255)
                                                                                 ggggAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGc      >  2:1615260/1‑93 (MQ=255)
                                                                                 ggggAACGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGc      <  1:1725350/93‑1 (MQ=255)
                                                                                     aaCGGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGaa  <  1:1246442/93‑1 (MQ=255)
                                                                                     aaCAGACTAAAGGCCGGTGCCTTTAACTCCTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGaa  <  1:3323931/93‑1 (MQ=255)
                                                                                        |                                                                                         
GGGGCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAA  >  NC_000913/2071510‑2071687

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