Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 999,571 A→G Q378R (CAG→CGG)  ycbS → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913999,5710AG99.6% 782.3 / ‑4.6 240Q378R (CAG→CGG) ycbSputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (101/137);  total (101/139)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.12e-01

Only 100 of 254 total aligned reads displayed.

CGTTCTCCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCAGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACC  >  NC_000913/999481‑999660
                                                                                          |                                                                                         
 gTTCTCCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGga                                                                                        <  2:2305383/93‑1 (MQ=255)
 gTTCTCCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGga                                                                                        >  2:1439641/1‑93 (MQ=255)
  ttCTCCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGgag                                                                                       >  2:3669258/1‑93 (MQ=255)
   tctcCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGgaga                                                                                      <  2:4297757/93‑1 (MQ=255)
     tcCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGc                                                                                    <  1:688376/93‑1 (MQ=255)
       cAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCGATAGTAATGAACAGCGGGAGGGCaa                                                                                  >  1:823663/1‑93 (MQ=255)
       cAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCaa                                                                                  >  1:294657/1‑93 (MQ=255)
        aGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCaaa                                                                                 <  2:2891404/93‑1 (MQ=255)
         gtgtGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAAt                                                                                >  1:6055285/1‑93 (MQ=255)
          tgtgCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAAtt                                                                               <  2:1928610/93‑1 (MQ=255)
           gtgCCATTACTCCAGCGTCAGGGGCGTATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAAttt                                                                              <  1:1426962/93‑1 (MQ=255)
           gtgCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAAttt                                                                              <  2:33762/93‑1 (MQ=255)
           gtgCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAAttt                                                                              <  1:5228628/93‑1 (MQ=255)
            tgCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACTGAACCAATAGTAATTAACAGCGGGAGAGCATATTTg                                                                             >  1:6621192/1‑93 (MQ=255)
             gCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGc                                                                            >  2:300196/1‑93 (MQ=255)
             gCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGc                                                                            <  1:3120962/93‑1 (MQ=255)
              ccATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGcc                                                                           >  2:2175725/1‑93 (MQ=255)
               cATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGccc                                                                          <  2:1201430/93‑1 (MQ=255)
               cATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGccc                                                                          <  2:3943525/93‑1 (MQ=255)
                aTTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCa                                                                         <  2:6623683/93‑1 (MQ=255)
                aTTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCa                                                                         <  1:5724/93‑1 (MQ=255)
                aTTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCa                                                                         <  2:5524951/93‑1 (MQ=255)
                 ttACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAg                                                                        <  2:3896600/93‑1 (MQ=255)
                   aCTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGc                                                                      <  2:4609501/93‑1 (MQ=255)
                   aCTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGc                                                                      >  2:447142/1‑93 (MQ=255)
                     tCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCa                                                                    >  2:4101153/1‑93 (MQ=255)
                     tCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCa                                                                    <  1:5000330/93‑1 (MQ=255)
                        aGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGt                                                                 <  1:4577788/93‑1 (MQ=255)
                           gTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGc                                                              <  2:3299442/93‑1 (MQ=255)
                                gggCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTgg                                                         <  1:1608887/93‑1 (MQ=255)
                                gggCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTgg                                                         <  2:401889/93‑1 (MQ=255)
                                 ggCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTggg                                                        <  1:1400840/93‑1 (MQ=255)
                                 ggCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTggg                                                        >  1:4036961/1‑93 (MQ=255)
                                  gCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTgggg                                                       <  1:5098132/93‑1 (MQ=255)
                                  gCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTgggg                                                       <  2:142466/93‑1 (MQ=255)
                                    gAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCg                                                     <  1:4701302/93‑1 (MQ=255)
                                    gAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCg                                                     >  1:4242355/1‑93 (MQ=255)
                                    gAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCg                                                     <  1:5813133/93‑1 (MQ=255)
                                     aaTCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCgg                                                    >  1:3704468/1‑93 (MQ=255)
                                     aaTCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCgg                                                    >  1:5578037/1‑93 (MQ=255)
                                      aTCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGa                                                   <  2:2633307/93‑1 (MQ=255)
                                      aTCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGa                                                   >  1:2700252/1‑93 (MQ=255)
                                      aTCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGa                                                   >  1:6080227/1‑93 (MQ=255)
                                         aaaTACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCg                                                >  1:6356512/1‑93 (MQ=255)
                                         aaaTACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCg                                                >  1:1843229/1‑93 (MQ=255)
                                         aaaTACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCg                                                <  2:622626/93‑1 (MQ=255)
                                         aaaTACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCg                                                <  2:2012169/93‑1 (MQ=255)
                                           aTACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTg                                              >  1:1780489/1‑93 (MQ=255)
                                           aTACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTg                                              <  2:2070364/93‑1 (MQ=255)
                                            tACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTgg                                             <  2:4350333/93‑1 (MQ=255)
                                            tACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTgg                                             <  2:6086266/93‑1 (MQ=255)
                                            tACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTgg                                             <  1:99983/93‑1 (MQ=255)
                                            tACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTgg                                             >  1:2009116/1‑93 (MQ=255)
                                             aCGCGGTGACACTGGCGAACTACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTggg                                            >  1:6560848/1‑93 (MQ=255)
                                               gcgGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGa                                          <  1:2479721/93‑1 (MQ=255)
                                                cgGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGaa                                         >  2:3142869/1‑93 (MQ=255)
                                                   tGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGa                                      >  1:4537238/1‑93 (MQ=255)
                                                    gACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGAc                                     >  1:2713780/1‑93 (MQ=255)
                                                    gACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGAc                                     <  2:3759459/93‑1 (MQ=255)
                                                        cTGCCGAAATACAGAACCAATAGTAATGAGCAGCGGGAGAGCAGATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTgg                                 <  1:3653792/93‑1 (MQ=255)
                                                          ggCGAAATACAGAACCAATAGTAATGAACAGCTGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTa                               <  1:3270123/93‑1 (MQ=255)
                                                          ggCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTa                               <  2:6721464/93‑1 (MQ=255)
                                                          ggCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTa                               >  2:802914/1‑93 (MQ=255)
                                                          ggCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTa                               <  2:691454/93‑1 (MQ=255)
                                                           gCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAt                              >  1:3987035/1‑93 (MQ=255)
                                                           gCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAt                              <  2:1325275/93‑1 (MQ=255)
                                                           gCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAt                              >  1:2150058/1‑93 (MQ=255)
                                                             gAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtgg                            >  1:1348487/1‑93 (MQ=255)
                                                             gAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtgg                            >  2:5029051/1‑93 (MQ=255)
                                                             gAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtgg                            <  1:6207552/93‑1 (MQ=255)
                                                             gAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtgg                            >  2:3443022/1‑93 (MQ=255)
                                                             gAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGAGGACCGTGGGGAACGACCTGGTAtgg                            >  2:5254762/1‑93 (MQ=255)
                                                              aaaTACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGGACGACCTGGTAtggt                           >  1:6504924/1‑93 (MQ=255)
                                                              aaaTACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggt                           >  2:4192693/1‑93 (MQ=255)
                                                               aaTACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggtg                          >  2:2459142/1‑93 (MQ=255)
                                                                  aCAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggtggtg                       <  2:324996/93‑1 (MQ=255)
                                                                   cAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggtggtgg                      <  2:4846043/93‑1 (MQ=255)
                                                                   cAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggtggtgg                      <  1:4109246/93‑1 (MQ=255)
                                                                        ccAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAat                 >  2:1692720/1‑93 (MQ=255)
                                                                         cAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAata                >  2:5772612/1‑93 (MQ=255)
                                                                         cAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAata                >  1:2556019/1‑93 (MQ=255)
                                                                         cAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAata                <  1:1031473/93‑1 (MQ=255)
                                                                          aaTAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGTGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAatat               >  1:5312895/1‑93 (MQ=255)
                                                                          aaTAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAatat               <  1:496451/93‑1 (MQ=255)
                                                                          aaTAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAatat               <  1:5703617/93‑1 (MQ=255)
                                                                           aTAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATg              <  2:290620/93‑1 (MQ=255)
                                                                           aTAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATg              <  1:3370129/93‑1 (MQ=255)
                                                                           aTAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATg              >  2:3839541/1‑93 (MQ=255)
                                                                            tAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGc             >  2:1069712/1‑93 (MQ=255)
                                                                            tAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGc             <  1:2439887/93‑1 (MQ=255)
                                                                             aGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCt            <  1:5992738/93‑1 (MQ=255)
                                                                             aGTAATGAACAGCGAGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCt            <  1:1153642/93‑1 (MQ=255)
                                                                              gTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTg           >  1:4389346/1‑93 (MQ=255)
                                                                                aaTGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGaa         <  1:3081745/93‑1 (MQ=255)
                                                                                aaTGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGaa         <  2:4979797/93‑1 (MQ=255)
                                                                                aaTGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGGATGGTGGTGGACAATATGCTGaa         >  2:4179187/1‑93 (MQ=255)
                                                                                 aTGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAat        <  1:4798812/93‑1 (MQ=255)
                                                                                  tGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAata       <  1:5124644/93‑1 (MQ=255)
                                                                                    aaCAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATAtt     <  1:4648959/93‑1 (MQ=255)
                                                                                       agcggGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTAcc  >  2:1822224/1‑93 (MQ=255)
                                                                                          |                                                                                         
CGTTCTCCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCAGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACC  >  NC_000913/999481‑999660

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: