Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,720,240 A→T 6.0% N122I (AAT→ATT)  slyB → outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,720,2400AT 6.0% 288.2 / 7.6 83N122I (AAT→ATT) slyBouter membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (29/49);  new base T (2/3);  total (31/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.73e-01

TCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAT  >  NC_000913/1720093‑1720388
                                                                                                                                                   |                                                                                                                                                    
tCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTaa                                                                                                                                                      <  1:786521/148‑1 (MQ=255)
tCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTaa                                                                                                                                                      <  1:380876/148‑1 (MQ=255)
 ccTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATa                                                                                                                                                    <  3:863995/149‑1 (MQ=255)
   tGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATAcc                                                                                                                                                  <  3:126683/149‑1 (MQ=255)
    ggggAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATAcc                                                                                                                                                  >  1:857134/1‑148 (MQ=255)
      ggAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTATTACcat                                                                                                                                                <  2:74789/148‑1 (MQ=255)
      ggAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACcat                                                                                                                                                >  1:400271/1‑148 (MQ=255)
         aTACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCAtg                                                                                                                                            <  3:474022/149‑1 (MQ=255)
          tACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCAtg                                                                                                                                            <  2:837450/148‑1 (MQ=255)
           aCTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCAtgg                                                                                                                                           <  2:31649/148‑1 (MQ=255)
            cTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCAtggt                                                                                                                                          <  1:289142/148‑1 (MQ=255)
              gTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTATTACcatcgtc                                                                                                                                            >  2:858178/1‑141 (MQ=255)
              gTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTATTACcatcgtc                                                                                                                                            <  1:858178/144‑4 (MQ=255)
               ttGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATggtggt                                                                                                                                       <  1:832864/148‑1 (MQ=255)
               ttGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTa                                                                                                                                      <  3:202088/149‑1 (MQ=255)
                tggtggCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTa                                                                                                                                      <  1:364219/148‑1 (MQ=255)
                  gtggCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACa                                                                                                                                    >  2:694198/1‑148 (MQ=255)
                    ggCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGaa                                                                                                                                 <  3:643634/149‑1 (MQ=255)
                        gAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACa                                                                                                                              >  2:380242/1‑148 (MQ=255)
                        gAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACa                                                                                                                              <  1:694198/148‑1 (MQ=255)
                         aaCCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAg                                                                                                                            <  4:235361/149‑1 (MQ=255)
                              ggCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCaa                                                                                                                        >  1:737620/1‑148 (MQ=255)
                              ggCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCaa                                                                                                                        <  1:194169/148‑1 (MQ=255)
                                  ttCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACt                                                                                                                    >  2:768019/1‑148 (MQ=255)
                                  ttCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACt                                                                                                                    <  2:33316/148‑1 (MQ=255)
                                  ttCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTc                                                                                                                   <  4:754254/149‑1 (MQ=255)
                                    ctctGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCaa                                                                                                                        >  4:487453/1‑142 (MQ=255)
                                    ctctGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCaa                                                                                                                        <  3:487453/142‑1 (MQ=255)
                                    ctctGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGt                                                                                                                 <  3:362952/149‑1 (MQ=255)
                                      ctGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGtt                                                                                                                <  1:65127/148‑1 (MQ=255)
                                         gCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTtct                                                                                                             <  2:355270/148‑1 (MQ=255)
                                            aCTGCAGCAGGCGCTGTGGCAGGTGGCGTAGCTGGTCAGAGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTtctctc                                                                                                          >  1:265104/1‑148 (MQ=255)
                                                      gcgcTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACa                                                                                                                                    <  4:776190/112‑1 (MQ=255)
                                                      gcgcTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACgt                                                                                               <  4:179915/149‑1 (MQ=255)
                                                      gcgcTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACg                                                                                                <  1:341804/148‑1 (MQ=255)
                                                      gcgcTGTTGCAGGTGGCGTAGCTGGTCAGCGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACa                                                                                                                                    >  3:776190/1‑112 (MQ=255)
                                                                  gTGGCGTAGCTGGTCTGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGGACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGc                                                                                    <  2:265104/148‑1 (MQ=255)
                                                                            tGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGc                                                                          <  1:768019/148‑1 (MQ=255)
                                                                            tGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGc                                                                          <  2:574798/148‑1 (MQ=255)
                                                                            tGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCa                                                                         <  4:429102/149‑1 (MQ=255)
                                                                               tCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCg                                                                                                                  >  2:441130/1‑105 (MQ=255)
                                                                               tCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCg                                                                                                                  <  1:441130/105‑1 (MQ=255)
                                                                                  gggCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAg                                                                    <  2:90758/148‑1 (MQ=255)
                                                                                   ggCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGt                                                                  >  4:774400/1‑149 (MQ=255)
                                                                                      gTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGa                                                                <  2:857134/148‑1 (MQ=255)
                                                                                      gTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGa                                                                <  1:183367/148‑1 (MQ=255)
                                                                                      gTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGAc                                                               >  3:388114/1‑149 (MQ=255)
                                                                                       tACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGAcc                                                              >  3:368085/1‑149 (MQ=255)
                                                                                        aCAGATTGCAATGAACAAAACGCAGGGAGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGAcc                                                              >  2:647212/1‑148 (MQ=255)
                                                                                           gagTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGttt                                                          <  4:388114/149‑1 (MQ=255)
                                                                                             gTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTtc                                                         <  1:380242/148‑1 (MQ=255)
                                                                                                 aaTGAACGAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCTGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCcgc                                                    >  4:766613/1‑149 (MQ=255)
                                                                                                 aaTGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCcg                                                     >  2:438173/1‑148 (MQ=255)
                                                                                                  aTGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACa                                                                                                                                    >  1:345317/1‑68 (MQ=255)
                                                                                                  aTGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACa                                                                                                                                    <  2:345317/68‑1 (MQ=255)
                                                                                                  aTGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCcgcg                                                   >  4:511044/1‑149 (MQ=255)
                                                                                                  aTGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCcgc                                                    >  2:724739/1‑148 (MQ=255)
                                                                                                       cAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTc                                                                      >  1:848487/1‑125 (MQ=255)
                                                                                                       cAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTc                                                                      <  2:848487/125‑1 (MQ=255)
                                                                                                       cAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTaa                                               >  1:131814/1‑148 (MQ=255)
                                                                                                        aaaaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATg                                                                            >  4:781919/1‑118 (MQ=255)
                                                                                                        aaaaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATg                                                                            <  3:781919/118‑1 (MQ=255)
                                                                                                         aaaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAAt                                             >  1:206762/1‑148 (MQ=255)
                                                                                                          aaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAAt                                                                             >  2:645835/1‑115 (MQ=255)
                                                                                                          aaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATa                                            <  1:847640/148‑1 (MQ=255)
                                                                                                          aaCGCAGGGTGTCGAGATGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAAt                                                                             <  1:645835/115‑1 (MQ=255)
                                                                                                           aCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCgg                                                                                                       >  3:46599/1‑88 (MQ=255)
                                                                                                           aCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCgg                                                                                                       <  4:46599/88‑1 (MQ=255)
                                                                                                           aCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAg                                          <  3:542826/149‑1 (MQ=255)
                                                                                                             gCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACGCTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGtt                                        <  4:651309/149‑1 (MQ=255)
                                                                                                              cAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGttt                                       <  3:774400/149‑1 (MQ=255)
                                                                                                              cAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGtt                                        >  2:6962/1‑148 (MQ=255)
                                                                                                                      cGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGtgtg                                <  2:737620/148‑1 (MQ=255)
                                                                                                                            ggAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCgtgt                                 >  2:544949/1‑141 (MQ=255)
                                                                                                                            ggAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCgtgt                                 <  1:544949/141‑1 (MQ=255)
                                                                                                                            ggAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAgg                          <  2:400271/148‑1 (MQ=255)
                                                                                                                              aaaTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGtt                        >  1:246137/1‑148 (MQ=255)
                                                                                                                                aTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTAc                      <  1:438173/148‑1 (MQ=255)
                                                                                                                                aTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACt                     <  4:144747/149‑1 (MQ=255)
                                                                                                                                aTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACt                     <  3:511044/149‑1 (MQ=255)
                                                                                                                                     tAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGGCAGCAATGGCAGTCAGGTGACCGTTGCTCCGCGCTAAATATGTTTGCGTGTGGTCAGGTTACTGAcc                 <  1:647212/148‑1 (MQ=255)
                                                                                                                                     tAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACca                <  4:368085/149‑1 (MQ=255)
                                                                                                                                         catgATGGTATTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACcaca              >  4:307052/4‑147 (MQ=255)
                                                                                                                                         catgATGGTATTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACcaca              <  3:307052/144‑1 (MQ=255)
                                                                                                                                             aTGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACTACACGccc         >  2:202310/1‑148 (MQ=255)
                                                                                                                                                   aTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAt  >  4:787583/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAT  >  NC_000913/1720093‑1720388

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: